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Nn  number  of  points  making  up  streambed,  streambank,  or  floodplain 

Np  number  of  profiles  for  which  output  is  requested 

Nrs  number  of  segments  comprising  a rating  curve 

Ns  number  of  storm  events  for  which  output  is  requested 

Nse  number  of  segments  comprising  the  walls  of  a generic  structure 

Nt  time-step  reduction  factor  for  diffusion  equation 

Nts  number  of  time  series  for  which  output  is  requested 

Nx  number  of  cross  sections  in  a subreach 


Q discharge  m3/s 

R hydraulic  radius  m 

Re  equivalent  hydraulic  radius  m 

S shear  force  at  the  base  of  the  failure  block  per  unit  channel  length  N/m 
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Uppercase  variables  (continued) 


variable 

meaning 

unit 

sb 

bed  slope 

m/m 

Sc 

slope  of  culvert  barrel 

m/m 

Sf 

friction  slope 

m/m 

Su 

sediment  flux  from  subsurface  layers  to  surface  layer 

m2/s 

T 

time-of-rise  of  flood  wave  (Chapter  1 ) 

s 

T 

time  scale  representing  the  rate  of  adjustment  of  sediment  mass  to 
equilibrium  sediment  mass  (Chapter  2) 

s 

uw 

pore-water  force  acting  on  the  base  of  the  failure  block  per  unit 
channel  length 

N/m 

K 

volume  eroded  from  the  streambank  per  unit  channel  length 

m3/m 

w 

actual  rate  of  bank  erosion 

m/s 

w, 

initial  rate  of  bank  erosion 

m/s 

Ws 

weight  of  failure  block  per  unit  channel  length 

N/m 

weight  of  surface  water  on  the  failure  block  per  unit  channel  length 

Lowercase  variables 

N/m 

variable 

meaning 

unit 

a 

thickness  of  surface  or  active  layer 

m 

c 

celerity  of  flood  wave  (Chapter  1 ) 

m/s 

c 

point  sediment  concentration  (Chapter  2) 

ppmw 

c 

cohesion  (Chapter  3) 

Pa 

c ' 

effective  cohesion 

Pa 

c 

kinematic-wave  celerity 

m/s 

ca 

apparent  cohesion 

Pa 

cd 

developed  cohesion 

Pa 

d 

particle  diameter 

m 

d 

depth  below  soil  surface  (page  61) 

m 

d 

mean  sediment  size 

m 
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Lowercase  variables  (continued) 


variable 

meaning 

unit 

cho 

sediment  size  for  which  50%  of  the  sediment  mixture  is  finer 

m 

dc 

critical  sediment-particle  diameter 

m 

e 

erosion-rate  constant 

m/s 

f 

flow  variable  (page  24) 

fk 

ratio  between  sediment  mass  entering  the  stream  corridor  and  the 
sediment-carrying  capacity  of  the  flow  at  the  inlet 

ft 

ratio  between  tension  crack  height  and  thickness  of  soil  layer  with 
negative  earth  pressures 

g 

gravitational  acceleration 

m2/s 

g 

gravitational  acceleration  vector  (0,0,-g)7^ 

m2/s 

h 

flow  depth 

m 

h 

total  hydraulic  head  (page  60) 

m 

hd 

hydraulic  depth 

m 

i 

hydraulic  gradient 

m/m 

j 

space  index 

k 

sediment-size  class 

m 

ratio  between  actual  change  in  bed  elevation  at  the  outlet  and  that 
computed 

n 

Manning  roughness  coefficient 

s/m 1/3 

n 

time  index  (superscript) 

n ' 

Manning  grain-roughness  coefficient 

s/m 1/3 

ne 

equivalent  Manning  roughness  coefficient 

s/m 1/3 

P 

fraction  of  sediment  available  for  transport 

q 

lateral  inflow  of  water  into  the  stream  corridor 

m2/s 

q.s 

lateral  volumetric  inflow  of  sediment  into  the  stream  corridor 

m2/s 

r 

Meyer-Peter  and  Mueller’s  ratio  of  grain  roughness  and  bed 
roughness 

ru 

pore-water  pressure  ratio 

t 

time 

s 

ll 

flow  velocity 

m/s 
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Lowercase  variables  (continued) 


variable 

meaning 

unit 

ua 

pore-air  pressure 

Pa 

uw 

pore-water  pressure 

Pa 

ll  * 

shear  velocity 

m/s 

X 

space  coordinate  in  horizontal  direction 

m 

XP 

horizontal  coordinate  of  intersect  of  phreatic  surface  and  bank  face 

m 

y 

stage 

m 

z 

space  coordinate  in  vertical  direction 

m 

zp 

vertical  coordinate  of  intersect  of  phreatic  surface  and  bank  face 

m 

zt 

thickness  of  layer  with  negative  earth  pressure 

m 

ztc 

height  of  the  tension  crack 

m 

Symbols 

variable 

meaning 

unit 

a 

coefficient  of  power  function  representing  a flow  variable  (page  23) 

a 

angle  of  seepage  force  (Chapter  3) 

O 

a-L 

step-length  constant 

P 

exponent  of  power  function  representing  a flow  variable  (page  23) 

P 

angle  of  failure  plane  (Chapter  3) 

O 

P* 

fractional  content  by  volume  of  a sediment-size  class  in  the 
streambank 

P5 

fractional  content  by  volume  of  a sediment-size  class  in  the  surface 
layer 

P" 

fractional  content  by  volume  of  a sediment-size  class  in  the 
subsurface  layer 

y 

specific  weight  of  water 

N/m3 

y b 

soil  bulk  unit  weight 

N/m3 

Ys 

specific  weight  of  sediment 

N/m3 

A 

increment 

s 

rate  constant 

1/VPa 

xix 


Symbols  (continued) 

variable  meaning  unit 

Sy  erosion  rate  of  cohesive  floes  m/s 

9 temporal  weighting  factor  (Chapter  1 ) 

9 Shields  parameter  (Chapter  2) 

9t.  critical  Shields  parameter 

X porosity  of  streambed 

v kinematic  viscosity  nr/s 

p density  of  water  kg/m3 

ct  contraction  ratio 

0 contraction  ratio  1 - ct 

ct  total-stress  tensor  Pa 

ct„  total  stress  normal  to  the  base  of  the  failure  block  Pa 

1 shear  stress  Pa 

i ' bed  shear  stress  due  to  grain  resistance  Pa 

xb  bed  shear  stress  Pa 

xc  critical  tractive  stress  Pa 

xd  threshold  bed  shear  stress  below  which  deposition  of  cohesive  Pa 

sediments  occurs 

xe  shear  strength  of  bed  material  Pa 

<j)  angle  of  internal  friction  ° 

<j) 1 effective  angle  of  internal  friction  ° 

<\>b  angle  indicating  increase  in  shear  strength  for  an  increase  in  matric 
suction 

X hiding  coefficient 

4*  matric-suction  force  acting  on  slip  surface  per  unit  channel  length  N/m 

\\f  spatial  weighting  factor 

to  particle  fall  velocity  m/s 
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Introduction 


Background 

The  Conservational  Channel  volution  and  Pollutant  ransport  v ystem 
(CONCEPTS)  is  a computer  model  that  simulates  open-channel  hydraulics, 
sediment  transport,  and  channel  morphology.  CONCEPTS,  version  1.0, 
incorporates  unique  advances  on  alluvial  channel  evolution  made  by  scientists 
at  the  U.S.  Department  of  Agriculture,  Agricultural  Research  Service,  National 
Sedimentation  Laboratory  (NSL)  in  Oxford,  Mississippi.  Presently,  there  are 
two  versions  of  CONCEPTS  available  that  consider  different  scales  of  stream 
systems: 

• watershed-scale  stream  network,  and 

• reach-scale  stream  corridor. 

This  document  only  describes  the  technology  and  input  specifications  of  the 
stream-corridor  version  of  CONCEPTS.  You  may  find  information  on  the 
watershed-scale  version  of  CONCEPTS  on  the  NSL’s  website  at  URL: 

http://www.sedlab.olemiss.edu/AGNPS2001/ConceDts/concepts.html 

You  can  download  both  versions  of  CONCEPTS  from  the  above  URL. 

The  basic  components  of  CONCEPTS  are  channel  hydraulics,  channel 
morphology,  and  transport  of  sediments  and  contaminants.  CONCEPTS 
version  1 .0  does  not  include  the  transport  of  contaminants. 

Purpose  and  Capabilities 

The  overall  objective  of  the  ongoing  work  on  CONCEPTS  is  to  develop  an 
integrated  stream-riparian  corridor  model  to  be  used  in  conjunction  with 
watershed-scale  analysis,  for  evaluating  the  long-term  effectiveness  of 
proposed  corridor  rehabilitation  designs.  This  modeling  capability  is  a critical 
need  for  hydraulic  engineers,  fluvial  geomorphologists,  landscape  architects, 
and  stream  ecologists  involved  in  stream  rehabilitation  projects.  These  projects 
are  usually  challenged  by  the  lack  of  integrated  design  tools  to  evaluate  the 
long-term  stability  of  reconstructed  stream  corridors.  This  model  is  an  integral 
component  of  the  suite  of  modeling  tools  incorporated  into  the  AGNPS98 
technology  at  NSL  supported  by  the  U.S.  Department  of  Agriculture,  Natural 
Resources  Conservation  Service  (http.V/www. sedlab.olemiss.edu/ 
AGNPS.html). 

Version  1.0  of  CONCEPTS  simulates  unsteady,  one-dimensional  flow, 
graded-sediment  transport,  and  bank-erosion  processes  in  stream  corridors.  It 
can  predict  the  dynamic  response  of  flow  and  sediment  transport  to  instream 
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hydraulic  structures.  It  computes  channel  evolution  by  tracking  bed  changes 
and  channel  widening.  Bank  erosion  accounts  for  basal  scour  and  mass  wasting 
of  unstable  cohesive  banks.  CONCEPTS  simulates  transport  of  cohesive  and 
cohesionless  sediments,  both  in  suspension  and  on  the  bed,  and  selectively  by 
size  classes.  CONCEPTS  includes  channel-boundary  roughness  varying  along 
a cross  section,  for  example  due  to  varying  vegetation  patterns. 

CONCEPTS  can  be  used  to  evaluate  the  efficiency  of  instream  grade-control 
structures  to  reduce  sediment  yield  and  to  stabilize  streams.  CONCEPTS  can 
evaluate  location  and  sizing  alternatives  of  grade-control  structures. 
CONCEPTS  can  be  used  to  evaluate  the  design  of  specific  stream  corridor 
rehabilitation  measures  used  for  stream  stability  and  habitat  improvement. 

CONCEPTS  version  2.0,  already  under  development,  will  incorporate  the 
simulation  of  riparian  buffers,  vegetated  streambanks,  and  the  onset  of  channel 
meandering  due  to  the  deposition  of  alternate  bars. 

Limits  of  Application 

CONCEPTS  assumes  or  is  limited  to  the  following: 

• gradually-varying,  one-dimensional  flow; 

• four  types  of  instream  hydraulic  structures: 

pipe  and  box  culverts, 

bridge  crossings  with  a trapezoidal  cross  section, 

drop  structures  with  a trapezoidal  cross  section, 

generic  structures  with  a multi-segment  trapezoidal  cross  section. 

• straight  channels  or  channels  of  very  low  sinuosity; 

• different  hydraulic  roughness  for  bed,  banks,  and  floodplains; 

• total  load  sediment  transport; 

• homogeneous  cohesionless  or  cohesive  bed-material  in  transverse  direction; 

• 13  pre-defmed  sediment  particle-size  classes; 

• non-equilibrium  transport  of  cohesionless  bed-material; 

• linear  erosion  rate  of  cohesive  bed  material; 

• consolidation  of  deposited  cohesive  sediments  is  unavailable; 

• homogeneous  cohesive  bank  material; 

• linear  lateral-erosion  rate  of  the  bank  toe  based  on  excess  shear  stress; 

• slab  or  planar  type  of  bank  failure; 

• (quasi-)  hydrostatic  groundwater  pressure  in  streambank; 

• metric  system  of  units. 
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System  Requirements 


To  use  CONCEPTS  you  need  a PC  with  a Microsoft  Windows®  95  or  later 
operating  system.  Hard  disk  and  memory  usage  will  vary  depending  on  the  size 
of  the  problem.  The  CONCEPTS  executable  file  requires  approximately  600 
KB  of  hard  disk  space.  Input  files  for  typical  applications  require  one  to  two 
MB  of  hard  disk  space,  used  mainly  by  the  discharge  data  file.  Depending  on 
the  output  options  specified  by  the  user,  the  combined  hard  disk  usage  of  the 
output  files  may  vary  between  a couple  of  MB  to  a couple  of  GB.  A minimum 
of  32  MB  of  RAM  is  recommended. 

Overview  of  Manual 


This  manual  presents  the  computer  model  CONCEPTS,  specifically: 

• Part  1,  scientific  basis; 

• Part  2,  data  needs;  and 

• Part  3,  sample  applications. 

Chapters  1 through  3 include  the  scientific  basis  and  cover  hydraulics, 
sediment-transport  mechanics,  and  streambank  erosion.  Formulations  of  these 
components  and  their  numerical  implementation  are  presented. 

Chapters  4 and  5 include  the  format  of  the  input  and  output  data  files.  Chapter 
5 also  presents  a method  to  produce  graphs  from  the  output  data  using 
Microsoft®  Excel. 

Chapters  6 through  8 present  applications  of  CONCEPTS  to  channels  in 
northern  Mississippi  and  eastern  Nebraska. 
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Part  1 


Model 

Description 


Flow  Hydraulics 


Theory 

This  chapter  presents  the  flow  hydraulics  submodel  of 
CONCEPTS.  You  will  find  the  scientific  basis  of  the  mathematical 
formulations  of  one-dimensional  flow  in  open  channels  and  at 
hydraulic  structures  used  by  CONCEPTS.  Further,  you  will  find  the 
numerical  implementation  of  the  mathematical  model,  which  will 
help  you  in  selecting  suitable  user-defined  parameters  incorporated 
in  the  mathematical  formulations. 

CONCEPTS  assumes  the  flow  in  stream  systems  to  be  one  dimensional  along 
the  centerline  of  the  channel.  It  neglects  cross-stream  variations  induced  by 
cross-section  geometrical  features  such  as  riffles  and  pointbars  (Figure  1.1a), 
constrictions  or  expansions  (Figure  1.1b),  and  obstructions  such  as  woody 
debris  or  rocks  (Figure  1.1c). 

Various  one-dimensional  flow  routing  methods  exist.  They  are  commonly 
known  as  (Fread,  1996): 

• lumped  flow  routing  or  hydrologic  routing,  and 

• distributed  flow  routing  or  hydraulic  routing  models. 

Lumped  flow  routing  methods  compute  the  flow  as  a function  of  time  at  one 
location  along  the  stream,  e.g.  Muskingum  and  level-pool  routing  methods. 
Distributed  flow  routing  methods  compute  the  flow  as  a function  of  time 
simultaneously  at  a series  of  cross  sections  along  the  stream.  Lumped  flow 
routing  methods  are  computationally  more  efficient,  whereas  distributed  flow 
routing  methods  are  most  accurate.  Lumped  flow  routing  methods  cannot 
predict  backwater  upstream  of  hydraulic  structures  and  channel  confluences. 
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Flow  Hydraulics 


Figure  1.1  Cross-stream  variations  caused  by:  (a)  bars,  (b)  constrictions,  and  (c) 
obstructions. 


- 


(C) 


CONCEPTS  uses  a distributed  flow  routing  method  to  accurately  model  the 
effects  of  instream  hydraulic  structures  on  stream  corridor  rehabilitation 
designs. 
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Governing  Equations 


Theory 


Open-Channel  Flow 

The  Saint  Venant  equations  describe  unsteady,  gradually-varying  one- 
dimensional open-channel  flow.  Cunge  et  al.  (1980)  give  an  extensive 
overview  of  these  equations.  The  Saint  Venant  equations  consist  of  a 
continuity  equation  representing  mass  conservation  of  water  and  a momentum 
equation  representing  the  conservation  of  fluid  momentum. 

The  continuity  equation  reads 

+ ® = 9 <H 

ot  ox 


where  y is  stage,  Q is  discharge,  B is  flow  top  width,  q is  lateral  flow  into  the 
channel  per  unit  length  of  channel  (e.g.,  overland  flow  or  groundwater  return 
flow),  x is  distance  along  the  channel,  and  t is  time.  The  momentum  equation  is 


dt 


+1.(21 

dx  ' A 


+ gA 


& + S 
dx 


= 0 


(1.2) 


where  A is  flow  area,  g is  gravitational  acceleration,  and  5)  is  friction  slope. 
The  friction  slope  is  defined  as 


s =o: 

7 *2 


nQ2 

a2r4  } 


(1.3) 


where  K is  conveyance,  R is  hydraulic  radius,  and  n is  Manning  roughness 
coefficient.  The  cross-sectional  average  flow  velocity  used  hereafter  is  defined 
as  u = Q/A  . Figure  1 .2  shows  a definition  sketch  of  the  above  variables. 

The  variables  stage,  y,  and  discharge,  Q , are  the  dependent  variables.  The 
dependent  variables  are  those  that  are  determined  by  the  solution  method. 
Time,  t,  and  distance,  x,  are  the  independent  variables.  All  other  variables  are 
functions  of  the  independent  or  dependent  variables. 

The  Saint  Venant  equations  are  also  known  as  the  dynamic  wave  model.  If  we 
neglect  the  inertia  terms  in  ( 1 .2),  the  first  and  second  term  on  the  left-hand  side 
of  ( 1 .2),  for  example  when  Froude  number  is  small,  ( 1 .2)  simplifies  to 

o (i.4) 

Equations  (1.1)  and  ( 1 .4)  are  known  as  the  diffusion  wave  model.  The  diffusion 
wave  model  is  computationally  more  efficient  than  the  dynamic  wave  model. 
Also,  omitting  the  nonlinear  convective  acceleration,  the  second  term  in  (1.2), 
results  in  a more  robust  numerical  model.  The  numerical  approximation  of  the 
convective  acceleration  term  may  produce  erroneous  solutions  or  even  fail  to 
calculate  a solution. 

Ponce  et  al.  (1978)  postulated  the  following  applicability  criterion  for  the 
diffusion  wave  model: 
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Flow  Hydraulics 


Figure  1.2  Definition  sketch  of  hydraulic  variables  used  in  the  governing  equations:  a) 
cross-sectional  view  and  b)  longitudinal  view. 


energy 

line 
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reference 

datum 


j,>  lf>  jh 

~ SbNg 


(1.5) 


where  T is  time-of-rise  of  the  flood  wave,  h is  flow  depth,  and  Sb  is  bed  slope. 
The  flow  depth  and  bed  slope  are  reach-averaged  values.  Depending  on 
criterion  (1.5)  CONCEPTS  automatically  switches  between  the  dynamic  and 
diffusion  wave  models. 


Along  a watercourse  you  will  find  two  types  of  boundaries: 

• external  boundaries , the  upstream  and  downstream  extremities  of  the  stream 
corridor;  and 

• internal  boundaries , locations  where  the  flow  is  rapidly  varied  rather  than 
gradually  varied,  hence  the  St.  Venant  equations  are  not  applicable. 

Example  internal  boundaries  are  hydraulic  structures,  rapids,  etc.  The 
following  section  discusses  the  representation  of  hydraulic  structures  in 
CONCEPTS. 
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Hydraulic  Structures 

Structures  can  act  as  controls  or  transitions , changing  local  rating  curves. 
Since  sediment  transport  is  directly  related  to  the  local  velocities,  hydraulic 
structures  will  not  only  affect  discharge  and  stage  hydrographs  but  also 
sediment  movement  in  stream  channels.  For  example,  grade-control  structures 
reduce  channel  grade  by  creating  backwater  and  inducing  upstream  deposition, 
and  downstream  scour. 

Presently,  there  are  four  types  of  structures  included  in  CONCEPTS- 1.0: 

• culvert, 

• bridge  crossing, 

• drop  structure,  and 

• generic  structure. 

The  generic  structure  can  be  any  structure  for  which  a rating  curve  is  available. 

To  guarantee  an  efficient  solution  method,  the  mathematical  representation  of 
the  flow  at  hydraulic  structures  has  to  be  equivalent  to  that  of  open-channel 
flow,  that  is,  equations  (1.1)  and  (1.2).  To  model  the  flow  at  structures 
CONCEPTS  uses  continuity  and  dynamic  equations , which  establish  a relation 
between  discharges  and  stages  upstream  and  downstream  of  the  structure. 

The  continuity  equation 

f = 0 (1.6) 

OX 

states  that  discharge  across  the  structure  is  constant.  Thus,  the  discharge 
upstream  of  the  structure  equals  that  downstream  of  the  structure.  The 
continuity  equation  is  the  same  for  all  types  of  structures. 

The  dynamic  equation  describes  the  flow  at  the  structure  by  relating  the 
discharge  through  the  structure  to  the  stages  upstream  and  downstream  of  the 
structure.  Because  flow  pattern  varies  among  different  types  of  structures,  each 
structure  has  a different  dynamic  equation. 

Culvert 

CONCEPTS- 1.0  models  two  culvert  shapes:  pipe  culverts  and  box  culverts. 
Various  flow  types  characterize  the  flow  through  culverts.  Bodhaine  (1968) 
classifies  culvert  flow  into  six  flow  types  on  the  basis  of  location  of  the  control 
section  and  the  relative  heights  of  the  headwater  and  tailwater  elevations. 
CONCEPTS  distinguishes  three  flow  types  (see  also  Figure  1.3): 

• flow  controlled  by  the  culvert  inlet; 

• flow  controlled  by  both  culvert  inlet  and  outlet;  and 

• uniform  flow  in  the  culvert  barrel. 

The  latter  flow  type  is  a representation  of  culvert  flow  at  base  flow  conditions. 
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Figure  1.3  Flow  types  at  culverts:  inlet  control,  inlet  and  outlet  control,  and  uniform 
flow.  The  numbers  1, 2,  and  3 indicate  sections  used  in  the  structure’s 
dynamic  equation. 
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In  the  case  of  inlet  control  flow,  the  capacity  of  the  culvert  opening  limits  the 
capacity  of  the  culvert.  The  flow  in  the  culvert  entrance  is  critical.  The  US 
Federal  Highway  Administration  (USFHWA,  1985)  has  published 
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nomographs  that  estimate  the  headwater  upstream  of  the  culvert  for  a given 
discharge  at  the  culvert. 

The  discharge  at  the  inlet  (section  2,  referred  to  as  a subscript  in  the  equation 
below)  is  then 

e2  - a2J^  «i.7, 

where  the  flow  area  and  top  width  are  for  critical  depth. 

The  US  Federal  Highway  Administration  (1985)  distinguishes  two  headwater 
equations;  one  in  the  case  where  the  inlet  is  not  submerged  and  another  in  the 
case  where  the  inlet  is  submerged.  When  the  inlet  is  not  submerged  the 
headwater  equation  is 

u2 

H'  = hl  + 2i  + °( C|  Ql ' ~ C}Sc)  (1'8) 

and  when  the  inlet  is  submerged 

H,  = D(CAQl  + C5-C2Sc)  (1.9) 


where  H is  total  head,  u is  velocity,  D is  height  of  the  culvert  barrel,  Sc  is  slope 

of  the  culvert  barrel,  Q = Q/(AcjD)  in  which  Ac  is  the  area  of  the  culvert 
barrel,  and  Cj  through  C5  are  coefficients  depending  on  inlet  and  outlet 
characteristics,  that  is  the  type  of  culvert  (USFHWA,  1985). 

Equations  ( 1 .7)  through  ( 1 .9)  provide  a framework  to  relate  discharge  at  the 
culvert  to  stage  upstream  of  the  culvert  at  section  1 . 


If  the  tailwater  elevation  causes  the  culvert  barrel  to  run  full,  both  culvert  inlet 
and  outlet  control  the  flow  through  the  culvert.  The  dynamic  equation  then 
reads 


Q 2 ^2 


l2g(y  i -y3 ) 
V CL 


(1.10) 


where  CL  is  an  energy-loss  coefficient  that  is  a combination  of  entrance,  exit, 
and  friction  losses.  Table  4.20  shows  entrance  loss  coefficients  given  by 
USFHWA  (1985).  For  a sudden  expansion  of  flow,  such  as  in  a typical  culvert, 
the  exit  loss  coefficient  is  normally  set  to  unity.  In  general,  exit  loss 
coefficients  can  vary  between  0.3  and  1 .0.  The  exit  loss  coefficient  should  be 
reduced  as  the  transition  becomes  less  abrupt.  The  friction  loss  coefficient,  C y, 
is  calculated  as 

2 

Cf  = 2 gLr-^—  (1.11) 

/ S crA/7>  ' 1 


where  Lc  is  the  length  of  the  culvert  barrel,  and  R2  is  the  hydraulic  radius  of  the 
flow  in  the  culvert  barrel. 


13 


Flow  Hydraulics 


The  flow  depth  in  the  culvert  barrel  is  very  small  for  the  case  of  base  flow  in 
the  channel.  To  improve  the  robustness  of  the  solution  algorithm  CONCEPTS 
assumes  uniform  flow  in  the  culvert.  The  Manning  equation  then  gives 
discharge  as 

Qi  = C2R2JSc  (1.12) 

The  shape  of  the  culvert  is  defined  by  the  USFHWA  chart  and  scale  numbers. 
The  chart  and  scale  number  refer  to  a series  of  nomographs  in  USFHWA 
(1985).  Table  4.19  lists  the  chart  and  scale  numbers  for  pipe  and  box  culverts. 


Bridge  Crossing 

CONCEPTS  distinguishes  two,  ‘low-flow  conditions’  flow  types  at  bridge 
crossings.  Low-flow  conditions  exist  when  the  free  surface  of  the  water  is  not 
in  contact  with  the  lower  chord  of  the  bridge  deck.  Henderson  (1966)  terms 
these  flow  types  A and  B (see  Figure  1 .4). 


Class  A Flow  occurs  when  the  flow  through  the  bridge  opening  is  subcritical. 
Class  B Flow  or  “choked”  flow  occurs  when  the  flow  in  the  bridge  opening  is 
critical.  The  transition  from  class  A to  B and  vice  versa  is  a function  of  the 
contraction  ratio  o,  which  is  the  ratio  of  the  width  of  the  bridge  opening  to  the 
total  width  of  the  bridge  crossing.  Equating  specific  momentum  at  section  2 
and  section  3 yields  the  limiting  contraction  ratio 


(2+  l/a)V? 

a = 

0+2  Fj)1 


(1.13) 


where  F is  the  flow  Froude  number. 


From  experiments  on  Class  A flow,  Yamell  (1934a,  1934b)  obtained  the 
following  relation  between  the  change  in  water-surface  elevation  across  the 
bridge  and  the  downstream  Froude  number 

= C F32(C  +5f32-0.6)(a  + 15i4)  (1.14) 

"3 


where  Cp  is  a pier-shape  coefficient,  and  ct  = 1 - a . From  (1.14)  we  can 
derive  a relation  between  discharge  at  the  bridge  and  upstream  and 
downstream  stages 


Qi  A- 


0.\gh3(C _-0.6)+  |0.01(C  - 0.6)  + 


T1-T3 


(1.15) 


A|  5h3Cp(o  + 15a  ) 

Table  4.22  lists  values  of  the  pier  shape  coefficient  for  various  pier  shapes. 


Class  B Flow  is  similar  to  inlet-control  culvert  flow.  Therefore  the  governing 
equations  are  (1.7)  and  the  headwater  equation 
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Figure  1.4  Flow  at  bridge  crossings:  (a)  top  view  of  bridge  crossing  showing  a flow 
contraction  at  section  2;  and  lateral  views  of  (b)  Low  Flow  Class  A and  (c) 
Low  Flow  Class  B.  The  numbers  1,  2,  and  3 indicate  sections  used  in  the 
structure’s  dynamic  equation. 
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(1.16) 


where  CL  is  an  energy-loss  coefficient  depending  on  pier  shape.  For  square- 
ended  piers  C,  = 0.35  and  for  rounded  ends  C,  = 0. 1 8 . The  above  values 

. p ^ p 

of  energy  loss  are  for  a pier  length-width  ratio  of  four;  longer  piers  increase  the 
loss.  Hendersen  (1966)  states  that  the  amount  of  increase  is  approximately  5% 
and  10%  for  7:1  (length: width)  and  13:1  piers  respectively,  with  rounded  ends. 
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Figure  1.5  Flow  at  a drop  structure:  (a)  free  overfall  and  (b)  submerged  flow.  The 
numbers  1,  2,  and  3 indicate  sections  used  in  the  structure’s  dynamic 
equation. 
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For  the  square-ended  piers  the  backwater  was  slightly  less  for  longer  piers  than 
for  the  4: 1 pier. 

Drop  Structure 

CONCEPTS  distinguishes  two  flow  types  at  drop  structures  based  on  the 
degree  of  submergence  of  the  drop  structure  (see  Figure  1 .5).  The  flow  at  the 
brink  of  the  drop  structure  is  a free  (ventilated)  overfall  for  low  tailwater 
elevations,  and  is,  therefore,  critical.  The  drop  structure  becomes  submerged 
for  high  tailwater  elevations  resulting  in  a non-ventilated  overfall  condition. 

In  the  case  of  critical  flow  the  governing  equations  are  (1.7)  and  an  expression 
for  the  headwater 


where  C is  an  energy-loss  coefficient  representing  entrance  losses. 

In  the  case  of  a submerged  drop  structure,  both  upstream  and  downstream  flow 
conditions  control  the  flow  at  the  drop  structure.  CONCEPTS  uses  (1.10)  to 
compute  discharge.  Future  versions  of  CONCEPTS  will  incorporate  an 
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improved  description  of  this  flow  type  based  on  the  work  of  Wu  and 
Rajaratnam  (1996,1998). 

Generic  Structure 

A generic  structure  is  any  structure  for  which  a rating  curve  is  available.  The 
governing  equations  are  the  rating  curve 

q2  = Q(h2)  n-18) 

and  the  headwater  equation 


where  C,  is  an  energy-loss  coefficient  representing  entrance  losses  at  the 

g 

generic  structure. 
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Representation  of  Stream  Cross-Section  and  Stream 
Corridor 


Corridor 

The  stream  corridor  is  schematized  as  reaches  connecting  cross  sections 
(Figure  1 .6).  A reach  is  a stream  segment  that  transfers  information  between 
two  cross  sections.  A cross  section  is  a node  that  holds  hydraulic  information. 
The  information  transferred  between  cross  sections  tells  CONCEPTS  how  the 
flow  at  one  cross  section  relates  to  the  flow  at  its  upstream  and  downstream 
cross  section. 

CONCEPTS  assumes  that  the  channel  is  straight.  A sinuous  stream  therefore 
needs  to  be  mapped  to  a straight  line  (Figure  1 .7).  CONCEPTS  does  not 
require  any  data  to  be  entered  for  reaches.  CONCEPTS  automatically 
determines  reach  profile  and  length  from  river  kilometer  and  thalweg  elevation 
of  its  upstream  and  downstream  cross  sections. 

Cross  Section 

Cross  sections  specify  the  boundary  geometry  of  the  stream  corridor.  Cross 
sections  characterize  the  flow-carrying  capability  of  the  stream  and  its  adjacent 
floodplain.  They  should  extend  across  the  entire  floodplain  and  should  be 
perpendicular  to  the  anticipated  flow  lines. 

Cross  sections  are  required  at  representative  locations  throughout  the  stream 
corridor.  They  must  accurately  represent  the  stream  wise  changes  in  flow,  bed 
slope,  shape,  flow-resistance  characteristics,  and  bank-  and  bed-material 
geotechnical  properties  (i.e.,  cohesion,  angle  of  internal  friction,  soil-water 
properties,  etc.).  Where  abrupt  changes  occur,  several  cross  sections  should  be 
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Figure  1.6  Representation  of  a stream  corridor  as  reaches  and  nodes  (arrows  depict 
information  exchanged  among  adjacent  nodes). 


used  to  describe  the  change.  Fread  (1988)  gives  the  following  maximum 
computational  space  step  or  distance  between  cross  sections.  Ax: 

Ax  < — (1.20) 

20 

where  c is  celerity  of  the  flood  wave  propagating  through  the  stream  corridor. 
Basco  (1987)  states  that  the  ratio  of  flow  area  at  successive  cross  sections 
should  not  be  smaller  than  0.635  or  larger  than  1 .576.  This  condition  results  in 
the  following  approximation  for  the  maximum  space  step  (Fread,  1988): 

A*< - , (1.21) 

l+2\Aj-AJ+l\/A 

where  A = AJ+  x for  a contracting  reach  and  A = A:  for  an  expanding  reach, 

and  L is  the  original  distance  between  adjacent  cross  sections  at  space  indices  j 

and  j+ 1. 
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Figure  1.7  Projection  of  a meandering  stream  corridor  onto  a straight  line. 
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Each  cross  section  is  identified  by  a cross  section  label  and  a river  kilometer.  In 
CONCEPTS- 1.0  river  kilometer  has  to  increase  in  streamwise  direction. 
Station  and  elevation  describe  the  profile  of  the  cross  section  (Figure  1 .8). 
Station  values  must  increase  from  left  to  right.  The  user  needs  to  identify  left 
floodplain,  left  bank,  streambed,  right  bank,  and  right  floodplain. 
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Figure  1.8  Representation  of  cross  section  geometry. 


Flow  Resistance 

CONCEPTS  uses  Manning  n to  parameterize  the  resistance  to  flow  in  a 
watercourse.  Manning  n represents  the  effect  of  roughness  elements  on  the 
floodplain,  channel  banks,  and  bed  (Figure  1 .9).  Also,  small  eddy  losses  due  to 
mild  expansion  and  contraction  of  reaches  as  well  as  river  bend  losses  are  often 
included  as  components  of  the  Manning  n. 

Manning  n varies  with  the  magnitude  of  flow.  For  example,  as  flow  increases 
and  more  portions  of  the  bank  and  overbank  become  flooded,  the  vegetation  at 
these  elevations  causes  an  increase  in  the  resistance  to  flow.  There  are  some 
basic  references  available  for  estimating  n values:  e.g.,  Cowan  (1956),  Chow 
(1959),  USSCS  (1975),  USFHWA  (1979),  Coon  (1998).  Table  4.13  lists  n 
values  for  various  types  of  channels  from  Chow  (1959). 

CONCEPTS  allows  the  user  to  input  different  Manning  n values  for 
streambed,  left  and  right  banks,  and  left  and  right  floodplains.  It  uses  the 
method  of  Garbrecht  ( 1 990)  to  compute  an  equivalent  friction  factor  for  the 
wetted  section  of  the  cross  section.  The  flow  area  comprises  sections  affected 
only  by  either  left  floodplain,  left  bank,  channel  bed,  right  bank,  or  right 
floodplain  (see  Figure  1.10).  The  user  can  select  from  two  methods  to  divide 
the  sections: 

Method  1 . The  boundaries  between  the  sections  are  vertical. 

Method  2.  The  boundaries  between  floodplain  and  bank  sections  are  vertical, 
and  those  between  bank  and  bed  sections  bisect  the  angle  between 
bed  and  bank  profile. 
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Figure  1.9  Example  of  roughness  elements  on  a cross  section:  trees  and  grasses  on 
the  floodplain,  grasses  on  the  left  bank,  and  rocks  and  sediment  particles 
on  the  bed. 


bankfull 


right 

bank 


elevation 


left  bank 


Figure  1.10  Segmentation  of  flow  area  to  compute  effective  friction  factor.  Method  1 
uses  vertical  boundaries  between  the  subsections.  Method  2 uses  vertical 
boundaries  between  floodplain  and  bank  sections,  and  boundaries 
between  bed  and  bank  sections  that  bisect  the  angle  between  bed  and 
bank  profiles. 
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We  compute  the  equivalent  roughness,  ne,  as 

AR2/3 

K 


(1.22) 


where  flow  area  A and  conveyance  K are  calculated  as  a summation  of  the  flow 
area  and  conveyance  of  each  section  (left  floodplain,  left  bank,  bed,  right  bank, 
right  floodplain) 

A = Y,  Ai  (1.23a) 

sections 

K = £ K>  = X ,7 AiRi  (123b) 

sections  sections  i 

and  the  equivalent  hydraulic  radius  Re  is  a conveyance-weighted  average  over 
the  sections 

Re  = 2 I K,R,  (1.23c) 

sections 

Equations  (1.22)  through  (1.23c)  assure  continuity  and  a smooth  increase  in 
total  conveyance  as  the  water  surface  rises  above  overbank  elevation.  For  the 
case  of  uniform  flow,  the  equivalent  one-dimensional  representations  (1.22) 
and  (1 ,23c)  of  roughness  and  hydraulic  radius  yield  the  same  total  discharge  as 
that  obtained  by  summing  the  flow  rates  of  each  section  (two-dimensional 
approach).  The  equivalent  hydraulic  radius  is  not  necessarily  a monotonous 
increasing  function  of  flow  depth.  A decrease  may  occur  as  flow  moves  from 
predominantly  main  channel  to  predominantly  overbank  flow. 

Figure  1.11  shows  flow  area,  conveyance,  equivalent  hydraulic  radius,  and 
equivalent  Manning  n for  the  cross  section  shown  in  Figures  1 .9  and  1.10,  and 
using  ‘Method  2’  (cf.  Figure  1.10)  to  segment  the  cross  section.  The  following 
n values  were  assigned  to  the  sections:  left  floodplain,  0.20;  left  bank,  0.08; 
bed,  0.03;  right  bank,  0.03;  and  right  floodplain,  0.20.  In  this  case  equivalent 
hydraulic  radius  is  a monotonous  increasing  function  of  flow  depth  because  the 
floodplain  is  relatively  narrow. 


Computation  of  Flow  Variables 

The  flow  variables:  flow  area,  hydraulic  radius,  wetted  perimeter,  top  width, 
conveyance,  and  Manning  n vary  with  flow  depth,  and  have  to  be  updated 
when  flow  depth  changes  during  a runoff  event.  It  is  computationally  time- 
consuming  to  recalculate  these  flow  variables  at  each  time  step.  There  exist 
two  methods  to  represent  the  flow  variables: 

1 look-up  tables,  and 

2 power  functions. 

Look-up  tables  store  values  of  flow  variables  at  a series  of  flow  depths  that 
range  from  zero  at  the  channel  bed  to  the  flow  depth  at  the  valley  elevation. 
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Figure  1.11  Flow  area,  conveyance,  equivalent  hydraulic  radius,  and  equivalent 

Manning  n as  a function  of  flow  depth  for  the  cross  section  shown  in  Figure 
1 .9.  Manning  n values  assigned  to  the  sections  are:  floodplain  n= 0.20,  left 
bank  n=0.08,  bed  and  right  bank  n= 0.03. 
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Finding  the  depth  interval  in  the  look-up  table  that  contains  a given  flow  depth 
and  interpolating  provide  the  values  of  flow  variables. 

Power  functions  represent  the  flow  variables  as  (Garbrecht,  1 990) 

flow  variable  = a hP  (1-24) 

where  a and  P are  coefficient  and  exponent  of  the  power  function, 
respectively.  Power  functions  may  accurately  represent  flow  variables  if  they 
vary  monotonically  with  depth,  as  Figure  1.11  depicts.  However,  the  values  of 
flow  variables  may  locally  vary  markedly  if  the  cross-sectional  profile  is  quite 
irregular,  which  cannot  be  represented  by  power  functions.  Also,  hydraulic 
radius  may  decrease  when  the  water  surface  rises  above  the  overbank 
elevation,  before  increasing  further. 

Look-up  tables  are  more  accurate  than  power  functions  and  are,  therefore,  used 
by  CONCEPTS. 
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Figure  1.12  Preissmann  scheme  in  x-t  space:  •,  node  with  known  values  of  the 

dependent  variables  Q and  y;  o,  node  with  unknown  values  of  0 and  y; 
and  +,  point  at  which  the  dependent  variables  and  their  derivatives  are 
approximated. 


At 


Discretization 


Generalized  Preissmann  Scheme 

There  exist  many  approaches  to  approximate  the  system  of  governing 
equations:  (1.1)  and  ( 1 .2)  for  the  dynamic  wave  model  and  (1.1)  and  ( 1 .4)  for 
the  diffusion  wave  model.  CONCEPTS  uses  the  generalized  Preissmann 
scheme  (Abbott  and  Basco,  1 989).  The  advantages  of  the  Preissmann  scheme 
are  among  others: 

• It  has  been  extensively  used  by  many  scientists  with  good  results. 

• It  is  simple,  robust,  and  compact. 

• It  allows  a variable  spatial  grid. 

• It  is  implicit  in  time. 

This  scheme  is  also  known  as  the  box  implicit  scheme  because  of  the  way  it 
approximates  the  hydraulic  variables  (see  Figure  1.12).  Implicit  time 
integration  allows  large  time  steps,  reducing  the  computational  effort 
associated  with  long-term  simulations. 

The  generalized  Preissmann  scheme  approximates  any  variable / and  its  spatial 
and  temporal  derivatives  as  (see  Figure  1.12): 

/(*,  t ) = etiRZ-VV  +(\  + - 0)[v| tjf+  j + (1  - H r)Jf]  (1.25a) 

p i + l _ p i p i + l _ pi 

x,  t)  = v'l±±—l±±  + ( 1 - V fj  - '/  (1.25b) 

At  At 
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a pi  + \ _ fn+  \ 

JblXO  = ■' 

ox  Ax 


P'  - P1 
+ (l-9)-v+l  V 
Ax 


(1.25c) 


where  j is  space  index  denoting  the  location  of  a cross  section,  n is  time  index 
denoting  a time  level,  and  9 (9  < 9 < 1 ) and  \\i  (9  < \\i  < 1 ) are  temporal  and 
spatial  weighting  factors,  respectively  to  provide  model  flexibility  with  respect 
to  numerical  stability  and  convergence  of  the  solution.  Setting  the  spatial 
weighting  factor  \\i  = 9.5  reduces  the  scheme  to  Preissmann’s  original 
scheme  (Preissmann,  1961). 


Approximations  of  Governing  Equations  of  Open-Channel 
Flow 

It  is  convenient  to  express  temporal  variations  of  hydraulic  variables  in 
variations  of  the  dependent  variables  stage  and  discharge,  for  example  the 
flow  area  at  cross  section  j and  time  level  n+ 1 reads 


A"+  1 = A?  + (— 
J J ^dy 


A Vj  + 


8QJ, 


(1.26) 


where  the  differential  changes  in  stage  and  discharge  are  Ay  = yn + 1 —yn  and 
A Q = Qn  + 1 - Qn , respectively. 

Substituting  the  above  approximations  of  variables  and  their  derivatives  in  the 
governing  equations  of  open-channel  flow  (1.1)  through  ( 1 .3)  yields  a 
nonlinear  system  of  equations.  We  can  linearize  the  equations  by  assuming 
that  the  differential  changes  in  stage  and  discharge  are  much  smaller  than  the 
respective  values  of  stage  and  discharge.  Consequently,  we  obtain  a system  of 
finite  difference  continuity  and  momentum  equations  expressed  in  the  unknown 
quantities  Ay.,  Ay-+  , , A Qj,  and  AQJ+  ^ of  the  form: 

AjAyj  + BjAQj  + CjAyj  + , + DjAQj  + , = Ej  (1.27a) 

Fj*yj + GAQi + HAyj+ 1 + , = jj  ^ 

where  the  coefficients  Aj  through  are  flow-depending  coefficients.  The 
coefficients  Aj  through  Ej  are: 

A,  - — ^-B"  (1.28a) 

J At  J 

B,  = -f  (1.28b) 

J Ax 

ci  = lBP  i f1-28'* 

Dj  - ± (1.28d) 
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9 r ^-i* 

Ax 


-T-[Q}+i 


j j 


(1.28e) 


AxL^+1  ^ 

In  the  case  of  the  dynamic  wave  model  the  coefficients  Fj  through  Jj  are: 


f = 20(1  \\f)  llJ  rgn 

J Ax  (hd)fQj+1 


Q '/] 


(1.280 


+ 


Ax 


«?(¥(«,"+ i)2  + U-V)(«?)2) 


e 


+ g—[( i - V)W+ 1 -y'j)- (vAj+ 1 +0 -vM/)] 


gQ 


2(  i - v)(vm;+  i + (i  - 


7 Vy 


-(1  - v|/)^"[v|/(^)"+  , + (1-V)(S,)»] 


Fj 


G;  = 


1 -y  20 
A/  Ax 


(V«y  + 1 +(1  -V|/)Mf)-(l  -V}/) 


(<) 


Qr 


+ 1 


QL 


aj 


(1.28g) 


Ax  J 1 J Qn 
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+ 2g0(l  -v|/)(v|/^;+  , +(1  -y)A?)^LL 


Qj 
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//  = .ll.Ll  ~-(Qn  , - On) 

• A*(Arf,7+1(e'+'  G'> 


_0_ 

Aa' 

- 2 v|/ 
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77 


+ 2e 
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+ 2g0i|/(iM/  + i + ( 1 - vj/^j7)— ^ 
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(1.28h) 


(1.28i) 
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(1.28J) 


j,  = - - e," ) - f.(  e;+ , - e;+ , > 


At  J J At 


-T-[9(V“/+i  + 0 -\|/)w*)  + (1  _ 0)(M,Mf+  i + 0 -¥)“")] 


Ajc 


x[(i-0)(0;+1-e;)  + 0(0;+i-e;)] 


7+1  ^7 


+ [0<¥(«;+ , )2 + ( i - ¥>(»;)2) + < i - ex  V( »,"+ , >2 + ( i - ¥)(»,")')] 


X J-[(  1 - 0)(Af+  , - A ”)  + 0(Aj+  , -A‘)\ 


- g[6(  V-4*+  , + ( 1 - VK')  + ( 1 - 6 )(yA\  , + ( 1 - VM/)] 


x a! [<  1 “ e)°7+ 1 ~yJ)  + 0(f + 1 )] 


+ 


Q(y(Sf)J+]+(\-y)(Sf)*) 


+ (l-Q)(v(Sf)”+l+{l-y)(Sf)”)] 


f j 


where  hd  = A/ B is  hydraulic  depth. 

In  the  case  of  the  diffusion  wave  model  the  coefficients  Fj  through  J;  are: 


fj  - -0 


Ax  T'  KJ 


(.S',)" 

G,  = 20(1  — «/) — i-i- 

' Q'l 


Hj  = 0 


2 viz 

Ax  KJ\i  vayy+l 

/ = 20vi£A±i 
~/  + 1 


(1.28k) 


(1.281) 


(1.28m) 


(1.28n) 
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JJ  - - byb  T7(y"+ 1 -y"}  (128o) 

-*[v{Sf)*+l+(l-v)(Sf)J] 

-(i-0)[vl/(S/);+I  +(\-y){sf)p 

We  must  solve  (1.27a)  and  (1.27b)  iteratively  to  advance  the  solution  from 

time  f to  time  C+1.  The  superscript  * in  (1.28e),  (1.28j),  and  (1.28o)  denotes 
that  the  flow  variable  is  evaluated  using  the  solution  obtained  in  the  latest 
iteration  step. 

Solution  Method 


Double  Sweep  Algorithm 

The  system  of  finite  difference  equations  ( 1 .27a),  ( 1 .27b),  ( 1 .39a),  and  ( 1 ,39b) 
forms  a pentadiagonal  matrix.  CONCEPTS  successively  applies  the  double 
sweep  method  or  Thomas  algorithm  to  solve  the  pentadiagonal  matrix  (Roache, 
1977).  The  computed  increments  Ay  and  A Q are  used  to  update  the  flow 

variables  at  time  t , and  consequently  the  coefficients  Ej  and  Jy  The  iteration 
method  is  stopped  when  the  solutions  of  y and  Q have  converged,  that  is 
Ah  — > 0 and  A Q — > 0 . The  paragraphs  below  briefly  discuss  the 
implementation. 

We  introduce  the  auxiliary  variables  S.-  and  T:  to  linearly  relate  the  unknown 
differential  increments  A Vj  and  A Qj  as: 

A Qj  = SjAyj+Tj  (1.29) 


Substituting  (1.29)  into  (1.27a)  and  (1.27b),  and  eliminating  A y - yields 

^Qj+  i = Sj+  i + i + Tj+  i 


(1.30) 


where 

1 i i b.sj/  i / ; • <;  .v.  >/> 
J + ' (AJ  + BJSJ)Ir(FpaJSJ)DJ 


(1.31a) 

(1.31b) 


In  the  first,  forward  sweep  we  apply  (1 .3  la)  and  (1 .3  lb)  recursively,  with  j 
varying  from  1 to  J- 1 with  J being  the  number  of  cross  sections  of  the 
schematized  stream  corridor.  The  upstream  boundary  conditions  specify  the 
values  of  S\  and  T\.  After  completion  of  the  first  sweep,  the  calculated  values 
of  Ay  and  T j and  the  downstream  boundary  condition  enable  the  computation  of 
Ay j and  A Qj.  In  the  second,  backward  sweep  we  compute  A y • from 
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(f,-g,7')-(CAVl+D,Ag,+  |) 


■ J 


(A]  + BjSj) 


(1.32) 


and  A Q-  from  (1 .29)  by  applying  the  equations  recursively,  with  j varying  from 
J- 1 to  1 . 


Treatment  of  External  Boundary  Conditions 

Discharge  hydrographs , Q]  = Q(t) , are  imposed  at  the  upstream  boundary 
of  the  stream  corridor.  The  coefficients  S)  and  Tj  read 

A,  = 0 and  Tx  = Q'{+]  -Q\  (1.33) 


At  the  downstream  boundary  a rating  curve  is  imposed.  The  user  can  specify  a 
rating  curve  when  one  is  available,  or  CONCEPTS  can  calculate  a loop-rating 
curve  based  on  local  flow  conditions.  The  user-specified  rating  curve  has  the 
form  of  a multi-segment  power  function: 


Q = Ah)  = 


a j h 1 h < h , 

a 2h^2  hx  <h<h2 


aNrh  h>hN  i 


(1.34) 


where  Nrs  is  the  number  of  segments  comprising  the  rating  curve,  and  oq,  (3,, 
and  hj  are  coefficient,  exponent,  and  breakpoint  of  rating  curve  segment  i.  If 
there  is  only  one  segment,  h\  should  be  set  to  a very  large  number. 


A loop-rating  boundary  condition  allows  the  unsteady  wave  to  pass  the  outlet 
with  minimal  disturbance  of  the  outlet  itself.  We  assume  the  following  relation 
between  discharge  and  stage  at  the  outlet 

q = m = K(h)jsf  (i.35) 


Rewriting  the  momentum  equation  (1.2)  yields  an  expression  for  the  friction 
slope  in  (1.35) 


.f 


'dv 
_ dx 


+ J-(a2+i/ehY| 

gA^  dt  dx^  A D _ 


(1.36) 


In  the  case  of  the  diffusion  wave  model  we  only  retain  the  water  surface  slope 
in  ( 1 .36).  We  approximate  ( 1 .36)  as 


Sf 


ynj 


-_y± 

Ax 


+ 


i (Q5-QT 


§A') 


At 


+ 


(: Q2/A)j~(Q2/A)A 
Ax 


(1.37) 


The  flow  depth  increment  at  the  outlet  is  then 
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Ay/ 


Tj+QJ-f 

X-Sj 

dy  J 


(1.38) 


Treatment  of  Internal  Boundary  Conditions 

Similar  to  the  open-channel  flow  equations,  we  express  the  approximated 
continuity  and  dynamic  equations  of  flow  at  hydraulic  structures  in  terms  of 
differential  changes  of  the  dependent  variables: 

A0,  -A03  = Q%-Q*x  (1.39a) 


A0, 


df  * dfA 

— Ay,  - Ti-Ay3 


dy 


dv 


y 3 


f-Q\ 


(1.39b) 


where/ represents  the  dynamic  equation  of  the  structure.  Equations  ( 1 ,39a)  and 
( 1 ,39b)  are  equivalent  to  ( 1 .27a)  and  ( 1 .27b). 


We  can  derive  analytical  expressions  of  the  derivatives  of / with  respect  to  rq 
andy3.  However,  if  the  structure  is  a control , that  is  the  flow  at  the  structure  is 
critical,  this  approach  is  not  very  robust.  In  this  case  we  approximate  the 
derivative  of / with  respect  tojq  as 


d£_  = 2A£> 

fy\  y\(Q  + &Q)-y\(Q-&Q) 


(1.40) 


where yx(Q)  is  the  headwater  equation  and  A Q is  a pre-set  discharge  increment. 


One-Equation  Diffusion  Wave  Model 


The  Preissmann  scheme  applied  to  the  dynamic  wave  or  diffusion-wave  model 
will  predict  unrealistic  solutions  of  the  dependent  variables  when  there  are 
large  increases  of  discharge  and  flow  depth  at  the  upstream  end  of  a 
computational  section  (e.g.,  when  a flood  wave  enters  the  computational 
section).  We  can  readily  show  this  through  the  discretized  continuity  equation 
from  which  we  obtain  the  following  expression  for  the  differential  change  of 
stage  at  the  downstream  end  of  a reach  AyJ+  j : 


= - ~LB’’ Ay,  + ±(AQ.-AQJ+])-  -L( Q« 


At 


Ax 


Ax 


'j+  l 


Qjn ) 


(1.41) 


When  a flood  wave  enters  the  reach.  Ay,  and  AQj  are  positive  and  may  be  quite 
large.  Assuming  AQf ■ + x * 0 and  Qj\  t - Q"  ~ 0 , A yy-+1  is  negative  and  quite 
large.  As  a consequence,  the  predicted  flow  depth  hJ+  j may  become  negative. 
If  this  condition  occurs,  CONCEPTS  switches  to  the  ‘one-equation’  diffusion 
wave  model. 


We  can  reformulate  the  diffusion  wave  model  as  a single  equation  in  either 
discharge  or  stage  by  eliminating  one  of  the  dependent  variables.  We  can 
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eliminate  stage  by  taking  the  derivative  of  ( 1 . 1 ) with  respect  to  x and 
subtracting  the  derivative  of  ( 1 .4)  with  respect  to  t.  This  yields  the  one- 
equation  diffusion  wave  model : 


d£ 


dt 


-D& 

dx2 


cq 


(1.42) 


where  the  celerity  c is 


- _ Q dK 
C BKdh 


(1.43) 


which  is  also  known  as  the  kinematic  wave  celerity,  the  celerity  c is 

- DdB 

c = c + 

B dx 


(1.44) 


and  the  diffusion  coefficient 


D 


= JL 

2 BS 


(1.45) 


CONCEPTS  uses  a fractional  step  method  (e.g.,  see  Yanenko,  1971)  to 
integrate  ( 1 .42).  In  the  first  step  we  employ  the  method  of  characteristics 
(Abbott,  1966)  to  integrate  the  advection  term  and  add  lateral  inflow 

dO  + dQ  - 
dt  dx 


(1.46) 


the  solution  of  which  is  the  initial  condition  for  the  integration  of  the  diffusion 
term 

= Dlo 

8‘  dx2 

for  which  we  use  an  explicit,  central  finite-difference  scheme  (e.g.,  see  Abbott 

2 

and  Basco,  1989).  If  At  > 0.5  Ax  /D , we  iteratively  integrate  (1.47)  using  a 

2 

time  step  At/Nt,  where  the  coefficient  Nt  > 2AtD/ Ax“ . The  above  method 

ensures  that  discharge  is  always  positive  and  its  peak  attenuates  in  streamwise 
direction. 


Finally,  differential  changes  in  water-surface  elevation  are  computed  as 

Ay  = -t±2- 

dy 


§,Jsf 


(1.48) 


CONCEPTS  also  switches  to  the  one-equation  diffusion  wave  model  if 
A Q/Q  >0(  1 ) or  Ay/y  >0(1). 
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This  chapter  presents  the  one-dimensional  sediment  transport  and 
river-bed  evolution  submodel  used  in  CONCEPTS- 1.0.  You  will 
find  a brief  description  of  the  conceptual  basis  for  the  mathematical 
formulations  of  sediment  transport  in  alluvial  channels  and  how 
streamwise  variations  in  sediment  transport  affect  the  evolution  of 
the  streambed.  This  chapter  also  gives  a detailed  presentation  of  the 
numerical  approximation  of  the  mathematical  model  of  streambed 
evolution. 


Theory 


Sediment  Transport 

Sediment  transport  rates  are  a function  of  flow  hydraulics,  bed  composition, 
and  upstream  sediment  supply.  The  composition  of  the  channel  bed  may 
change  as  particles  are  eroded  from  or  deposited  on  the  bed,  thereby  changing 
flow  hydraulics  and  fractional  transport  rates. 

To  model  these  physical  processes,  one  commonly  distinguishes  two  layers 
across  the  water  column  and  one  or  more  layers  covering  the  streambed  (see 
Figure  2.1;  e.g.,  Armanini  and  Di  Silvio,  1988;  van  Niekerk  et  al.,  1992).  The 
two  layers  across  the  water  column  are  (see  Figure  2.2): 

1 a layer  near  the  bed  where  sediment  particles  roll,  slide,  or  saltate,  and  are 
transported  as  bed  load , and 

2 a layer  spanning  the  remainder  of  the  water  column  where  particles  are  in 
suspension  and  are  transported  as  suspended  bed-material  load  or  wash 
load. 
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Figure  2.1  Multi-layer  approximation  of  sediment  transport. 


Figure  2.2  Modes  of  sediment  transport:  bed  load  and  suspended  load. 


grains  held  in  suspension 

o 

by  turbulent  eddies  in  the 

o° 

suspended  load  layer 

Sediment  particles  exchange  between  the  bed  and  suspended  load  layers  and 
the  bed.  Governing  equations  for  each  layer  and  empirical  expressions  of  the 
fluxes  between  the  layers  and  the  bed  enable  us  to  compute  the  sediment 
transport  rates  in  each  layer.  Conversely,  we  can  combine  the  suspended  and 
bed  load  layers  into  a single  total  load  layer.  The  fluxes  between  the 
transporting  layers  disappear  and  only  the  sediment  flux  between  the  bed  and 
fluid  remains. 

This  two-layer  approach  may  be  a more  accurate  description  of  sediment 
transport  mechanics  because  of  the  difference  in  suspended  and  bed  load 
transport  modes.  However,  the  expressions  for  the  sediment  fluxes  between  the 
layers  are  empirical  and  are  mainly  derived  from  data  acquired  from  sand-bed 
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rivers  and  flume  experiments  (e.g.,  Garcia  and  Parker,  1991).  The  one-layer 
approach  cannot  distinguish  the  different  behavior  of  sediment  particles  in  the 
suspended  and  bed  loads.  However,  it  is  computationally  more  efficient 
because  of  its  simplicity,  which  is  important  for  long-term  simulations  of  the 
evolution  of  stream  systems.  CONCEPTS,  therefore,  characterizes  the 
sediment  transport  as  total  load. 

The  sediment  flux  between  streambed  and  overlying  fluid  not  only  depends  on 
flow  conditions  and  size  fractions  available  in  the  bed  or  transported  by  the 
flow,  but  also  on  the  interparticle  forces  in  the  presence  of  clay  particles  on  the 
streambed.  The  entrainment  rate  of  sediment  particles  on  cohesive  streambeds 
markedly  differs  from  that  of  cohesionless  streambeds. 


Governing  Equations 


Mass  Balance 

Mass  conservation  of  sediment  by  size  fraction  is 


dCk  + duCk 
dt  dx 


= Ek-Dk  + qs 


(2.1) 


where  t is  time,  x is  distance  along  the  channel  centerline,  u is  flow  velocity,  E 
is  entrainment  rate  of  particles  from  the  bed,  D is  deposition  rate  of  particles 
onto  the  bed,  qs  is  rate  of  sediment  inflow  from  streambanks  and  fields 
adjacent  to  the  channel,  the  subscript  k denotes  the  k-ih  size  class,  and  C is  the 
sediment  mass  with  dimension  length  . The  sediment  mass  C is  defined  as: 


C = 


1 

1 xlO6  — 
7 


(2.2) 


where  c is  point  concentration  in  ppm  by  weight,  ys  is  specific  weight  of 
sediment,  and  y is  specific  weight  of  water. 


Figure  2.3  shows  a definition  sketch  of  the  variables  used  to  determine 
sediment  transport.  Temporal  variations  in  bed-material  area  are  given  by 


dAb 

U-M-f  = Dk-Ek  (2.3) 


where  X is  porosity  and  Ab  is  cross-sectional  area  of  the  mixing  layer. 

The  entrainment  and  deposition  rates  differ  for  cohesive  and  cohesionless  bed 
material  and,  thus,  are  computed  using  different  methods. 


Non-Equilibrium  Adjustment  of  Cohesionless  Bed-Material 
Transport 

For  cohesionless  streambeds  we  use  the  formulation  proposed  by  Bennett 
(1974)  analogous  to  that  of  Foster  and  Meyer  (1972)  who  assume  that  the  local 
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Figure  2.3  Definition  sketch  of  sediment  transport  variables. 


Figure  2.4  Progress  of  sediment  mass  C toward  equilibrium  sediment  mass  C for  the 
case  of  constant  flow  conditions. 


erosion  or  deposition  rate  is  proportional  to  the  difference  between  the 
sediment  transport  rate  and  sediment  transport  capacity: 

Ek~Dk=  Uct-Ck)  (2.4) 

*/< 

where  Tis  a time  scale  representing  the  adjustment  rate  of  the  sediment  mass  C 

to  C,  which  is  the  equilibrium  sediment  mass  and  is  a function  of  local, 
instantaneous  flow  conditions.  Figure  2.4  shows  an  example  of  how  C 

approaches  C for  the  case  of  erosion  and  deposition  under  constant  flow 
conditions.  Various  formulations  of  T exist  in  literature.  For  particles 
transported  in  suspension , Armanini  and  Di  Silvio  (1988)  suggest 
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1/6 


CO 


(2.5) 


where  co  is  particle  fall  velocity,  h is  flow  depth,  a is  thickness  of  the  surface  or 
active  layer,  and  u*  is  shear  velocity.  For  sediment  particles  transported  as  bed 
load , Phillips  and  Sutherland  ( 1 989)  propose 


where  aL  is  step-length  constant,  0 is  Shields  parameter,  0C  is  critical  Shields 
parameter,  and  d is  particle  diameter.  Phillips  and  Sutherland  (1989) 
experimentally  found  4000  < aL  < 9000 . 

CONCEPTS  employs  formulation  (2.5).  Although  (2.5)  was  not  derived  for 
sediment  particles  transported  as  bed  load,  it  correctly  tends  to  C = C 
because  T « Ax/u  for  coarse  particles,  where  Ax  is  the  space  step  used  by 
CONCEPTS. 

Entrainment  and  Deposition  of  Cohesive  Bed-Material 

The  transport,  deposition,  and  erosion  processes  of  cohesive  sediments  are 
extremely  complex  due  to  their  highly  varying  properties  and,  therefore, 
behavior  when  their  environment  changes.  Cohesive  beds  in  incised  channels 
may  be  very  firm  and  highly  consolidated.  Ariathurai  and  Arulanandan  (1978) 
give  erosion  rate  as 


where  e is  erosion-rate  constant,  B is  wetted  width  of  streambed,  xb  is  bed 
shear  stress,  and  xe  is  shear  strength  of  the  bed  material  (see  Figure  2.5).  The 
erosion-rate  constant  and  bed  shear-strength  vary  with  type  of  sediment,  water 
content,  total  salt  concentration,  ionic  species  in  the  water,  pH,  and 
temperature  (Mehta  et  al.,  1989).  For  soft  (water  content  well  above  100%), 
partially  consolidated  beds,  Parchure  and  Mehta  (1985)  found  the  erosion  rate 
to  be 


where  e^is  floe  erosion  rate  and  e is  a rate  constant.  When  the  excess  shear 
stress  becomes  large,  the  bed  may  fail  at  some  plane  below  the  surface  and 
clumps  of  material  are  mass  eroded.  Once  entrained,  CONCEPTS  breaks  up 
the  clumps  into  their  primary  particles,  which  are  then  transported  as  wash  load 
(clay  and  fine  silt  particles)  or  become  part  of  the  bed-material  load. 

The  deposition  rate  is  commonly  given  by  Krone’s  (1962)  formulation: 


T - oll(6  - 6c)d/u 


(2.6) 


(2.7) 


(2.8) 


(2.9) 
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Figure  2.5  The  erosion  rate  of  a cohesive  streambed  as  a function  of  bed  shear 
stress. 


bed  shear  stress  xb 


where  xd  is  the  shear  stress  below  which  sediment  particles  in  transport  begin 
to  deposit.  The  term  in  between  parentheses  in  (2.9)  is  the  probability  of 
particles  sticking  to  the  bed  and  not  being  re-entrained  by  the  flow. 
CONCEPTS  sets  the  deposition  rate  to  zero  if  xb  > xd . 

Once  deposited,  cohesive  particles  or  aggregates  undergo  consolidation  and 
thixotropic  effects.  This  results  in  a vertical  stratification  of  the  bed  with 
respect  to  density  and  shear  strength,  with  both  properties  typically  increasing 
with  depth  (Mehta  et  al. , 1982).  These  processes  are  complex  and  highly 
varying  in  space  and  time.  A simple  description  of  these  processes  is  not  yet 
available.  CONCEPTS,  therefore,  assumes  that  a deposited  layer  of  cohesive 
particles  has  the  same  characteristics  as  the  bed  surface. 

Mixing  Layer 

For  graded  bed  material  the  sediment-transport  rates  depend  on  the  bed- 
material  composition,  which  itself  depends  on  historical  erosion  and  deposition 
rates.  Hirano  (1971)  divided  the  bed  into  a surface  or  active  layer  and  a 
substrate  or  subsurface  layer.  These  layers  constitute  the  so-called  mixing 
layer  (see  also  Figure  2.1).  Sediment  particles  are  continuously  exchanged 
between  flow  and  the  surface  layer.  Sediment  particles  exchange  between 
surface  layer  and  substrate  when  the  bed  scours  or  fills.  The  volumetric 
fraction  content  by  size  class  in  the  surface  layer  is  determined  by  the 
following  mass  conservation  equation: 
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spf a 

dt 


(2.10) 


where  is  the  fractional  content  by  volume  of  size  class  k in  the  surface  layer 
(I.k  PA.  = 1 ),  As  is  the  area  of  the  surface  layer,  and  Su  is  the  sediment  flux 
from  the  substrate  to  the  surface  layer  given  by: 


The  term  within  parentheses  divided  by  the  active  bed  width  is  the  downward 
displacement  of  the  lower  boundary  of  the  surface  layer,  and 


where  (3£  is  the  fractional  content  by  volume  of  size  class  k in  the  subsurface 
layer. 

The  thickness  of  the  surface  or  active  layer  is  associated  with  the  time  scale 
under  consideration,  here  the  time  step  At  of  the  time-integration  scheme 
(Rahuel  et  al.,  1989).  For  very  small  time  scales,  the  surface  layer  can  be 
considered  as  a thin  layer  containing  particles  susceptible  to  entrainment  into 
the  flow  due  to  a momentary  increase  in  the  bed  shear  stress.  If  the  time  scale  is 
of  the  order  of  the  time  it  takes  for  a bed  form  to  traverse  its  own  wavelength, 
the  thickness  of  the  surface  layer  is  approximately  the  height  of  the  bed  form. 
For  very  long  time  scales  during  which  considerable  changes  in  bed  level 
elevation  may  occur,  the  surface  layer  thickness  is  related  to  the  thickness  of 
the  layer  of  material  eroded  or  deposited. 

Scientists  have  considered  the  surface  layer  thickness  to  be  a function  of  dune 
height  or  water  depth  (Armanini  and  Di  Silvio,  1988;  Rahuel  et  al .,  1989),  or 
grain  size  (Borah  et  al.,  1982;  van  Niekerk  et  al.,  1992;  Cui  et  al.,  1996). 
CONCEPTS  sets  the  thickness  of  the  surface  layer  to  10%  of  the  flow  depth. 
Rahuel  et  al.  (1989)  propose  a thickness  of  the  surface  layer  between  10%  and 
20%  of  the  flow  depth.  Armanini  and  Di  Silvio  (1988)  suggest  a minimum 
thickness  of  the  surface  layer  of  5%  of  the  flow  depth. 

Sediment  Transport  Capacity 

The  flow  sediment  transport  capacity  is  the  transported  sediment  load  under 
equilibrium  conditions  (uniform  flow  and  E = D , that  is,  no  net  erosion  or 
deposition).  CONCEPTS  uses  a modification  of  the  sediment  transport 
capacity  predictor  SEDTRA  developed  by  Garbrecht  et  al.  (1996).  SEDTRA 
calculates  the  total  sediment  transport  by  size  fraction  for  twelve  predefined 
size  classes  (0.01  mm  < d < 50  mm ) with  a suitable  transport  equation  for 


(2.11) 


Pit  if^SO 
Pi' 


(2.12) 
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Table  2.1  Sediment  size  fractions  and  corresponding  transport  equations. 


upper  representative 
size  bound  diameter 

class  (mm)  (mm) 

1 0.010 


clay  - very  fine  silt 
fine  - medium  silt 
medium  - coarse  silt 
fine  sand 

medium  - coarse  sand 
very  coarse  sand 
very  fine  gravel 
fine  gravel 
fine  gravel 
medium  gravel 
coarse  gravel 
coarse  gravel 
very  coarse  gravel 


description 


transport  equation 

wash  load 


2 0.025  0.016 

3 0.065  0.040 

4 0.250  0.127 

5 0.841  0.458 

6 2.000  1.297 

7 3.364  2.594 

8 5.656  4.362 

9 9.514  7.336 

10  16.000  12.338 

1 1 26.909  20.749 

12  38.055  32.000 

13  50.000  43.713 


0.016 

0.040 

0.127 

0.458 

1.297 

2.594 

4.362 

7.336 


Laursen 

Laursen 

Laursen 

Yang 

Yang 


Meyer-Peter  & Mueller 
Meyer-Peter  & Mueller 
Meyer-Peter  & Mueller 
Meyer-Peter  & Mueller 
Meyer-Peter  & Mueller 
Meyer-Peter  & Mueller 
Meyer-Peter  & Mueller 


each  size  fraction:  Laursen  (1958)  for  silt  size  classes,  Yang  (1973)  for  sand 
size  classes,  and  Meyer-Peter  and  Mueller  (1948)  for  gravel  size  classes.  We 
have  added  a thirteenth  sediment  size  class  for  sediment  particles  smaller  than 
10  pm.  Once  entrained,  the  flow  carries  this  size  class  as  wash  load  without 
deposition.  Table  2.1  lists  the  13  size  classes  and  their  transport  equations. 
SEDTRA  takes  into  account  the  interdependence  between  size  fractions  for 
initiation  of  motion  by  a critical  sediment  diameter  for  each  size  fraction  that  is 
a function  of  the  sediment  mixture. 

The  Shields  curve  is  a reliable  predictor  of  the  flow  strength  necessary  for  the 
initiation  of  cohesionless  particles  with  a narrow  size  range.  For  sediments 
with  a widely  graded  size  distribution,  however,  the  differences  in  the  critical 
flow  strength  tend  to  be  significantly  reduced.  Widely  graded  sediment  beds 
tend  to  increase  the  critical  flow  strength  for  initiation  of  the  sizes  finer  than 
the  mean  size  and  decrease  the  critical  flow  strength  of  the  sizes  coarser  than 
the  mean  size.  To  account  for  the  effect  of  the  mixture  on  the  critical  flow 
strength  of  the  individual  size  fractions,  SEDTRA  defines  the  critical  diameter 
for  initiation  of  each  of  the  twelve  size  fractions  as  (Kuhnle  et  al.,  1996) 


where  d is  the  mean  size  of  the  bed  material,  and  % is  a hiding  coefficient.  The 
hiding  coefficient  ranges  from  0 to  1.  For  x = 1 , the  mean  size  of  the 
sediment  is  the  critical  diameter  for  all  size  fractions  and  all  fractions  tend  to 
move  at  the  same  flow  strength.  For  x = 0 , each  size  fraction  behaves 
independently  of  the  others  and  the  dk  for  each  size  fraction  is  used  to  calculate 
flow  strength  at  which  motion  begins. 


(2.13) 
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Figure  2.6  Measured  (Brownlie,  1981 ) versus  computed  (SEDTRA)  sediment 

transport  capacity.  The  black  line  is  the  line  of  perfect  agreement  and  the 
magenta  line  is  the  regression  line. 

• Laursen  (1958)  ©Yang  (1973)  ® Meyer-Peter  & Mueller  (1948) 
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Extensive  testing  of  SEDTRA  found  it  to  be  applicable  to  channels  with 
widely  graded  sediment  distributions  and  to  channel  networks  with  variable 
sediment  characteristics  (Garbrecht  et  al.,  1996).  Garbrecht  et  al.  (1996)  tested 
SEDTRA  against  Brownlie’s  (1981)  data  set.  These  data  consist  of  5,263 
laboratory  and  1,764  field  measurements  made  under  equilibrium  or  near- 
equilibrium conditions  and  a wide  range  of  channel,  flow,  and  bed-material 
characteristics.  Garbrecht  et  al.  ( 1 996)  did  not  use  data  having  one  or  more  of 
the  following  values  for  testing:  water  temperature  above  35  Celsius,  sediment 
specific  gravity  other  than  between  2.4  and  2.8,  energy  slopes  less  than 
0.00001 , gradation  greater  than  1 .5  (laboratory  data)  or  2.0  (field  data),  and 
channel  width-to-depth  ratios  less  than  5. 

Figure  2.6  shows  total  measured  versus  total  computed  sediment  transport. 
About  80%  of  the  computed  laboratory  data  are  within  a factor  of  two  of  the 
measured  values,  and  about  90%  within  a factor  of  three.  About  55%  of  the 
computed  field  data  are  within  a factor  of  three,  and  about  70%  within  a factor 

of  5.  The  regression  line  plots  closely  to  the  line  of  perfect  agreement,  r2 
equals  0.85. 
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Figure  2.7  Function  3(u*/u>)  in  Laursen’s  approach  (Laursen,  1958). 


The  following  sections  present  the  transport  equations  in  SEDTRA.  The  point 
concentration  c is  related  to  sediment  mass  Cby  (2.2).  SEDTRA  applies 
transport  equations  (2.15),  (2.17),  and  (2.19)  to  individual  particle-size  classes. 
The  total  sediment  mass,  Ct,  is  defined  as: 

C,  = Jj>kCk  (2-14) 

k 

where  pk  is  the  fraction  of  sediment  in  the  &-th  size  class  available  for 
transport.  The  fraction  pk  depends  on  (3|  and  the  fraction  of  sediment  in  the 
&-th  size  class  entering  the  reach  from  upstream  and  the  streambanks. 

Laursen  (1958) 

SEDTRA  uses  Laursen’s  (1958)  transport  equation  in  the  form  as  given  in 
Vanoni  (1975) 


where  the  shear  velocity  w*  = JghSj -,  g is  gravitational  acceleration,  Sp  is 

friction  or  energy  slope,  xc  is  critical  tractive  force  as  given  by  (2. 13),  and  7 is 

a function  of  the  parameter  w*/ co , see  Figure  2.7.  Laursen  (1958)  used  the 

Manning-Strickler  equation  to  express  the  bed  shear  stress  x'  due  to  grain 
resistance  as 
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T'  = &L 

58 


d \ 1/3 
_50 ' 

h 


(2.16) 


where  d 50  is  the  sediment  size  for  which  50%  of  the  mixture  is  finer. 


Yang  (1973) 

Yang  (1973)  assumed  sediment  concentration  to  be  related  to  unit  stream 
power  uSf.  From  analysis  of  laboratory  data,  Yang  found 


logc  = 5.435  - 0.2861og—  - 0.4571og(  — 


v 


03 


(2.17) 


co  d 


+ 1 .799  - 0.4091og  — - 0.3 141og—  log  — 4 - 
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( uSf  (w5y) 


^ co 


03 


where  v is  kinematic  viscosity,  and  (uSj)c  is  critical  unit  stream  power.  The 

critical  dimensionless  unit  stream  power  is  the  product  of  the  friction  slope  and 
critical  dimensionless  velocity  (ic^co),  where 


_£  _ 


03 


2.5 

log (u*dc/\)  - 0.06 
2.05 


u*d. 

for 

1.2  < 

c 

V 

u*dn 

for 

70  < 

c 

< 70 


(2.18) 


v 


Meyer-Peter  and  Mueller  (1948) 

Meyer-Peter  and  Mueller  (1948)  proposed  the  following  relation  for  bed  load 
transport 


c 


8x10  y M*dr  3/2 
^ c-(r  Q - 0.047) 

p q 


(2.19) 


where  p is  water  density,  q is  unit  discharge,  r is  the  ratio  of  bed  roughness  and 
grain  roughness  Strickler  coefficients,  and  0 is  Shields  parameter: 


= r RSf 

(y  ,-y)dc 


(2.20) 


The  ratio  r is  determined  as 


(2.21) 


where  Mueller  determined  the  grain  roughness  as 

d]/6 
...  _ “90 


26 


(2.22) 
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Fall  Velocity 

Equations  (2.5),  (2.9),  (2.15),  and  (2.17)  require  the  fall  velocity  co. 
CONCEPTS  uses  the  fall  velocities  provided  by  the  U.S.  Interagency 
Committee  on  Water  Resources,  Subcommittee  on  Sedimentation  (USICWR, 
1963).  USICWR  (1963)  determined  fall  velocity  for  various  sediment  sizes 
and  shapes  under  different  water  temperatures.  CONCEPTS,  however,  only 
uses  the  fall  velocities  for  well-rounded  natural  sediment  particles. 

USICWR  (1963)  established  these  fall  velocities  for  particles  falling  by 
themselves  in  still  water.  However,  the  presence  of  other  particles,  turbulence 
within  the  flow,  and  particle  density  will  all  affect  fall  velocity,  but 
CONCEPTS  version  1.0  does  not  account  for  them. 
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Computation  of  Fractional  Sediment  Concentration 


Noncohesive  Streambed 

The  sediment  transport  equation  (2.1)  is  an  advection  equation  with  a source 
term  representing  a net  erosion  from  or  deposition  on  the  bed  of  sediment 
particles  and  a lateral  inflow  of  eroded  bank  material  and  eroded  soil  carried  by 
lateral  runoff.  Bank  erosion  rates  are  computed  by  the  streambank  component 
of  CONCEPTS  (see  Chapter  3 on  “Bank  Erosion  and  Channel  Widening”). 
Runoff  sediment  yield  is  user-supplied  input  data. 


There  exists  much  literature  on  methods  that  have  been  developed  to  integrate 
the  advection  equation.  Though  the  character  of  the  advection  equation  with 
constant  coefficients  is  well  understood,  its  numerical  integration  may  be  quite 
difficult,  especially  if  the  solution  displays  large  gradients  (e.g.,  Hirsch,  1988). 
Also,  if  the  coefficient  in  the  advection  term  is  varying,  as  is  velocity  in  (2.1), 
the  complexity  of  the  numerical  integration  increases.  Further,  one  may  use  a 
stable  scheme  to  integrate  the  advection  equation,  however,  to  obtain  a 
diverging  solution  because  of  the  behavior  of  the  source  term.  If  the  ratios 
A t/T  and  A x/uT,  where  At  is  time  step,  are  much  larger  than  unity,  the 
discretized  source  term  is  much  larger  in  magnitude  than  the  discretized 
unsteady  and  advection  terms.  Therefore,  the  source  term  controls  the 
sediment  concentration,  and  any  solution  method  of  the  advection  equation 
will  be  either  unstable  or  inaccurate.  The  source  term  has  to  be  approximated 
by  an  implicit  method  to  assure  a stable  solution  (LeVeque  and  Yee,  1990). 
Alternatively,  we  can  employ  a fractional  step  method  (Yanenko,  1971),  and 
solve  the  following  system  of  equations: 


dCk  duCk 
dt  dx 


(2.23a) 
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Figure  2.8  Solution  of  advection  equation  in  x-t  space. 

▲ 


dCk  i - 

— - = — (Ck-C.)  (2.23b) 

dt  Tk  k kJ 

The  advection  equation  (2.23a)  is  solved  first.  Its  solution  is  used  as  initial 
condition  of  the  initial  value  problem  (2.23b). 

The  behavior  of  sediment  mass  C can  be  studied  in  x-t  space  (see  Figure  2.8). 

The  characteristic  projected  backwards  from  point  P intersects  the  t = t"  line 
at  point  O.  Thus,  sediment  particles  travel  from  point  O to  point  P in  a time 

step  At  = t"  + - tn . Neglecting  the  source  term  we  simply  have  Cp  = C0 . 
Including  the  source  term,  C will  adjust  itself  from  CG,  approaching  a value  of 
Cp  at  a rate  given  by  timescale  T.  In  the  above  analysis  we  neglected  lateral 
inflow  of  sediments. 

CONCEPTS  solves  (2.23a)  by  the  method  of  characteristics  to  obtain  an 
intermediate  solution  C : 


Ch  = 


( 1 - Nc) CkJ  + NcClj_  , + qsAt  for  < 1 

1 “ ch  - i + 37  ch  - * + for  > 1 

c 


(2.24) 
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where  Nc  = uAt/Ax  is  the  Courant  number.  The  discretization  scheme  (2.24) 

is  first  order  in  space  and  time,  and  therefore  may  be  quite  diffusive  if 
longitudinal  gradients  of  sediment  concentration  are  large. 

The  intermediate  solution  C*k  • is  the  initial  condition  for  solving  the  ordinary 

differential  equation  (2.23b).  Many  methods  are  available  to  solve  (2.23b) 
numerically,  e.g.,  see  Gear  (1971)  or  Lambert  (1973).  However,  if  we  assume 

that  both  Ck  and  T are  constant  during  the  time  step  At,  we  obtain  the 
following  analytical  solution 

qj  1 = Ckj  + (Clj-CkJ)eM-^/T)  (2.25) 

Equation  (2.25)  yields  1 = Ckj  for  coarse  sediment  particles  (T « At), 
and  C£  +j  1 = C*k  j for  very  fine  particles  (say  wash  load,  T ->  go). 


Cohesive  Streambed 

In  the  case  of  a cohesive  streambed,  (2.7)  and  (2.9)  comprise  the  source  and 
sink  terms,  respectively.  Equation  (2.23b)  then  reads 

dC, 

nr  = E*~D*  (226) 


which  is  solved  by  the  well-known  forward  Euler  method  (e.g.,  Celia  and 
Gray,  1992) 


(Ck)f  1 = (Ck)*  + A t((Ek)J  - (Dk)J)  (2.27) 


where 


r 1 


{EX  - B 


2[(^);+(^);+1] 

(Te)/ 


- 1 


{E>k)j  E(okck 


1 - 


^[(Ti,);  + <Tf>);+']- 
(x d)i 


(2.28a) 


(2.28b) 


Computation  of  Variations  in  Streambed  Elevation 


The  relation  describing  change  in  bed  material  area  (2.3)  is  equivalent  to 
(2.23b)  and  (2.26).  Therefore,  the  change  in  bed  material  storage  per  unit 
length  of  channel,  , equals  the  difference  between  the  intermediate 

sediment  concentration  and  that  at  r"+1 : 

= (2-29> 
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Figure  2.9  Changes  in  elevation  of  the  streambed:  a)  erosion  of  a partly  wetted  bed, 
b)  erosion  of  a fully  wetted  bed,  c)  deposition  on  a partly  wetted  bed,  and 
d)  deposition  on  a fully  wetted  bed.  The  darker,  shaded  area  denotes 
deposited  or  eroded  sediments. 
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Cross  Section  Evolution 

CONCEPTS  uses  the  change  in  bed-material  storage  (2.29)  to  update  the 
elevation  of  cross-sectional  nodes  along  the  wetted  part  of  the  bed,  see  Figure 
2.9.  Figures  2.9a  and  2.9b  show  the  change  in  bed  elevation  for  the  case  of 
erosion  of  a partly  wetted  bed  and  a fully  wetted  bed,  respectively.  Figures  2.9c 
and  2.9d  show  the  change  in  bed  elevation  for  the  case  of  deposition  on  a partly 
wetted  bed  and  a fully  wetted  bed. 

CONCEPTS  shifts  the  bed  profile  up  and  down  parallel  to  the  original  bed 
profile  without  changing  the  bank  profile.  In  the  case  of  a partly  wetted  bed,  a 
new  point  at  the  intersection  of  water  surface  and  bed  profile  is  inserted. 


Mixing  Layer  Composition 


CONCEPTS  vertically  divides  the  bed  into  a surface  layer  and  several 
subsurface  layers,  which  reflect  historical  deposition  of  sediment  particles  on 
the  streambed  or  the  undisturbed  streambed.  CONCEPTS  keeps  track  of  the 
different  compositions  of  surface  and  subsurface  layers.  Rewriting  the  mass 
conservation  equation  of  the  surface  layer  (2.10)  using  (2.3)  and  (2.1 1)  yields 


m*s 

dt 


(1  -A.) 


(2.30) 
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Figure  2.10  Sorting  of  graded  bed-material  within  the  mixing  layer:  a)  erosion, 

sediment  particles  from  surface  layer  mix  with  material  from  substrate 
layers  1 and  2;  and  b)  deposition,  part  of  the  surface  layer  becomes  new 
substrate  layer  1 . 
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We  assume  that  (3j*  is  constant  during  a time  step  At,  thus 

(PP"+I  = -77Tt(P^/  + (l-X)A^t+(P;)"(A^-A^)]  (2.31) 

AS 

Figure  2.10  shows  the  mixing  of  sediment  particles  with  those  in  subsurface 
layers  if  erosion  occurs  and  the  initiation  of  new  subsurface  layers  if  deposition 
occurs.  Entrainment  of  sediment  particles  from  the  surface  layer  and  its 
ensuing  downward  displacement  causes  particles  from  subsurface  layers  to  be 
mixed  with  those  in  the  surface  layer  (2.10a). 

Deposition  of  sediment  particles  on  the  bed  leads  to  an  upward  displacement  of 
the  surface  layer  and  the  initiation  of  new  subsurface  layers  (2. 10b). 
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Boundary  Conditions 

CONCEPTS  requires  boundary  conditions  at  the  upstream  and  downstream 
boundaries  of  the  channel,  and  at  hydraulic  structures.  The  upstream  and 
downstream  boundaries  are  also  known  as  external  boundaries , whereas 
boundaries  such  as  hydraulic  structures  are  known  as  internal  boundaries. 


External  Boundary  Conditions 

At  the  upstream  boundary  the  user  must  supply  the  sediment  discharge  by  size 
class  entering  the  channel.  This  is  done  by  entering  (see  also  sections  “Run 
Control  Data”  and  “Dynamic  Upstream  Boundary  Conditions”  in  Chapter  4, 
“Input  Data”): 


1 sediment  load  as  a fraction,  fk,  of  the  local  sediment-carrying  capacity  of  the 
flow: 


fk 


£k 

Ck 


(2.32) 


or, 

2 a time  series  of  sediment  load  (in  kg/s),  Ck  = Ck(t) , in  the  file  with 
upstream  boundary  conditions. 

The  solution  method  (2.24),  (2.25),  and  (2.29)  does  not  require  a downstream 
boundary  condition.  However,  this  means  that  only  processes  occurring  in  the 
modeling  reach  determine  the  evolution  of  the  outlet.  This  may  lead  to 
inaccurate  results.  The  user  has  the  option  to  adjust  the  predicted  change  (2.29) 
in  bed  elevation  as 

A Abk  = {\-m)/±Abk  (2.33) 

where  ra(0<w<l)isa  user-specified  modifier.  If  m = 0 , no  control, 
CONCEPTS  does  not  adjust  the  predicted  change  in  bed  elevation.  If  m — 1 , 
full  control,  CONCEPTS  keeps  the  bed  elevation  at  the  outlet  constant 
throughout  the  simulation. 

Internal  Boundary  Conditions 

At  hydraulic  structures  the  sediment  transport  rate  is  a function  of  the  supply  of 
sediment  and  the  elevation  of  the  upstream  invert  of  the  structure.  If  the 
upstream  invert  of  the  structure  is  located  above  the  channel  thalweg, 
CONCEPTS  deposits  that  part  of  the  sediment  load  transported  as  bed  load  just 
upstream  of  the  structure.  Based  on  the  following  criterion: 

n* 

— < 0.47  (2.34) 

CONCEPTS  determines  the  size  classes  k that  are  carried  mainly  as  bed  load 
and  are  deposited.  The  sediment  particles  in  the  other  size  classes  pass  through 
the  structure  without  depositing.  This  process  stops  once  the  thalweg  elevation 
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reaches  the  elevation  of  the  upstream  invert  of  the  structure.  At  this  point  all 
sediment  is  passed  through  the  structure  without  depositing. 


. 
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Bank  Erosion  and  Channel  Widening 


This  chapter  presents  the  stream-width  adjustment  submodel  of 
CONCEPTS- 1.0.  The  conceptual  basis  for  the  mathematical 
characterization  of  streambank-erosion  processes  are  discussed, 
followed  by  the  numerical  implementation  of  the  mathematical 
model. 


Theory 

Channel-width  adjustment  occurs  in  a wide  variety  of  geomorphic  contexts 
and  is  usually  accompanied  by  changes  in  other  morphological  parameters 
such  as  channel  depth,  roughness,  bed-material  composition,  riparian 
vegetation,  energy  slope,  and  channel  planform.  The  processes  responsible  for 
width  adjustment  are  diverse,  and  the  adjustment  process  itself  displays  a wide 
variety  of  spatial  and  temporal  patterns.  The  ASCE  Task  Committee  on 
Hydraulics,  Bank  Mechanics,  and  Modeling  of  River  Width  Adjustment  lists 
(ASCE,  1998): 

• widening  by  erosion  of  one  or  both  banks  without  substantial  incision; 

• widening  in  sinuous  channels  by  erosion  of  the  outer  bank  exceeding 
advancement  of  opposite  bank; 

• widening  in  braided  rivers  by  flows  deflected  around  growing  braid  bars; 

• rapid  widening  in  degrading  streams  by  increasing  bank  height  and 
steepness; 

• narrowing  by  formation  of  in-channel  berms  or  benches  at  the  margins; 

• narrowing  in  sinuous  channels  by  advancing  point  bars  exceeding  the  retreat 
of  the  opposite  cut  bank; 

• narrowing  in  braided  channels  by  the  abandoning  of  a marginal  anabranch. 
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It  is  unlikely  that  equilibrium  approaches  such  as  regime  theory,  extremal 
hypothesis,  or  tractive  force  methods  can  accurately  predict  width  adjustment 
over  time.  CONCEPTS  simulates  channel  width  adjustment  by  incorporating 
the  fundamental  physical  processes  responsible  for  bank  retreat: 

• fluvial  erosion  or  entrainment  of  bank  material  particles  by  the  flow,  and 

• mass  bank  failure,  for  example  due  to  channel  incision. 

Fluvial  Erosion  Process 

Water  flowing  in  an  alluvial  channel  exerts  forces  of  drag  and  lift  on  the 
boundaries  that  tend  to  detach  and  entrain  surface  particles.  The  origin  of 
resisting  forces  varies  according  to  the  grain  size  and  the  nature  of 
electrochemical  bonding  that  may  exist  between  cohesive  particles.  Alluvial 
bank  materials  are  formed  primarily  by  fluvial  deposition  and  are  often 
stratified.  Therefore,  the  characteristics  and  erodibility  of  the  bank  may  vary 
with  elevation.  Also,  floodplain  deposits  typically  include  alluvial  sand  and 
gravels,  clay  plugs,  and  strongly  cohesive  backswamp  deposits,  so  that  bank 
material  properties  vary  spatially  over  relatively  short  distances. 

In  the  case  of  noncohesive  sands  and  gravels,  the  forces  resisting  erosion  are 
generated  mainly  by  the  immersed  weight  of  the  particles.  Generally,  the 
mobility  of  noncohesive  bank  materials  can  be  predicted  using  a Shields-type 
entrainment  function,  but  this  must  be  modified  to  take  into  account  the 
destabilizing  effect  of  the  channel  side  slope. 

Fine-grained  bank  materials,  containing  significant  amounts  of  silt  and  clay, 
are  to  some  degree  cohesive  and  resist  entrainment  primarily  through 
interparticle,  electrochemical  bonding  rather  than  through  the  immersed 
weight  of  the  particles.  When  cohesive  bank  materials  are  entrained  by  the 
flow,  it  is  aggregates  of  grains  (such  as  soil  crumbs  or  peds  that  have  been 
produced  by  soil-forming  processes)  that  are  detached.  Fluvial  entrainment, 
therefore,  requires  that  the  local  boundary  shear  stresses  exceed  the  critical 
value  to  initiate  motion  of  aggregates  or  peds  rather  than  that  related  to  the 
primary  soil  particles.  Ped  size,  stability,  and  interped  bonding  strength  are  not 
conservative  soil  properties  as  they  depend  to  some  degree  on  the  local  history 
of  soil  development,  in  general,  and  recent  antecedent  conditions  of  wetting 
and  drying,  in  particular.  It  follows  that  the  conditions  of  incipient  motion  for 
cohesive  bank  materials  are  complex,  time-dependent,  and  difficult  to  define. 
Critical  boundary  shear  stresses  for  cohesive  bank  soils  tend  to  be  higher  than 
for  noncohesive  bank  materials.  As  a result,  erosion  rates  for  cohesive  banks 
are  generally  lower  than  for  noncohesive  banks.  Once  entrained,  aggregates 
disintegrate  rapidly  due  to  corrasion  at  the  channel  boundaries  and  turbulent 
buffeting  in  the  flow,  so  that  most  fine  sediment  derived  from  bank  erosion  is 
transported  in  suspension  and  is  conventionally  classified  as  wash  load. 

The  erodibility  of  bank  soils  can  be  increased  markedly  by  weakening 
processes  such  as  weathering.  The  processes  responsible  for  loosening  and 

52  CONCEPTS  Manual 


Theory 


detaching  grains  and  aggregates  are  closely  associated  with  soil  moisture 
conditions  at  and  beneath  the  bank  surface.  Swelling  and  shrinking  of  soils 
during  repeated  cycles  of  wetting  and  drying  can  contribute  to  cracking  that 
significantly  increases  erodibility  and  reduces  soil  shear  strength.  Heaving  and 
the  growth  of  needle  ice  crystals  at  the  bank  surface,  followed  by  collapse  of 
ice  wedges  and  needles  during  thawing  of  soil  moisture,  are  highly  effective  in 
increasing  the  susceptibility  of  cohesive  bank  materials  to  fluvial  erosion. 

Temporal  variability  in  the  erodibility  of  bank  soils  due  to  the  operation  of 
weathering  processes  means  that  the  effectiveness  of  a given  flow  event  in 
eroding  the  bank  depends  not  only  on  the  magnitude  and  duration  of  a 
particular  event  but  also  on  the  antecedent  conditions. 

Bank  vegetation  locally  increases  the  roughness  of  the  boundary,  increasing 
flow  resistance,  retarding  the  near-bank  flow,  and  damping  turbulence  (see, 
among  others,  the  early  work  by  Kouwen  et  al. , 1969,  and  the  more  recent 
work  by  Lopez  and  Garcia,  1997).  The  type  and  density  of  vegetation  are  very 
important.  Flexible  vegetation,  such  as  grasses  and  shrubs , are  usually  only 
effective  in  retarding  the  flow  at  low  velocities.  Rigid,  woody  species  continue 
retarding  the  flow  up  to  very  high  velocities,  but  may  generate  serious  bank 
scour  through  local  acceleration  of  flow  around  their  trunks.  For  trees  to  be 
effective  in  reducing  flow  attack  on  the  bank  they  must  be  spaced  sufficiently 
closely  (Thome,  1990). 

Vegetation  can  also  increase  the  erosion  resistance  of  the  soil.  The  roots  and 
rhizomes  of  plants  bind  the  soil  and  introduce  an  ‘added’  cohesion  over  the 
intrinsic  cohesion  of  the  bank  material.  The  critical  condition  for  erosion  of  a 
vegetated  bank  is  the  threshold  of  failure  of  the  plant  stems  by  snapping,  stem 
scour,  or  uprooting  rather  than  that  for  detachment  of  bank  material  itself.  This 
is  usually  associated  with  much  higher  levels  of  flow  intensity.  Further, 
vegetated  banks  are  generally  better  drained  and  drier  than  unvegetated  banks, 
so  that  the  impact  of  moisture-related  processes  that  weaken  and  loosen  the 
soil  is  reduced.  Gray  and  Leiser  (1982)  and  Coppin  and  Richards  (1990)  have 
reviewed  the  effects  of  herbaceous  vegetation  in  reducing  flow  erosivity  and 
bank  erodibility  and  concluded  that  major  effects  include  the  following: 

• Foliage  and  plant  residues  intercept  and  absorb  rainfall  energy  and  prevent 
soil  compaction  by  raindrop  impact. 

• Root  systems  physically  retain  soil  particles. 

• Near-bank  velocities  are  retarded  by  increased  roughness. 

• Plant  stems  dampen  turbulence  to  reduce  instantaneous  peak  shear  stresses. 

• Roots  and  humus  increase  permeability  and  reduce  excess  pore-water 
pressures. 

• Depletion  of  soil  moisture  reduces  water-logging. 
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Figure  3.1  Critical  shear  stress  xe  versus  SAR  for  different  soil  salt  concentrations 
and  dielectric  dispersion  values  (modified  from  Arulanandan  et  al.,  1980). 


Osman  and  Thorne  Method 

Osman  and  Thome  (1988)  developed  a method  to  predict  the  lateral-erosion 
rate  at  the  bank  toe.  Osman  and  Thome  follow  the  excess  shear  stress  approach 
set  forth  by  Arulanandan  et  al.  (1980).  Excess  shear  stress  is  defined  as  the 
difference  between  the  shear  stress  exerted  by  the  flowing  water  on  the  bank,  x, 
and  the  critical  shear  stress,  xe,  at  which  the  soil  particles  are  entrained. 

Arulanandan  et  al.  (1980)  determined  the  critical  shear  stress  as  a function  of 
sodium  adsorption  ratio  (SAR),  pore-fluid  salt  concentration,  and  the  dielectric 
dispersion  as  shown  in  Figure  3.1.  Once  you  obtain  the  value  of  xe,  the  method 

consists  of  the  following  steps  if  x > xe : 
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Step  1 . An  initial  rate  of  lateral  bank  erosion  Wl  is  given  by 

w,  = xp(_,  3 ) (3.1) 

Pi 

where  p5  is  mass  density  of  the  bank  material. 

Step  2.  The  rate  of  soil  erosion  is  assumed  to  have  an  approximately  linear 

increase  with  shear  stress  once  the  critical  shear  stress  is  exceeded.  The 

actual  erosion  rate  W is  then 

W = wl—-  0 (3.2) 

The  lateral  erosion  distance  during  a time  step  At  is 

A W = WAt  (3.3) 

A lower  bound  of  xe  for  the  application  of  Osman  and  Thome’s  method  is 
0.6  Pa.  Soils  with  lower  critical  shear  stresses  are  only  weakly  cohesive  and 
will  generally  behave  as  cohesionless  materials. 

Mass  Wasting  Process 

Fluvial  erosion  drives  bank  retreat  directly  by  the  removal  of  material  from  the 
bank  toe,  and  indirectly  causes  bank  retreat  by  triggering  mass  instability.  The 
stability  of  the  bank  with  regard  to  mass  failure  depends  on  the  balance 
between  gravitational  forces  that  tend  to  drive  the  soil  mass  downwards  and 
the  forces  of  friction  and  cohesion  that  resist  movement.  Failure  of  the  bank 
occurs  when  driving  forces  exceed  resisting  forces  on  the  most  critical 
potential  failure  surface , and  the  bank  collapses  in  a gravity-induced,  mass 
failure.  This  may  occur  as  fluvial  erosion  of  the  bed  next  to  the  bank  toe 
increases  the  bank  height,  or  undercutting  increases  the  bank  angle. 

The  analysis  of  slope  stability  with  respect  to  mass  failure  has  been  the  topic  of 
considerable  research  effort,  primarily  by  geotechnical  engineers,  but  also  by 
geomorphologists  and  geophysicists.  The  treatment  of  river  bank  stability  is 
generally  based  on  research  for  engineered  slopes  and  embankments. 

There  is  a clear  contrast  in  failure  mechanics  between  noncohesive  and 
cohesive  materials  because  of  significant  differences  in  their  soil  mechanics.  In 
a noncohesive  bank,  shear  strength  increases  more  rapidly  with  depth  than 
does  shear  stress,  so  that  critical  conditions  are  more  likely  to  occur  at  shallow 
depths.  In  a cohesive  bank,  however,  shear  stress  increases  more  quickly  than 
shear  strength  with  increasing  depth  so  that  critical  surfaces  tend  to  be  located 
deep  within  the  bank. 

Noncohesive  materials  usually  fail  by  dislodgement  and  avalanching  of 
individual  particles  or  by  shear  failure  along  shallow,  very  slightly  curved  slip 
surfaces.  Deep-seated  failures  occur  in  cohesive  materials  with  a block  of 
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disturbed,  but  more  or  less  intact,  bank  material  sliding  into  the  channel  along 
a curved  failure  surface.  In  banks  with  shallow  slope  angles  (<60°),  the  failure 
surface  is  curved  and  the  block  tends  to  rotate  back  toward  the  bank  as  it  slides, 
in  a rotational  slip  (Figure  3.2a).  Steep  banks  characteristically  fail  along 
almost  planar  surfaces , with  the  detached  block  of  soil  sliding  downward  and 
outward  into  the  channel  in  either  a planar  slip  or  a toppling  failure  (Figure 
3.2b). 

Cantilevered  or  overhanging  banks  are  generated  when  erosion  of  an  erodible 
layer  in  a stratified  or  composite  bank  leads  to  undermining  of  overlying, 
erosion-resistant  layers  (Figure  3.2c).  Frequently,  the  strength  of  cantilever 
blocks  is  significantly  increased  by  root  reinforcement  due  to  riparian  and 
floodplain  vegetation  and  decreased  by  tension  cracks. 

Tensile  forces  (negative  earth  pressure)  may  develop  in  the  upper  soil  layer 
producing  cracks  as  a consequence.  Terzaghi  (1943)  gives  the  thickness  of  this 
layer  as: 

z = — tan(45  + (f)/2)  (3.4) 

Is 

In  material  with  little  ability  to  withstand  tensile  stresses,  a small  cut  will  result 
in  time  in  the  development  of  a tension  crack  (with  height  ztc)  deep  enough  to 

produce  instability.  Commonly,  assuming  ztc  = 0 may  lead  to  better 

agreement  between  observed  and  predicted  failures  (Lohnes  and  Handy,  1 968). 
Tension  cracks  may  develop  at  the  instant  bank  failure  is  about  to  occur, 
therefore  not  significantly  contributing  to  the  instability  of  the  bank. 

Hagerty  (1990)  describes  streambank  collapse  by  exfiltrating  seepage,  which  is 
called  piping  or  sapping  (Figure  3.2dl-3).  The  outflowing  water  removes  soil 
particles  in  the  exfiltration  zone.  Undercut  strata  located  above  the  zone  of  soil 
loss  become  unstable  and  collapse. 

Whether  bank  failure  occurs  by  rotational  slip,  toppling,  or  cantilever  collapse, 
the  primary  force  tending  to  move  the  failure  block  is  the  weight  of  the  block. 
Fluvial  erosion  can  increase  the  driving  force  by  increasing  the  bank  height  or 
by  increasing  the  bank  slope  angle.  The  weight  of  bank  material  also  increases 
with  the  moisture  content  of  the  soil  and  failure  often  follows  the  change  from 
submerged  to  saturated  conditions  that  develops  when  drawdown  occurs  in  the 
channel  (e.g.,  Sands  and  Kapitzke,  1998). 

The  role  of  vegetation  in  affecting  mass  bank  failure  and  width  adjustment  is 
complex  and  poorly  understood.  Although  vegetation  generally  reduces  soil 
erodibility,  its  impact  on  bank  stability  with  respect  to  mass  failure  may  be 
either  stabilizing  or  destabilizing.  Hence,  depending  on  the  geomorphic 
context  and  dominance  of  either  fluvial  processes  or  mass  failure,  vegetation 
may  produce  either  a net  increase  or  a decrease  in  bank  stability. 

Gray  and  Leiser  (1982)  and  Coppin  and  Richards  (1990)  reviewed  the  ways 
that  woody  vegetation  may  affect  the  balance  of  forces  promoting  and  resisting 
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Figure  3.2  Bank  failure  mechanisms:  a)  rotational  failure,  b)  planar  failure,  c) 
cantilever  failure,  and  d)  piping  or  sapping  failure. 
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mass  failure.  Roots  mechanically  reinforce  soil  by  transferring  shear  stresses  in 
the  soil  to  tensile  stresses  in  the  roots,  which  root  strength  is  able  to  resist. 
However,  this  effect  operates  only  to  the  rooting  depth  of  the  vegetation,  and  it 
does  not  reinforce  potential  failure  planes  that  pass  beneath  the  plant  rootballs. 
Hence,  root  reinforcement  is  negated  when  bank  height  significantly  exceeds 
rooting  depth. 

Soil  moisture  levels  are  decreased  by  canopy  interception  and 
evapotranspiration,  reducing  the  frequency  of  occurrence  of  saturated 
conditions  conducive  to  bank  collapse,  or  increased  by  stemflow  (Simon  and 
Collison,  2001).  Anchored  and  embedded  stems  can  act  as  buttress  piles  or 
arch  abutments  in  a slope,  counteracting  downslope  shear  stresses  and 
increasing  bank  stability.  However,  roots  may  also  invade  cracks  and  fissures 
in  a soil  or  rock  mass  and  thereby  cause  local  instability  by  their  wedging  or 
prying  action.  The  surcharge  of  weight  of  vegetation  may  significantly 
increase  motivating  forces,  causing  destabilization  of  the  bank,  and  wind 
loading  on  tall  vegetation  may  exert  an  additional  and  potentially  critical 
destabilizing  moment  on  the  bank. 

Streambank  Stability  Analysis 

Streambank  stability  can  be  analyzed  using  methods  developed  for  engineered 
slopes  and  embankments,  e.g.,  Bishop  (1955),  Morgenstem  and  Price  (1965), 
Terzaghi  and  Peck  (1967),  and  Fredlund  and  Krahn  (1977).  These  limit 
equilibrium  methods  are  based  upon  static  equilibriums  of  forces  and/or 
moments.  Simon  et  al.  (1999)  have  presented  a comprehensive  stability 
analysis  for  layered  streambanks  based  on  the  limit  equilibrium  method.  They 
included  forces  due  to  pore-water  pressures  and  confining  pressures.  We  have 
extended  their  stability  analysis  by  explicitly  accounting  for  the  seepage  force 
in  the  factor  of  safety  and  implementing  an  alternative  formulation  of 
confining  pressures,  though  presently  for  homogeneous  bank  material. 

Following  Huang  (1983)  the  surface  water  on  the  failure  block  is  modeled  by 
assuming  it  is  a material  with  no  strength  (cohesion  and  angle  of  internal 
friction  are  set  to  zero).  Hence,  the  slip  surface  is  extended  vertically  through 
the  water,  and  a horizontal  hydrostatic  force  is  applied  on  the  vertical  portion 
of  the  slip  surface  (see  Figure  3.3). 

The  forces  acting  on  the  failure  block  are  (Figure  3.3): 

1 the  weight  of  the  failure  block , Ws; 

2 the  weight  of  surface  water  on  the  failure  block,  Ww; 

3 the  hydrostatic  force  exerted  by  surface  water  on  the  vertical  slip  surface, 

F ■ 

4 the  hydrostatic  force  exerted  by  water  in  the  tension  crack,  Ft; 

5 the  seepage  force,  Fs\ 

6 the  shear  force  at  the  base  of  the  failure  block,  S\  and 


u a 1 
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Figure  3.3  Definition  sketch  of  forces  acting  on  a streambank. 
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7 the  total  normal  force  at  the  base  of  the  failure  block,  N. 

The  shear  force  acting  at  the  base  of  the  failure  block  can  be  written  using  the 
shear  strength  equation  for  either  saturated  or  unsaturated  soils  (Fredlund  and 
Rahardjo,  1993): 

5=  ~( c'  + ( <j  n — u q)  tan  ({>'  + ( ua  — uw)  tan  (f>^)  (3.5) 

b 

where  cr,;  is  total  stress  normal  to  the  base  of  the  failure  block,  ua  is  pore-air 
pressure,  uw  is  pore-water  pressure,  c'  is  the  effective  cohesion,  (J)'  is  the 
effective  angle  of  internal  friction  associated  with  the  normal  stress  state 

variable  (a/;  - uci) , <\>h  is  an  angle  indicating  the  increase  in  shear  strength  for 

an  increase  in  matric  suction  (ua  - uw) , L is  the  length  of  the  inclined  slip 

surface,  and  F is  the  factor  of  safety  which  is  defined  as  the  factor  by  which  the 
shear  strength  parameters  must  be  reduced  in  order  to  bring  the  soil  mass  into  a 
state  of  limiting  equilibrium  along  the  assumed  slip  surface. 

Integrating  the  expression  of  static  equilibrium  of  a two-dimensional  soil 
element 

V - a + Jbg  = 0 (3-6) 

where  a is  the  total  stress  tensor,  yb  is  bulk  unit  weight  of  the  bank  material, 
and  g is  the  gravitational  acceleration  vector,  over  the  failure  block  yields  the 
following  force  balances  in  horizontal  and  vertical  directions: 

Fw-Ft  + ScosP  - NsinP  = 0 (3.7a) 

Ws  + Ww  - Ssin (3  - TVcos P = 0 (3.7b) 


59 


Bank  Erosion  and  Channel  Widening 


where  (3  is  the  failure-plane  angle.  Setting  pore-air  pressures  to  zero 
(approximately  true  in  natural  soil  deposits),  and  using  an  = N/ L and 

i'w  = Uw/L  ’ where  Uw  is  pore-water  force  acting  on  the  slip  surface,  the 
shear  force  becomes 

S = ~(Lc'  + Man(|)'  + U tan<j>*)  (3.8) 


Eliminating  the  shear  force  S and  nonnal  force  N from  (3.7a),  (3.7b),  and  (3.8) 
yields  the  following  expression  for  the  factor  of  safety : 

_ c'L  + [(Ws  + lEw,)cos(3  + (Fw-F?)sinP]tan(j)'  - £/wtan<|>/> 


(Ws+  Ww) sin p - (Fw -Ft) cos p 


(3.9) 


For  saturated  soils  <j>  = (j)' , and  (3.9)  simplifies  to 

F _ c'L+  [(Ws+  JVw)  cos  (3  - Uw  + ( Fw  - F, ) sin  [3  ] tan  4>* 
(Ws+  Ww) sin P - (Fw - F t) cos P 


Pore-Water  Pressure 

The  pore-water  pressure  distribution  within  a saturated  zone  of  a two- 
dimensional,  semi-infinite  soil  mass  is  hydrostatic  if  the  groundwater  is  at  rest 
and  the  water  density  is  homogeneous.  We  can  readily  calculate  pore-water 
pressures  if  we  know  the  position  of  the  phreatic  surface.  If  there  is 
groundwater  flow  the  pore-water  pressure  distribution  is  no  longer  hydrostatic 
because  pressure  gradients  are  needed  to  maintain  the  flow.  In  the  following 
we  assume  a steady,  uniform  Darcian  flow  field  throughout  the  streambank. 
This  is  a considerable  approximation  because  of  the  curvature  and  fluctuation 
of  the  phreatic  surface  near  the  bank  surface  (Figure  3.3).  However,  it  allows 
us  to  examine  the  effect  of  seepage  on  the  factor  of  safety  of  a streambank  in  a 
simplified  manner. 

The  seepage  force  Fs  per  unit  area  is  defined  as 

fs  = <>»■  i3-11) 

where  i is  hydraulic  gradient  with  i = J(dh/dx)2  + ( dh/dzf , i.e.  i is  equal  to 
the  magnitude  of  the  change  in  head  per  unit  length  along  a flow  line,  h is  total 
hydraulic  head,  and  yw  is  unit  weight  of  water.  The  angle  a (see  Figure  3.3) 
defines  the  direction  of  the  seepage  force.  The  horizontal  component  of  the 
seepage  force  is  Fs  = |FS|  cos  a and  the  vertical  component  of  the  seepage 

force  is  Fs  = -\FS\  sin  a . Equilibrium  of  the  groundwater  requires 

du 

~L  + FS  = 0 (3.12a) 

dr 
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duU! 

— - + y + F = 0 
dz  r w s-- 


(3.12b) 


The  pore-water  distribution  is  linear  in  x and  z for  the  case  of  uniform,  Darcian 
seepage: 


uw  ~ c\x+  C2z+  C3 


(3.13) 


where  C1?  C2,  and  C3  are  constants.  From  (3.12a)  and  (3.12b)  it  follows  that 
C]  = -F  and  C9  = -(yu,  + Fs  ) . We  can  determine  C3  from  the  position 
(xp7zp)  (cf.  Figure  3.3)  of  the  phreatic  surface  on  the  face  of  the  streambank 
where  the  pore- water  pressure  equals  zero:  C3  = Fs  xp  + (Y  w + Fs)Zp-The 
pore-water  distribution  then  reads 

uw  = Yw(zp-z)  + YwliK(x/7-x)cosa-(zp-z)sina)  (3-14) 


where  the  first  term  on  the  right-hand  side  of  (3.14)  is  the  hydrostatic  part  of 
the  distribution  and  the  second  term  on  the  right-hand  side  of  (3.14)  is  the  non- 
hydrostatic part  of  the  distribution  due  to  seepage.  Equating  (3.14)  to  zero 
gives  the  position  of  the  phreatic  surface 


/cos a / x 

= zp+  i i'-i"1 — (XP~X) 
p 1 - / sina  1 


(3.15) 


The  coordinate  pair  (xp^p)  can  be  selected  as  the  top  of  the  seepage  face  on  the 
streambank. 


Factor  of  Safety  without  Seepage 

Herein,  we  derive  an  expression  for  the  factor  of  safety  for  a streambank 
assuming  the  groundwater  is  at  rest.  Hence,  the  first  term  on  the  right-hand 
side  of  (3.14)  represents  the  pore- water  pressure  distribution.  Commonly,  the 
pore-water  pressure  ratio  ru  is  used  instead  of  pore-water  pressure  Uw.  Bishop 
and  Morgenstem  (1960)  defined  the  pore-water  pressure  ratio  as 


ru  = 


(3.16) 


where  d is  depth  below  the  soil  surface.  Integrating  (3.16)  we  have 
Uw  = ru  W s sec  P > an<3  consequently 

^ c'L  + [(( 1 - r(/sec“P)  W s + fTu,)cos(3  + (Fw  - F^sinP]  tan ({>’  + ^tanij)^ 
F ( + Ww)  sin P - (Fw  - F f)  cos  P 

(3.17) 

where  T1  is  the  total  negative  pore-water  pressure  or  matric  suction  along  the 
slip  surface  in  the  unsaturated  zone.  Huang  (1983)  argued  that  definition  (3.16) 

of  pore- water  pressure  ratio  is  not  reasonable  as  the  term  1 - /-;(sec~P  becomes 

negative  for  larger  failure  plane  angles.  The  pore-water  pressure  ratio  ru  ~ 0.5 
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for  a fully  saturated  failure  block;  if  (3  > 45°  the  effective  force  normal  to  the 
failure  plane  is 

N = (1  -r((sec2p)^cosp<0  (3.18) 

Huang  (1983)  and  others  (e.g.,  Fredlund,  1995)  use  the  concept  of  buoyant 
weight  to  overcome  this  problem.  The  pore-water  pressure  is  used  to  reduce  the 
effective  stress  along  the  failure  surface  by  a factor  ( 1 - ru) , so 

N = ( 1 - ru)  ff^cosp . Huang  (1983)  defined  the  pore-water  pressure  ratio  as 

the  ratio  between  the  total  pore-water  pressure  and  the  total  overburden 
pressure: 

volume  of  sliding  mass  under  water  x yyA> 

r = — (3.19) 

volume  of  sliding  mass  x yb 

We  adjust  (3.19)  to  account  for  the  presence  of  surface  water  on  the  bank  as 
follows: 

(vol.  sliding  mass  under  water  + vol.  surface  water  on  bank)  x y 

r = (3.20) 

vol.  sliding  mass  xyb  + vol.  surface  water  on  bank  x y 


where  the  text  “vol.”  means  “volume  of.”  Introducing  (3.20)  into  (3.9)  yields 

c'L+[{  1 -r  )(W  + W w)cosf3  + (FM;  — Fisinpltanf  + ^ tanij)^ 

F = (3.21) 

( fVs  + Ww)  sin  P - ( Fw  - F t)  cos  P 

Defining  apparent  or  total  cohesion  as  (Fredlund  and  Rahardjo,  1993) 

ca  = c'  - ~(ru(  Ws  + lTvv)cosptan(j)'  - THan^)  (3.22) 

L 


the  factor  of  safety  is 

F = caL  + [{Ws  + Ww ) cos  P + ( Fw  - Fr)  sin  P ] tan  <j)' 
(JVS+  Ww)  sinp  - (Fw  -Ft)  cos  P 


(3.23) 


Factor  of  safety  with  Seepage 

Seepage  affects  the  factor  of  safety  through  apparent  cohesion.  The  horizontal 
component  of  the  seepage  force,  Fs  , influences  the  pore-water  pressure 

distribution  by  tilting  the  phreatic  surface,  consequently  changing  the  saturated 
area  of  the  failure  block.  The  vertical  component  of  the  seepage  force,  F , 

changes  the  vertical  gradient  of  pore-water  pressure.  We  can  account  for  the 
seepage  force  by  adjusting  ru,  T1,  and  introducing  an  excess  positive  pore- 
water  force  in  (3.21): 

Cq  = c'--[(r;(  + f;()(lTs,+ !TJcosptan4),  + t7M,tan(}),-(vT  + T>)tan(t)6] 

Z; 

(3.24) 

where  the  overbar  denotes  a correction  of  the  variable  due  to  seepage. 
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Factor  of  safety  for  Layered  Streambanks 

We  can  also  use  (3.23)  to  determine  the  stability  of  layered  streambanks  by 
using  a weighted,  average  cohesion  and  angle  of  internal  friction  for  the 
heterogeneous  bank  material.  However,  we  cannot  readily  extend  (3.23)  to 
account  for  the  strength  of  each  soil  in  a layered  streambank.  In  this  case  we 
have  to  revert  to  the  method  of  slices  (Fellenius,  1936). 

Inclination  of  the  Failure  Surface 

Darby  and  Thome  (1996)  followed  Taylor  (1948),  who  stated  that  the  most 
dangerous  failure  surface  is  that  with  the  maximum  stability  number 
Nb  = y hH/cd , where  H is  height  of  failure  block  and  cd  is  developed 
cohesion,  or  equivalently  the  failure  surface  for  which  cohesion  is  fully 
mobilized,  that  is,  dcd/d$  = 0 . However,  Darby  and  Thome  (1996)  assumed 

cd  = c' . Effective  cohesion  is  primarily  determined  by  soil  structure  and  soil 

moisture,  and  therefore  may  not  vary  much  with  failure  plane  angle.  Hence, 
Darby  and  Thome’s  approach  of  rewriting  (3.21)  to  yield  an  expression  for 
effective  cohesion,  and  equating  the  derivative  of  this  expression  with  respect 
to  failure  plane  angle  to  zero,  is  unrealistic.  Moreover,  the  factor  of  safety  for 
the  resulting  failure  plane  angle  is  not  a minimum. 

Here,  the  inclination  of  the  failure  plane  is  that  for  which  the  factor  of  safety  is 
minimum,  that  is,  5F/d(3  = 0 . Taking  the  derivative  of  (3.23)  with  respect  to 
the  failure  plane  angle  gives 


(3.25) 


/[(Ws+  Ww)  sin  (3  - (Fw  - Ffcos^f 


Equating  (3.25)  to  zero  yields  the  failure  plane  angle. 
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Figure  3.4  Representation  of  fluvial  erosion.  The  highlighted  area  is  scoured  by  the 
flow,  steepening  the  bank. 


Implementation  of  the  Channel-Widening  Submodel 

Fluvial  Erosion 


CONCEPTS  uses  (3. 1 ) — (3.3)  to  compute  lateral  erosion  distance  at  the  toe  of 
the  bank,  which  is  used  to  adjust  the  bank  profile  as  shown  in  Figure  3.4.  Like 
Osman  and  Thome  (1988),  CONCEPTS  assumes  that  lateral  erosion  produces 
a vertical  cut  on  the  streambank  and  a horizontal  cut  on  the  streambed. 
CONCEPTS  calculates  the  eroded  volume  per  unit  channel  length,  Vw,  by 
determining  the  intersection  of  the  vertical  cut  and  bank  profile  and  using  AW. 
We  convert  the  eroded  volume  of  bank  material  to  a lateral  flux  of  sediment 
into  the  channel 


nLi 

At 


(3.26) 


in  which  is  fractional  content  by  volume  of  size  class  k in  the  streambank. 

Equation  (3.26)  shows  that  the  eroded  volume  is  uniformly  partitioned  over  the 
time  step  At.  The  lateral  flux  is  added  to  the  conservation  of  sediment  mass 
carried  by  the  flow  (2.1). 


Mass  Wasting 

The  slip  surface  does  not  need  to  intersect  the  streambank  at  the  toe  of  the 
bank.  The  bank-stability  algorithm  evaluates  the  factor  of  safety  at  Ne 
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Figure  3.5  Possible  failure  block  geometries  for  the  case  where  the  number  of 
elevations  at  which  the  factor  of  safety  is  evaluated,  N?,  equals  four. 
Sample  values  of  factor  of  safety  are  given  for  each  failure  block 
geometry.  Here,  factor  of  safety  is  smallest  for  block  #2,  and  smaller  than 
unity.  Hence,  the  streambank  fails  along  the  respective  slip  surface. 


elevations,  which  are  equidistantly  distributed  along  the  bank  profile.  The  user 
must  specify  Ne  as  an  input  parameter  to  the  model. 

The  bank-stability  analysis  is  a two-step  process.  The  algorithm  firstly 
determines  the  angle  of  the  slip  surface  that  minimizes  factor  of  safety  at  each 
of  the  Ne  elevations  (Figure  3.5),  and  consequently  computes  factor  of  safety. 
Finally,  the  algorithm  selects  the  smallest  of  the  Ne  number  of  factor  of 
safeties.  The  bank  fails  if  this  factor  of  safety  is  smaller  than  unity. 

If  bank  failure  occurs,  CONCEPTS  computes  the  volume  of  the  failure  block 
and  converts  it  to  a lateral  flux  of  sediment  into  the  channel  (Simon  et  al., 
1991)  as  in  (3.26).  Field  studies  (Simon  et  al.,  1999)  show  that  this  is  not  very 
realistic;  the  failure  block  will  break  into  smaller  blocks  and  aggregates  in 
which  cohesion  holds  sediment  particles  together.  A deterministic  model  of 
this  process,  however  is  not  yet  available. 

The  user  specifies  the  depth  of  the  tension  crack  as  a fraction  ft  of  the  depth  of 
the  soil  layer  with  negative  earth  pressure  (3.4)  as: 

= U,  (3-271 
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Data 

Specification 


Input  Data 

This  chapter  presents  the  input  data  structure  of  CONCEPTS. 
CONCEPTS  requires  the  following  input  files  created  by  the  user: 
(1)  one  input  file  with  run  control  data,  (2)  one  input  file  with  dis- 
charge data  at  the  upstream  boundary  of  the  modeling  reach,  (3)  an 
input  file  for  each  cross  section  in  the  modeling  reach,  and  (4)  an 
input  file  for  each  hydraulic  structure  in  the  modeling  reach. 

Run  Control  Data 


The  run  control  data  input  file  consists  of  three  main  blocks  of  data: 

Block  1 . project  related  data,  general  data,  and  makeup  of  modeling  reach 
Block  2.  names  of  cross  section  and  hydraulic  structure  files 
Block  3.  output  options 

Block  1 


Block  1 consists  of  31  + input  lines,  see  Table  4.1,  where  Nj  is  the  number 

of  links  comprising  the  simulated  stream  corridor.  A link  is  defined  as  a 
subreach  or  a hydraulic  structure.  A subreach  is  a channel  segment  between 
two  structures.  If  there  are  no  hydraulic  structures  along  the  simulated  channel, 
the  number  of  links  equals  one  because  there  is  just  one  reach.  If  there  is  one 
hydraulic  structure  in  the  modeling  reach,  say  a box  culvert,  the  number  of 
links  equals  three:  (1)  a subreach  between  the  upstream  boundary  and  the  box 
culvert,  (2)  the  box  culvert  itself,  and  (3)  a subreach  between  the  box  culvert 
and  the  downstream  boundary. 

Line  5 contains  the  run  identifier,  which  is  a twelve-character  string. 
CONCEPTS  uses  the  run  identifier  to  name  the  output  files  as:  run  identifier  + 
‘ ’ + number  + ‘ .txt’.  For  example,  if  the  run  identifier  is  simulation,  then 
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Table  4.1  Input  data  in  block  #1  of  run  control  data. 


line  number 

line  description 

data  type 

1-4 

comment  lines 

5 

run  identifier 

string 

6 

comment  line 

7 

project  title 

string 

8-9 

comment  lines 

10 

name  of  discharge  file 

string 

11 

comment  line 

12 

flow-related  data 

real/integer 

13 

comment  line 

14 

upstream  sediment  discharge  option  and  rates 

integer/real 

15 

comment  line 

16 

fraction  of  fines  in  a cohesive  bed  and  bed 
control  at  downstream  boundary 

real 

17 

comment  line 

18 

bank-failure  analysis  options 

integer 

19 

comment  line 

20 

type  of  bed  resistance  formulation 

integer 

21 

comment  line 

22 

water  temperature 

real 

23 

comment  line 

24 

included  submodels 

integer 

25-26 

comment  lines 

27 

simulation  times 

string/integer 

28-29 

comment  lines 

30 

number  of  links  N, 

31 

comment  line 

32  - 32  + N,~  1 

link  IDs 

the  output  filenames  are  ‘simulation  050.txt’,  ‘simulation  051.txt’,  etc. 
The  numbering  of  the  output  files  starts  at  50.  The  project  title  on  line  7 is  an 
80-character  string. 

Line  10  contains  a 32-character  string  that  is  the  name  of  the  file  with  the  flow 
and  sediment  discharge  time-series  at  the  upstream  boundary  of  the  modeling 
reach  (see  page  88).  Line  12  contains:  (1)  baseflow  discharge  at  the  upstream 

boundary  (m3/s);  (2)  the  rate  of  lateral  inflow  ( Cq ) as  a fraction  of  the 

discharge  at  the  upstream  boundary  (nf1),  that  is 

q = CqQx  (4.1) 

where  Q j is  the  discharge  imposed  at  the  upstream  boundary;  and  (3)  a flag 
indicating  the  type  of  downstream  boundary  condition  used  in  the  flow 
simulation.  If  the  flag  is  zero,  CONCEPTS  uses  the  loop-rating  curve  (1.35).  If 
the  flag  is  set  to  one,  CONCEPTS  uses  the  user-specified  rating  curve  (1.34). 


70 


CONCEPTS  Manual 


Run  Control  Data 


Table  4.2  Values  assigned  to  the  different  processes  that  can  be 
accounted  for  by  the  bank  stability  analysis. 


process  value 

positive  pore-water  pressures  1 

confining  pressures  2 

matric  suction  4 

seepage  8 


In  the  latter  case  line  12  also  contains  the  number  of  segments  Nrs  and  values 
for  hh  ah  and  (3,  for  each  segment  by  segment. 

Line  14  contains  a flag,  which  can  be  set  to  0 or  1 , indicating  how  sediment 
load  is  imposed  at  the  upstream  boundary.  If  the  user  sets  the  flag  to  zero,  line 
14  also  contains  the  fraction  fk  of  the  local  sediment  transport  capacity  for  each 

of  the  13  size  classes  (see  (2.32)).  Thus,  if  fk  = 0. 1 for  a certain  size  class,  the 

flow  will  carry  a sediment  load  equalling  10%  of  the  local  sediment  transport 
capacity  for  that  size  class.  If  the  user  sets  the  flag  to  one,  the  sediment  load  is 
given  in  the  file  with  upstream  boundary  conditions  (line  10). 

The  first  number  on  line  16  informs  CONCEPTS  when  the  streambed  can  be 
assumed  cohesive  or  cohesionless.  If  the  fraction  of  fines  in  the  bed  is  larger 
than  this  input  entry,  the  bed  is  assumed  to  be  cohesive.  CONCEPTS  defines 
fines  as  clay  and  silt  particles,  that  is,  size  classes  1 through  3 (see  Table  2.1). 
The  second  number  is  the  modifier  m that  adjusts  the  predicted  change  in  bed 
elevation  at  the  downstream  boundary  of  the  simulated  stream  corridor  (see 
(2.33)). 

Line  18  contains  two  integers  used  by  the  bank-stability  analysis  algorithm. 
The  first  number  denotes  the  level  of  complexity  of  the  stability  analysis.  The 
user  specifies  which  physical  processes  need  to  be  taken  into  account.  The 
analysis  distinguishes  the  following  processes:  positive  pore-water  pressures, 
confining  pressures,  matric  suction  or  negative  pore-water  pressures,  and 
seepage.  Each  process  has  an  assigned  value,  see  Table  4.2.  The  user  specifies 
the  complexity  of  the  analysis  by  adding  these  values.  For  example,  if  the  user 
wants  to  account  for  positive  pore-water  pressures  (a  value  of  1 ) and  for 
confining  pressures  (a  value  of  2),  the  user  needs  to  enter  a value  of  3 in  the 
input  file.  If  the  user  enters  zero,  none  of  the  four  processes  will  be  taken  into 
account.  The  second  number  on  line  1 8 is  the  number  of  possible  elevations  at 
which  the  slip  surface  may  intersect  the  bank  profile  (see  section  on  the 
implementation  of  “mass  wasting”). 

Line  20  contains  a flag  that  sets  the  type  of  streambed  resistance  formulation. 
If  the  entered  number  is  two  (2),  CONCEPTS  uses  the  friction-factor  relation 
proposed  by  Karim  (1995).  Karim  obtained  satisfactory  results  within  the 
following  range  of  variables:  flow  depths  from  0.03  to  16.7  m,  velocities  from 
0.32  to  3.41  m/s,  energy  slopes  from  0.0000183  to  0.0243,  median  sediment 
sizes  from  0.08  to  28.6  mm,  and  gradation  coefficients  of  bed  sediments  from  1 


71 


Input  Data 


Table  4.3  Link  types  to  be  used  with  CONCEPTS. 


link  type 

ID 

channel  reach 

1 

pipe  culvert 

11 

box  culvert 

13 

bridge  crossing 

21 

drop  structure 

26 

generic  structure 

31 

to  2.  If  the  user  enters  a value  other  than  two,  CONCEPTS  uses  a constant 
Manning  n,  which  is  specified  in  the  cross  section  data  file. 

Line  22  contains  the  water  temperature  (°C),  which  CONCEPTS  uses  to 
compute  water  density  and  viscosity. 

Line  24  contains  three  flags,  which  can  be  set  at  either  0 or  1.  These  flags 
determine  whether  CONCEPTS  takes  into  account  sediment  routing  and  bed 
adjustment,  streambank  fluvial  erosion,  and  streambank  mass-wasting, 
respectively.  If  the  user  sets  a flag  to  zero  the  corresponding  process  is  not 
taken  into  account. 

Line  27  contains  the  start  time  and  end  time  of  the  simulation,  and  the  time  step 
(sec).  The  start  and  end  time  are  19-character  strings  formatted  as  “mm/dd/ 
yyyy  hh:mm:ss.”  For  example,  two  o’clock  in  the  afternoon  on  June  7,  1999  is 
represented  as:  06/07/1999  14:00:00. 

CONCEPTS  uses  the  integer  value  entered  for  the  time  step  as  initial  time  step. 
During  the  simulation  the  time  step  is  automatically  adjusted  based  on 
temporal  variations  in  discharge  and  flow  depth. 

Line  30  contains  the  number  of  links.  Lines  32  through  31  + Nt  consist  of  the 

link  IDs  for  each  link  in  the  modeling  reach.  Table  4.3  lists  the  available 
linktypes. 

Block  2 

Block  2 identifies  the  files  that  contain  input  data  for  cross  sections  and 
hydraulic  structures.  It  consists  of  3 + Nx  lines  for  each  link  in  the  modeling 
reach,  where  Nx  is  the  number  of  cross  sections  if  the  link  is  a reach  and 
N - 0 if  the  link  is  a hydraulic  structure,  see  Table  4.4. 

If  the  link  is  a reach,  line  3 contains  an  integer  representing  the  number  of  cross 
sections  in  that  reach.  If  the  link  is  a hydraulic  structure,  line  3 contains  a 32- 
character  string  representing  the  name  of  the  file  with  the  input  data  of  the 
structure. 
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Table  4.4  Input  data  in  block  #2  of  run  control  data.  The  user  has  to  repeat  this  block 
for  each  link. 


line  number 

line  description 

data  type 

1 -2 

comment  lines 

O 

number  of  cross  sections,  Nx , if  link  is  a subreach 

integer 

o 

name  of  input  file  if  link  is  a hydraulic  structure 

string 

4-3  + NX 

name  of  input  file  for  each  cross  section  if  link  is  a 
subreach 

string 

Table  4.5  Input  data  to  request  output  at  a certain  cross  section  and  for  a 
certain  runoff  event.  Lines  2 through  4 + Ns  are  repeated  NIc 


times. 

line  number 

line  description 

data  typ 

1 

number  of  locations  Nlc 

integer 

2 

type  of  outputted  data 

integer 

3 

location  reference 

integer 

4 

number  of  storm  events  Ns 

integer 

5-  5 + NS-  1 

dates  of  storm  events 

string 

If  the  link  is  a reach,  lines  4 through  3 + Nx  contain  the  names  (32-character 
string)  of  the  files  with  data  related  to  cross  sections. 

Block  3 


Block  3 comprises  output  options.  There  are  three  output  categories: 

1 output  at  a certain  location  and  for  a certain  runoff  event, 

2 time-series  output  at  a certain  location,  and 

3 output  for  a certain  runoff  event  along  a section  of  the  modeling  reach. 

The  user  needs  to  specify  the  options  following  the  above  sequence. 

Output  at  a Certain  Location  and  for  a Certain  Runoff  Event 

Table  4.5  lists  the  data  the  user  has  to  enter  to  request  output  from  the  first 
category.  Line  1 contains  the  number  of  locations  (N/c)  at  which  output  is 
requested,  followed  by  the  type  of  data  on  line  2.  The  latter  is  a summation  of 
values  assigned  to  various  variables  (see  Table  4.6).  For  example,  if  the  user 
wants  to  output  peak  discharge  and  stage,  and  sediment  yield  for  a certain 
runoff  event,  the  user  needs  to  enter  21  (=1+4+16). 

Line  3 contains  information  on  the  location  of  the  cross  section  within  the 
modeling  reach.  The  user  has  to  enter  the  link  number  and  in-link  cross  section 
number,  for  example  cross  section  4 in  link  1 . 
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Table  4.6  Parameters  that  can  be  output  for  a certain  runoff  event  and  cross 


section. 

outputted  parameter  value 

peak  discharge  1 

peak  flow  depth  2 

peak  stage  4 

peak  friction  slope  8 

sediment  yield  16 

cumulative  sediment  yield  over  all  runoff  events  thus  far  32 

change  in  bed  elevation  64 

cumulative  change  in  bed  elevation  over  all  runoff  events  128 

thus  far 

lateral  erosion  256 

cumulative  lateral  erosion  over  all  runoff  events  thus  far  512 

cross-sectional  geometry  1 ,024 

in-bank  top  and  bottom  width  of  cross  section  2,048 

bank  height  4,096 

not  used 

characteristic  particle  sizes  16,384 

particle  size  distribution  32,768 


Table  4.7  Input  data  to  request  time-series  output  at  a certain  cross  section. 
Lines  2 through  4 + Nts  are  repeated  Nlc  times. 


line  number 

line  description 

data  type 

1 

number  of  locations  Nlc 

integer 

2 

type  of  outputted  data 

integer 

3 

location  reference 

integer 

4 

number  of  time  series  Nts 

integer 

5 - 5 + Nts  - 1 

start  and  end  times  of  each  series 

string 

Line  4 contains  the  number  of  runoff  events  (Ns)  for  which  output  is  requested 
at  that  particular  cross  section.  Lines  5 through  Ns  + 4 contain  the  dates  of 
occurrence  of  the  runoff  events. 

Lines  2 through  Ns  + 4 are  repeated  N/c  times. 

Time-Series  Output  at  a Certain  Location 

Table  4.7  lists  the  data  the  user  has  to  enter  to  request  time-series  output.  Line 
1 contains  the  number  of  locations  (N/c)  at  which  output  is  requested,  followed 
by  the  type  of  data  on  line  2.  The  latter  is  a summation  of  values  assigned  to 
various  variables  (see  Table  4.8).  For  example,  if  the  user  wants  to  output 
discharge,  stage,  and  factor  of  safety,  the  user  needs  to  enter  524,297 
(=1+8+524,288). 

Line  3 contains  information  on  the  location  of  the  cross  section  within  the 
modeling  reach.  The  user  has  to  enter  the  link  number  and  in-link  cross  section 
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Table  4.8  Parameters  that  can  be  output  as  time  series. 

outputted  parameter 

value 

discharge 

1 

velocity 

2 

flow  depth 

4 

stage 

8 

flow  area 

16 

flow  top  width 

32 

wetted  perimeter 

64 

hydraulic  radius 

128 

conveyance 

256 

friction  slope 

512 

energy  head 

1,024 

Froude  number 

2,048 

bed  shear  stress 

4,096 

sediment  discharge  (silt/sand/gravel/total) 

8,192 

cumulative  sediment  yield  (silt/sand/gravel/total) 

16,384 

cumulative  change  in  bed  elevation 

32,768 

thalweg  elevation 

65,536 

cumulative  lateral  erosion 
not  used 

131,072 

factor  of  safety 

524,288 

apparent  cohesion 

1,048,576 

pore-water  force 

2,097,152 

matric  suction  force 

4,194,304 

weight  of  failure  block 

8,388,608 

weight  of  water  on  the  bank 

16,777,216 

horizontal  component  of  the  confining  force 

33,554,432 

groundwater  elevation 

67,108,864 

location  of  bank  top 

134,217,728 

number,  for  example  cross  section  4 in  link  1 . Line  4 contains  the  number  of 
time  series  (Nts)  for  which  output  is  requested  at  that  particular  cross  section. 

Lines  5 through  Nts  + 4 contain  the  start  and  end  dates  of  the  time  series. 
Lines  2 through  N[s  + 4 are  repeated  N/c  times. 


Output  for  a Certain  Runoff  Event  along  a Section  of  the 
Modeling  Reach 

Table  4.9  lists  the  data  the  user  has  to  enter  to  request  output  for  a certain 
runoff  event  along  a section  of  the  modeling  reach,  hereafter  referred  to  as  a 
pro  file.  Line  1 contains  the  number  of  profiles  (Np)  at  which  output  is 
requested,  followed  by  the  type  of  data  on  line  2.  The  latter  is  a summation  of 
values  assigned  to  various  variables  (see  Table  4.10).  For  example,  if  the  user 
wants  to  output  peak  discharge  and  stage,  and  sediment  yield  for  a certain 
storm  event,  the  user  needs  to  enter  67  (=1+2+64). 
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Table  4.9  Input  data  to  request  output  for  a certain  runoff  event  along  a 

section  of  the  modeling  reach.  Lines  2 through  4 + Ns  are  repeated 


NP 

times. 

line  number 

line  description 

data  type 

1 

number  of  profiles  N 

integer 

2 

type  of  outputted  data 

integer 

3 

location  reference 

integer 

4 

number  of  storm  events  Ns 

integer 

5 - S + Ns-  1 

dates  of  storm  events 

string 

Table  4.10  Parameters  that  can  be  output  for  a certain  runoff  event  along  a 


section  of  the  modeling  reach, 
outputted  parameter  value 

peak  discharge  1 

peak  stage  2 

thalweg  elevation  4 

cumulative  change  in  bed  elevation  over  all  runoff  8 

events  thus  far 

in-bank  top  width  16 

bank  height  32 

sediment  yield  64 

characteristic  particle  sizes  128 


Line  3 contains  information  on  the  locations  of  the  first  and  last  cross  section 
of  the  profile.  The  user  has  to  enter  the  link  number  and  in-link  cross  section 
number,  for  example  cross  section  4 in  link  1 for  the  first  cross  section  and 
cross  section  5 in  link  3 for  the  last  cross  section. 

Line  4 contains  the  number  of  storm  events  (Ns)  for  which  output  is  requested 
for  that  particular  section.  Lines  5 through  Ns  + 4 contain  the  dates  of 
occurrence  of  the  storm  events. 

Lines  2 through  Ns  + 4 are  repeated  Np  times. 


Cross  Section  Data 


The  cross  section  data  file  contains  the  cross-sectional  geometry  and 
parameters  that  may  vary  from  cross  section  to  cross  section.  CONCEPTS 
divides  the  cross  section  into  a streambed,  left  and  right  banks,  and  left  and 
right  floodplains.  Data  blocks  related  to  these  geometric  elements  make  up  the 
input  file,  see  Table  4.11. 
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Table  4.11  Makeup  of  cross  section  data  input  file. 


data  block 

1 

2 

3 

4 

5 

6 


description 

general  data 
left  floodplain 
left  bank 
streambed 
right  bank 
right  floodplain 


Table  4.12  Data  block  #1  of  the  cross  section  data  input  file. 


line  number 

line  description 

data  type 

unit 

1 -4 

comment  lines 

5 

name  of  cross  section 

string 

6 

river  kilometer 

real 

km 

7 

comment  line 

8 

Manning  n 

real 

s/m 1/3 

9 

comment  line 

10 

tributary  inflow 

real 

General  Data  (Block  1) 

Data  block  1 consists  of  10  input  lines,  see  Table  4.12.  Line  5 contains  a 40- 
character  string,  which  is  the  name  of  the  input  file.  Line  6 contains  the  river 
kilometer  (km)  of  the  cross  section.  CONCEPTS  assumes  river  kilometer 
increases  in  streamwise  direction.  Line  8 contains  the  Manning  n of  the  total 
cross  section.  During  the  simulation,  the  equivalent  friction  factor  (see 
subsection  on  Flow  Resistance)  replaces  this  value.  Table  4. 1 3 lists  example  n 
values  for  various  types  of  channels  from  Chow  (1959).  Line  10  contains  the 
discharge  rate  of  a tributary  entering  the  reach  at  this  cross  section.  The 
discharge  rate  is  defined  as  a fraction  of  the  discharge  at  this  cross  section.  For 
example,  if  the  tributary  contributes  30  percent  of  the  discharge  at  the  cross 
section,  the  user  needs  to  enter  0.3. 

Streambed  Data  Block  (Block  4) 


The  block  containing  streambed  data  consists  of  1 5 + + 1 9 x N/a  lines  of 

input  (see  Table  4.14),  where  Nn  is  the  number  of  points  in  the  bed  profile  and 
N/a  is  the  number  of  soil  layers  in  the  bed.  Line  3 contains  the  number  of  points 
in  the  bed  profile  (see  Figure  1 .8).  Lines  5 through  4 + Nn  contain  the  station 
and  elevation  values  (m)  of  these  points. 
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Table  4.13  Values  of  the  roughness  coefficient  n (Chow,  1959). 

Type  of  Channel  and  Description 

Minimum  Normal 

Maximum 

Natural  streams  - minor  streams  (top  width  at  flood 
stage  < 100  ft) 

Streams  on  plain 

clean,  straight,  full  stage,  no  rifts  or  deep  pools 

0.025 

0.030 

0.033 

same  as  above,  but  more  stones  and  weeds 

0.030 

0.035 

0.040 

clean,  winding,  some  pools  and  shoals 

0.033 

0.040 

0.045 

same  as  above,  but  some  weeds  and  stones 

0.035 

0.045 

0.050 

same  as  above,  lower  stages,  more  ineffective 
slopes  and  sections 

0.040 

0.048 

0.055 

same  as  two  up,  more  stones 

0.045 

0.050 

0.060 

sluggish  reaches,  weedy,  deep  pools 

0.050 

0.070 

0.080 

very  weedy  reaches,  deep  pools,  or  floodways 
with  heavy  stand  of  timber  and  underbrush 

0.075 

0.100 

0.150 

Mountain  streams,  no  vegetation  in  channel,  banks 
usually  steep,  trees  and  brush  along  banks 
submerged  at  high  stages 

bottom:  gravels,  cobbles,  and  few  boulders 

0.030 

0.040 

0.050 

bottom:  cobbles  with  large  boulders 

0.040 

0.050 

0.070 

Natural  streams  - floodplains 
Pasture,  no  brush 

short  grass 

0.025 

0.030 

0.035 

high  grass 
Cultivated  areas 

0.030 

0.035 

0.050 

no  crop 

0.020 

0.030 

0.040 

mature  row  crops 

0.025 

0.035 

0.045 

mature  field  crops 
Brush 

0.030 

0.040 

0.050 

scattered  brush,  heavy  weeds 

0.035 

0.050 

0.070 

light  brush  and  trees,  in  winter 

0.035 

0.050 

0.060 

light  brush  and  trees,  in  summer 

0.040 

0.060 

0.080 

medium  to  dense  brush,  in  winter 

0.045 

0.070 

0.110 

medium  to  dense  brush,  in  summer 
Trees 

0.070 

0.100 

0.160 

dense  willows,  summer,  straight 

0.110 

0.150 

0.200 

cleared  land  with  tree  stumps,  no  sprouts 

0.030 

0.040 

0.050 

same  as  above,  but  with  heavy  growth  of  sprouts 

0.050 

0.060 

0.080 

heavy  stand  of  timber,  a few  down  trees,  little 
undergrowth,  flood  stage  below  branches 

0.080 

0.100 

0.120 

same  as  above,  but  with  flood  stage  reaching 
branches 

0.100 

0.120 

0.160 

Natural  streams  - major  streams  (top  width  at  flood 
stage  > 100  ft);  the  n value  is  less  than  that  for 
minor  streams  of  similar  description  because  banks 
offer  less  effective  resistance 

regular  section  with  no  boulder  or  brush 

0.025 

0.060 

irregular  and  rough  section 

0.035 

0.100 

Excavated  or  dredged 
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Table  4.13  Values  of  the  roughness  coefficient  n (Chow,  1 959).  (continued) 


Type  of  Channel  and  Description 

Minimum 

Normal 

Maximum 

Earth,  straight,  and  uniform 

clean,  recently  completed 

0.016 

0.018 

0.020 

clean,  after  weathering 

0.018 

0.022 

0.025 

gravel,  uniform  section,  clean 

0.022 

0.025 

0.030 

with  short  grass,  few  weeds 

0.022 

0.027 

0.033 

Earth  winding  and  sluggish 

no  vegetation 

0.023 

0.025 

0.030 

grass,  some  weeds 

0.025 

0.030 

0.033 

dense  weeds  or  aquatic  plants  in  deep  channels 

0.030 

0.035 

0.040 

earth  bottom  and  rubble  sides 

0.028 

0.030 

0.035 

stony  bottom  and  weedy  banks 

0.025 

0.035 

0.040 

cobble  bottom  and  clean  sides 

0.030 

0.040 

0.050 

Dragline-excavated  or  dredged 

no  vegetation 

0.025 

0.028 

0.033 

light  brush  on  banks 

0.035 

0.050 

0.060 

Rock  cuts 

smooth  and  uniform 

0.025 

0.035 

0.040 

jagged  and  irregular 

0.035 

0.040 

0.050 

Channels  not  maintained,  weeds  and  brush  uncut 

dense  weeds,  high  as  flow  depth 

0.050 

0.080 

0.120 

clean  bottom,  brush  on  sides 

0.040 

0.050 

0.080 

same  as  above,  highest  stage  of  flow 

0.045 

0.070 

0.110 

dense  brush,  high  stage 

0.080 

0.100 

0.140 

Lined  or  built-up  channels  - metal 

Smooth  steel  surface 

unpainted 

0.011 

0.012 

0.014 

painted 

0.012 

0.013 

0.017 

Lined  or  built-up  channels  - nonmetal 

Cement 

neat  surface 

0.010 

0.011 

0.013 

mortar 

0.011 

0.013 

0.015 

Wood 

planed,  untreated 

0.010 

0.012 

0.014 

planed,  creosoted 

0.011 

0.012 

0.015 

unplaned 

0.011 

0.013 

0.015 

plank  with  battens 

0.012 

0.015 

0.018 

lined  with  roofing  paper 

0.010 

0.014 

0.017 

Concrete 

trowel  finish 

0.011 

0.013 

0.015 

float  finish 

0.013 

0.015 

0.016 

finished,  with  gravel  on  bottom 

0.015 

0.017 

0.020 

unfinished 

0.014 

0.017 

0.020 

gunite,  good  section 

0.016 

0.019 

0.023 

gunite,  wavy  section 

0.018 

0.022 

0.025 

on  good  excavated  rock 

0.017 

0.020 

on  irregular  excavated  rock 

0.022 

0.027 
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Table  4.13  Values  of  the  roughness  coefficient  n (Chow,  1959).  (continued) 


Type  of  Channel  and  Description 

Minimum 

Normal 

Maximum 

Concrete  bottom  float  finish  with  sides  of: 

dressed  stone  in  mortar 

0.015 

0.017 

0.020 

random  stone  in  mortar 

0.017 

0.020 

0.024 

cement  rubble  masonry,  plastered 

0.016 

0.020 

0.024 

cement  rubble  masonry 

0.020 

0.025 

0.030 

dry  rubble  or  riprap 

0.020 

0.030 

0.035 

Gravel  bottom  with  sides  of: 

formed  concrete 

0.017 

0.020 

0.025 

random  stone  mortar 

0.020 

0.023 

0.026 

dry  rubble  or  riprap 

0.023 

0.033 

0.036 

Brick 

glazed 

0.011 

0.013 

0.015 

in  cement  mortar 

0.012 

0.015 

0.018 

Masonry 

cemented  rubble 

0.017 

0.025 

0.030 

dry  rubble 

0.023 

0.032 

0.035 

Dressed  ashlar 

0.013 

0.015 

0.017 

Asphalt 

smooth 

0.013 

0.013 

rough 

0.016 

0.016 

Vegetal  lining 

0.030 

0.500 

Table  4.14  Input  block  #4  with  streambed  data. 

r~ 

3 

CD 

CO 

-fc. 

+ 

through 

32  + N„  are 

repeated  N,a  times. 

line  number  line  description 

data  type 

unit 

1-2  comment  lines 


3 

number  of  points  in  profile  Nt 

4 

comment  line 

5 - 5 + Nn  - 1 

station  and  elevation 

5 + /V„ 

comment  line 

6 + N„ 

elevation  of  bed  rock 

7 + Nn 

comment  line 

8 + 

porosity  of  streambed 

9 + /V„ 

comment  line 

1 0 + 7V„ 

hiding  factors 

11  +NU  & 12  + Nn 

comment  lines 

13  + N„ 

number  of  soil  layers  Nla 

14  + Af„ 

comment  line 

15  + N„ 

depth  below  bed  surface 

\6  + Nn 

comment  line 

1 7 + Nn  - 29  + Nn 

bed  composition 

30  + N„  & 31  + Na 

comment  lines 

32  + Nn 

cohesive  bed  parameters 

14  + N„  + 1 9 x N/a 

comment  line 

15  + /V„  + 19  xNla 

Manning  n 

integer 

real 

real 

real 

real 

integer 

real 

real 

real 

real 


m 

m 


m 

% 

Pa,  Pa,  m/S'Pa 
s/m 1/3 
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Table  4.15  Input  block  #3  and  #5  with  streambank  data. 

line  number 

line  description 

data  type 

unit 

1 -2 

comment  lines 

3 

number  of  points  in  profile  Nn 

integer 

4 

comment  line 

5 - 5 + Nn  - 1 

station  and  elevation 

real 

m 

5 + Nn 

comment  line 

6 + N„ 

bank  material  properties 

real 

Pa,  °,  N/m3 

1 + Nn 

comment  line 

8 + Nn  - 20  + Nn 

bank  material  composition 

real 

% 

21  + N„ 

comment  line 

22  + Nn 

critical  shear  stress 

real 

Pa 

23  + A rn 

comment  line 

24  + N„ 

bank  moisture  parameters 

real 

m,  0 

25  + Nn 

comment  line 

26  + Nn 

Manning  n 

real 

s/m 1/3 

Lines  6 + Nn , 8 + Nn , and  10  + Nn  contain  bedrock  elevation  (m),  streambed 

porosity,  and  hiding  factors  for  silt,  sand,  and  gravel  in  the  surface  layer, 
respectively.  The  channel  cannot  incise  below  bedrock  elevation. 

Line  13  + N contains  the  number  of  soil  layers  comprising  the  bed.  The  next 

lines  contain  information,  for  each  layer,  on  depth  below  the  bed  surface  (m), 
composition  (%),  and  erodibility  of  cohesive  material.  This  block  is  repeated 
N/a  times.  The  user  has  to  enter  the  fractional  content  of  each  of  the  1 3 size 
classes  as  percentages.  The  parameters  of  a cohesive  soil  layer  are  entered  as: 

( 1 ) critical  shear  stress  to  deposit  sediment  particles  (xd)  in  Pa,  (2)  critical  shear 
stress  to  entrain  sediment  particles  (xe)  in  Pa,  and  (3)  erodibility  coefficient  in 
m/s-Pa. 

Line  15  + Nn  + 19  x N{a  contains  the  Manning  n of  the  streambed. 

Streambank  Data  Block  (Blocks  3 and  5) 


Blocks  3 and  5 containing  streambank  data  consist  of  26  + Nn  lines  of  input 
(see  Table  4. 1 5),  where  Nn  is  the  number  of  points  in  the  bank  profile.  Line  3 
contains  the  number  of  points  in  the  bank  profile  (see  Figure  1 .8).  Lines  5 
through  4 + Nn  contain  the  station  and  elevation  values  (m)  of  these  points. 

Line  6 + Nn  contains  effective  cohesion  (Pa),  effective  angle  of  internal 
friction  (°),  the  angle  indicating  the  increase  in  shear  strength  for  an  increase  in 
matric  suction  (°),  and  the  bulk  weight  of  the  bank  material  (N/m3).  Lines 
8+tV„  through  20  + Nn  contain  the  fractional  content  of  each  of  the  1 3 size 
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Table  4.16  Input  block  #2  and  #6  with  floodplain  data. 


line  number 

line  description 

data  type 

1 -2 

comment  lines 

3 

number  of  points  in  profile  Nn 

integer 

4 

comment  line 

5 - 5 + Nn  - 1 

station  and  elevation 

real 

5 + Nn 

comment  line 

6 + Nn 

Manning  n 

real 

classes  as  percentages.  Line  22  + Nn  contains  the  critical  shear  stress  to  entrain 
bank  material  particles  (Pa).  Line  24  + Nn  contains  groundwater  table  (m), 
hydraulic  gradient  (m/m),  and  angle  of  seepage  force  (°).  Line  26  + Nn 
contains  Manning  n of  the  streambank. 

Floodplain  Data  Block  (Blocks  2 and  6) 


Blocks  2 and  6 containing  floodplain  data  consist  of  6 + N n lines  of  input  (see 
Table  4.16),  where  Nn  is  the  number  of  points  in  the  floodplain  profile.  Line  3 
contains  the  number  of  points  in  the  floodplain  profile  (see  Figure  1 .8).  Lines  5 
through  4 + Nn  contain  the  station  and  elevation  values  (m)  of  these  points. 

Line  6 + Nn  contains  the  Manning  n of  the  floodplain. 

Hydraulic  Structure  Data 

The  four  types  of  hydraulic  structures  (culvert,  bridge  crossing,  drop  structure, 
and  generic  structure)  have  common  parameters  and  structure  specific 
parameters.  The  common  parameters  comprise  the  first  part  of  each  structure’s 
input  file,  followed  by  the  structure  specific  data. 

Table  4.17  lists  the  common  data:  a name  field  on  line  5 of  the  input  file,  river 
kilometer  (km)  on  line  7,  Manning  n on  line  9,  length  of  the  structure  (m)  on 
line  1 1,  the  upstream  and  downstream  inverts  (m)  on  line  13,  and  the  upstream 
and  downstream  elevations  of  the  structure  above  the  streambed  (m)  on  line 
15. 

Culvert 


The  flow  computation  at  culverts  is  based  on  the  U.S.  Federal  Highway 
Administration’s  (1985)  nomographs.  CONCEPTS  can  simulate  the  flow  at 
box  and  pipe  culverts.  Presently,  we  have  not  implemented  other  culvert 
shapes.  All  input  parameters  for  pipe  and  box  culverts  are  the  same,  except  for 
their  geometry. 
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Table  4.17  Input  parameters  common  to  the  different  types  of  structures. 


line  number 

line  description 

data  type 

unit 

1 -4 

comment  lines 

5 

name  of  structure 

string 

6 

comment  line 

7 

river  kilometer 

real 

km 

8 

comment  line 

9 

Manning  n 

real 

s/m 1/3 

10 

comment  line 

11 

length 

real 

m 

12 

comment  line 

13 

inverts 

real 

m 

14 

comment  line 

15 

elevations  above  streambed 

real 

m 

Table  4.18 

Input  data  for  culverts. 

line  number  line  description 

data  type 

unit 

16-  17 

comment  lines 

18 

chart  and  scale  numbers 

integer 

19 

comment  line 

20 

entrance  loss  coefficient 

real 

21 

comment  line 

22 

number  of  barrels 

integer 

23-24 

comment  lines 

25 

dimensions 

real 

m 

Table  4.18  lists  the  input  data  required  for  culverts.  Line  18  contains  the 
USFHWA  (1985)  chart  and  scale  numbers  representing  the  shape  of  the 
culvert  (see  Table  4.19).  Line  20  contains  the  entrance  loss  coefficient  used 
when  the  culvert  flow  is  controlled  by  the  outlet  (see  Table  4.20).  Line  22 
contains  the  number  of  culvert  barrels  in  the  road  crossing.  Line  25  contains 
the  dimensions  of  the  culvert  barrel:  diameter  (m)  for  a pipe  culvert,  and  span 
(m)  and  rise  (m)  for  a box  culvert. 

Bridge  Crossing 

CONCEPTS  assumes  that  the  shape  of  the  bridge  crossing  is  trapezoidal  with  a 
horizontal  bed.  Table  4.21  lists  the  input  data  required  for  bridge  crossings. 
Line  1 8 contains  the  bottom  width  (m)  and  side  slope  of  the  cross  section.  Line 
20  contains  the  total  pier  width  (m),  pier  shape  coefficient  (Table  4.22),  and 
pier  loss  coefficient. 
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Table  4.19  USFHWA  (1985)  chart  and  scale  numbers  for  pipe  and  box  culverts. 


Chart 

Number 

Scale 

Number 

Description 

1 

Concrete  Pipe  Culvert 

1 

square-edged  entrance  with  headwall 

2 

groove  end  entrance  with  headwall 

3 

groove  end  entrance,  pipe  projecting  from  fill 

2 

Corrugated  Metal  Pipe  Culvert 

1 

headwall 

2 

mitered  to  conform  slope 

3 

pipe  projecting  from  fill 

3 

Concrete  Pipe  Culvert,  beveled  ring  entrance 

1 

small  bevel 

2 

large  bevel 

8 

Box  Culvert  with  Flared  Wingwalls 

1 

wingwalls  flared  30  to  75  degrees 

2 

wingwalls  flared  90  to  15  degrees 

3 

wingwalls  flared  0 degrees  (sides  extended  straight) 

9 

Box  Culvert  with  Flared  Wingwalls  and  Inlet  Top  Edge  Bevel 

1 

wingwall  flared  45  degrees,  inlet  top  edge  bevel  = 0.43D 

2 

wingwall  flared  18  to  33.7  degrees,  inlet  top  edge  bevel  = 0.083 D 

10 

Box  Culvert,  90-degree  Headwall,  Chamfered  or  Beveled  Inlet 
Edges 

1 

inlet  edges  chamfered  3/4-inch 

2 

inlet  edges  beveled  1 /2-in/ft  at  45  degrees 

3 

inlet  edges  beveled  1 -in/ft  at  33.7  degrees 

11 

Box  Culvert,  Skewed  Headwall,  Chamfered  or  Beveled  Inlet  Edges 

1 

headwall  skewed  45  degrees,  inlet  edges  chamfered  3/4-in 

2 

headwall  skewed  30  degrees,  inlet  edges  chamfered  3/4-in 

3 

headwall  skewed  15  degrees,  inlet  edges  chamfered  3/4-in 

4 

headwall  skewed  15  to  45  degrees,  inlet  edges  beveled 

12 

Box  Culvert,  Non-Offset  Flared  Wingwalls,  3/4-inch  Chamfer  at  Top 
of  Inlet 

1 

wingwalls  flared  45  degrees,  inlet  not  skewed 

2 

wingwalls  flared  18.4  degrees,  inlet  not  skewed 

3 

wingwalls  flared  18.4  degrees,  inlet  skewed  30  degrees 

13 

Box  Culvert,  Offset  Flared  Wingwalls,  Beveled  Edge  at  Top  of  Inlet 

1 

wingwalls  flared  45  degrees,  inlet  top  edge  bevel  = 0.042D 

2 

wingwalls  flared  33.7  degrees,  inlet  top  edge  bevel  = 0.083 D 

3 

wingwalls  flared  18.4  degrees,  inlet  top  edge  bevel  = 0.083 D 

Drop  Structure 

CONCEPTS  assumes  that  the  cross  section  of  the  drop  structure  is  trapezoidal 
with  a horizontal  bottom.  Table  4.23  lists  the  input  data  required  for  drop 
structures.  Line  1 8 contains  the  bottom  width  (m)  and  side  slope  of  the  drop 
structure.  Line  20  contains  the  entrance  loss  coefficient. 
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Table  4.20  Entrance  loss  coefficients  for  pipe  and  box  culverts. 


Type  of  Structure  and  Design  of  Entrance 

Concrete  Pipe  Projecting  from  Fill  (no  headwall) 
socket  end  of  pipe  (grooved  end) 
square  cut  end  of  pipe 

Concrete  Pipe  with  Headwall  or  Headwall  and  Wingwalls 
socket  end  of  pipe 
square  cut  end  of  pipe 
rounded  entrance 

Concrete  Pipe 

mitered  to  conform  to  fill  slope 
end  section  conformed  to  fill  slope 
beveled  edges,  33.7  and  45  degree  bevels 
side  slope  tapered  inlet 

Corrugated  Metal  Pipe  or  Pipe-Arch 
projected  from  fill  (no  headwall) 
headwall  or  headwall  and  wingwalls  square  edge 
mitered  to  conform  to  fill  slope 
end  section  conformed  to  fill  slope 
beveled  edges,  33.7  and  45  degree  bevels 
side  slope  tapered  inlet 

Concrete  Box,  Headwall  Parallel  to  Embankment  (no  wingwalls) 
square-edged  on  three  sides 

three  edges  rounded  to  radius  of  1/12  barrel  dimension 

Concrete  Box,  Wingwalls  at  30  to  75  degrees  to  Barrel 
square-edge  on  crown 

top  corner  rounded  to  radius  of  1/12  barrel  dimension 

Concrete  Box,  Wingwalls  at  10  to  25  degrees  to  barrel 
square-edge  on  crown 

Concrete  Box,  Wingwalls  parallel  (extension  of  sides) 
square-edge  on  crown 
side  or  slope  tapered  inlet 


Loss  Coefficient 

0.2 

0.5 

0.2 

0.5 

0.2 

0.7 

0.5 

0.2 

0.2 

0.9 

0.5 

0.7 

0.5 

0.2 

0.2 

0.5 

0.2 

0.4 

0.2 

0.5 

0.7 

0.2 


Table  4.21  Input  data  for  bridge  crossings. 


line  number 

line  description 

data  type 

unit 

16-17 

comment  lines 

18 

bridge  crossing  geometry 

real 

m,  - 

19 

comment  line 

20 

pier  parameters 

real 

m,  -,  - 
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Table  4.22  Yarnell’s  pier  shape  coefficient  for  various  pier  shapes. 

Pier  Shape  Yarnell  Cp  Coefficient 


semi-circular  nose  and  tail  0.90 

lens-shaped  nose  and  tail  0.90 

twin-cylinder  piers  with  connecting  diaphragm  0.95 

twin-cylinder  piers  without  diaphragm  1 .05 

90-degree  triangular  nose  and  tail  1 .05 

square  nose  and  tail  1 .25 


Table  4.23  Input  data  for  drop  structures. 

line  number  description  data  type  unit 

16-17  comment  lines 

18  drop  structure  geometry  real  m,  m/m 

19  comment  line 

20  entrance  loss  coefficient  real 


Table  4.24  Input  data  for  generic  structures. 


line  number 

16-  17 
18 
19 

20  - 20  + Nse  - 1 

20  + A 'se 

2\+Nse-22  + Nst 

23  + ^e 

24  + Nse 

25  + Nse- 
25  + N+N.-  1 


line  description  data  type 

comment  lines 

number  of  wall  segments  Nse  integer 

comment  line 

segment  elevation  and  slope  real 

top  elevation  real 

comment  lines 

number  of  rating  curve  segments  Nrs  integer 

comment  line 

rating  curve  parameters  real 


unit 


m,  m/m 
m 


m,  -, 


Generic  Structure 


A generic  structure  is  any  structure  for  which  a rating  curve  is  available.  The 
main  input  data  are  the  shape  of  the  structure  and  the  rating  curve.  CONCEPTS 
assumes  that  the  cross  section  of  the  walls  of  the  structure  consist  of  linear 
elements  characterized  by  an  elevation  and  slope  (see  Figure  4.1a).  The  rating 
curve  may  comprise  up  to  4 segments  (Figure  4.  lb).  Each  segment  is  a power 
function. 

Table  4.24  lists  the  input  data  required  for  generic  structures.  Line  18  contains 
the  number  of  segments,  Nse, comprising  the  walls  of  the  structure.  Lines  20 

through  19  + Nse  contain  the  bottom  elevation  (m)  and  slope  of  each  segment 
(Figure  4.1a).  Line  20  + N contains  the  elevation  of  the  top  of  the  structure 
(m). 
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Figure  4.1  Shape  of  a generic  structure  (a)  and  rating  curve  (b).  S,  is  side  slope  and  z, 
is  starting  elevation  of  wall  segment  /.  /i,  is  breakpoint  depth  in  the  rating 
curve. 


zi 


(a) 


Line  23  + Nse  contains  the  number  of  segments  comprising  the  rating  curve, 

Nrs.  Lines  25  + Nse  through  24  + Nse  + Nrs  contain  the  starting  depth  (m)  (see 

Figure  4.  lb)  and  coefficient  and  exponent  of  the  power  function  for  each 
segment. 
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Table  4.25  Layout  of  each  record  in  the  upstream  boundary  conditions  file. 


width 
data  type 
unit 

example 


date  field 

runoff  status 

flow  discharge 

sediment  load 

column  1 

column  2 

column  3 

columns  4-16 

19 

3 

10 

10 

string 

integer 

real 

real 

m* 1 2 3 4/s 

kg/s 

06/07/1999  14:00:00 

1 

10.45 

1.45E+01 

Dynamic  Upstream  Boundary  Conditions 

CONCEPTS  simulates  unsteady  flow.  The  user  has  to  specify  breakpoint  flow 
and  sediment  discharge  data  at  the  upstream  boundary  of  the  modeling  reach. 
Each  record  of  the  discharge  input  file  consists  of  up  to  16  fields  (Table  4.25): 

1 a 19-character  date  field, 

2 a 3-character  runoff  status  field, 

3 a 1 0-character  flow  discharge  field,  and 

4 13  10-character  sediment  load  fields. 

The  first  three  columns  are  mandatory.  If  the  user  sets  the  flag  indicating  the 
upstream  boundary  condition  option  regarding  sediment  load  to  one  (line  14  of 
the  run  control  file),  columns  4 through  16  contain  the  sediment  load  entering 
the  stream  corridor  for  each  size  fraction  (kg/s). 

The  date  is  formatted  as  “mm/dd/yyyy  hh:mm:ss,”  for  example  “06/07/1999 
14:00: 00”.  The  storm  status  informs  CONCEPTS  on  the  beginning  and  end  of 
runoff  events.  It  takes  the  following  values: 

• a value  of  1 , beginning  of  runoff  event; 

• a value  of  2,  end  of  runoff  event;  and 

• a value  of  0,  indicates  data  between  beginning  and  end  of  a runoff  event  or 
data  between  runoff  events. 

The  flow  discharge  is  represented  by  a real  number  with  two  significant  digits 
after  the  decimal  point.  The  sediment  loads  are  represented  by  real  numbers  in 
exponential  form  with  two  significant  digits  after  the  decimal  point. 
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This  chapter  presents  the  output  data  structure  of  CONCEPTS  and 
a way  to  process  these  data.  Depending  on  the  output  options 
selected  by  the  user,  CONCEPTS  creates  three  types  of  output  files: 
(1)  output  at  a certain  location  and  for  a certain  runoff  event,  (2) 
time-series  output  at  a certain  location,  and  (3)  output  for  a certain 
runoff  event  along  a section  of  the  modeling  reach. 

Output  Data  File  Specification 

CONCEPTS  can  generate  three  types  of  output  files,  see  section  on  “Run 
Control  Data”  in  Chapter  4.  There  are  three  output  categories: 

1 output  at  a certain  location  and  for  a certain  runoff  event, 

2 time-series  output  at  a certain  location,  and 

3 output  for  a certain  runoff  event  along  a section  of  the  modeling  reach. 

The  names  of  the  output  files  follow  the  following  convention:  run  identifier  + 

‘ ’ + number  + ‘ .txt’.  For  example,  if  the  run  identifier  is  simulation,  then  the 
output  file  names  are  ‘simulation  oso.txt’,  ‘simulation  051.txt’,  etc.  (see 
also  Chapter  4,  page  69).  The  layout  of  the  output  file  differs  for  each  output 
category. 

Output  at  a Certain  Location  and  for  a Certain  Runoff 
Event 


If  the  user  requests  output  at  a certain  location  and  for  a certain  runoff  event, 
CONCEPTS  generates  an  output  file  for  each  requested  cross  section.  Table 
4.6  lists  the  parameters  that  the  user  can  request  as  output.  The  output  file 
comprises  initial  data  (header)  and  data  for  each  runoff  event. 
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example  header 


example  initial 
geometry  of  cross 
section 


example  initial  in- 
bank channel 
widths 


Initial  Data 

The  initial  data  consist  of  a header  and  data  that  depend  on  the  requested 
output  parameters.  The  header  comprises  the  name  of  the  cross  section,  its 
river  kilometer,  type  of  data,  and  dates  of  the  runoff  events  for  which  output  is 
requested.  For  example: 

Output  file  for  cross  section  "First  cross  section"  at  river  kilometer 
0.676 


Output  type  = 17536 

Output  dates  are:  12/27/1982  16:00:00 
05/18/1983  16:45:00 
09/20/1983  08:30:00 
11/19/1983  19:30:00 


CONCEPTS  also  produces  initial  data  for  the  following  parameters  (cf.  Table 
4.6): 

• cross-sectional  geometry  (data  type  1024), 

INITIAL  CROSS  SECTION  GEOMETRY 

STATION  ELEVATION 


(M) 

(M) 

0. 

000 

76. 

950 

1. 

000 

75. 

950 

101. 

400 

75. 

950 

102. 

590 

74  . 

4 60 

107. 

470 

73. 

910 

119. 

230 

73. 

150 

123. 

800 

73. 

150 

130. 

910 

73. 

050 

135. 

810 

72. 

970 

138. 

220 

72. 

570 

144. 

160 

72. 

170 

144. 

960 

71. 

900 

146. 

360 

71. 

840 

151. 

330 

71. 

650 

152. 

300 

71. 

810 

153. 

430 

72. 

600 

156. 

290 

75. 

830 

158. 

700 

75. 

860 

258. 

700 

75. 

860 

259. 

700 

76. 

860 

• in-bank  top  and  bottom  widths  of  cross  section  (data  type  2048), 

INITIAL  CHANNEL  WIDTHS 

BOTTOM  FLOODPLAIN 
(M)  (M) 

33.070  54.890 
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example  initial 
bank  heights 


example  initial 
characteristic  sed- 
iment sizes 


example  initial 
sediment  compo- 
sition 


example  storm 
date  output  line 


Output  Data  File  Specification 

• height  of  left  and  right  banks  of  cross  section  (data  type  4096), 

INITIAL  BANK  HEIGHTS 

LEFT  RIGHT 

(M)  (M) 

2.800  4.020 

• characteristic  sediment  particle  sizes  of  surface  layer  of  streambed  (data  type 
16384),  and 

INITIAL  BED  MATERIAL  CHARACTERISTICS 


D16  D50  D84  D90  DMEAN 

(MM)  (MM)  (MM)  (MM)  (MM) 

0.227  3.317  16.740  20.889  3.245 


• composition  of  surface  layer  of  streambed  (data  type  32768). 

INITIAL  SIZE  FRACTION  DISTRIBUTION 

DIAMETER  FRACTION 


(MM) 

(-; 

i _ 

0. 

004 

0. 

.000 

0. 

016 

0. 

.000 

0. 

040 

0. 

.000 

0. 

127 

1. 

.479 

0. 

458 

32  . 

.143 

1. 

297 

8. 

.354 

2. 

594 

4 . 

.877 

4 . 

362 

6. 

,653 

7 . 

336 

10. 

.075 

12. 

338 

12. 

.073 

20. 

749 

14  . 

.220 

32. 

000 

8. 

.112 

43. 

713 

2. 

.014 

Runoff  Event  Related  Data 

CONCEPTS  writes  the  output  data  to  a file  following  the  sequence  shown  in 
Table  4.6,  but  grouping: 

• peak  discharge,  peak  flow  depth,  peak  stage,  and  peak  friction  slope; 

• sediment  yield  and  cumulative  sediment  yield;  and 

• change  in  bed  elevation,  cumulative  change  in  bed  elevation,  lateral  erosion, 
and  cumulative  lateral  erosion. 

The  first  line  in  the  output  file  for  each  runoff  event  is  the  date  of  the  storm 
event: 

OUTPUT  FOR  THE  STORM  EVENT  STARTED  ON  12/25/1982  11:58:  0 

This  line  is  followed  by  eight  groups  of  output: 
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1  peak  hydraulic  variables  (data  types  l,  2,  4,  and  8), 


example  output  of 
peak  discharge 
and  peak  flow 
depth 


TIME  TO  PEAK  DISCHARGE  TIME  TO  PEAK  DEPTH 

(CMS)  (M) 

12/25/1982  11:58:  0 10.460  12/25/1982  11:58:  0 2.451 


2  sediment  yield  (data  types  16  and  32), 


example  output  of  sediment  yield 


STORM  EVENT  GENERATED  SEDIMENT  YIELD  CUMULATIVE  SEDIMENT  YIELD 

SILT  YLD  SAND  YLD  GRAVEL  YLD  TOTAL  YLD  SILT  YLD  SAND  YLD  GRAVEL  YLD  TOTAL  YLD 
(TONS)  (TONS)  (TONS)  (TONS)  (TONS)  (TONS)  (TONS)  (TONS) 
4.93  71.24  6.52  82.69  1400.69  1419.18  144.25  2964.12 


3  cross-sectional  changes  (data  types  64,  128,  256,  and  512), 

example  output  of  cross-sectional  changes 

cross-sectional 

changes  bed  chng  cum  bed  lat  chnc-  cum  lat 

(M)  (M)  (M)  (M) 

-0.016  -0.208  0.384  5.495 


4  cross-section  geometry  (data  type  1024), 

example  output  of  cross  section  geometry 

cross  section 

geometry  station  elevation 


(M) 

(M) 

0. 

000 

76. 

950 

1. 

000 

75. 

950 

101. 

400 

75. 

950 

102. 

590 

74  . 

460 

107. 

470 

73. 

908 

119. 

230 

73. 

109 

123. 

800 

73. 

109 

130. 

910 

72. 

992 

135. 

810 

72. 

897 

138. 

220 

72. 

422 

144  . 

160 

71. 

961 

144  . 

960 

71. 

690 

146. 

360 

71. 

632 

151. 

330 

71. 

442 

152. 

300 

71. 

602 

157. 

795 

71. 

602 

157. 

795 

73. 

111 

160. 

520 

75. 

860 

258. 

700 

75. 

860 

259. 

700 

76. 

860 

5  in-bank  top  and  bottom  width  of  cross  section  (data  type  2048), 

example  output  of  channel  widths 
in-bank  channel 

widths  bottom  floodplain 

(M)  (M) 

38.565  59.120 
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example  output  of 
bank  heights 


example  output  of 
characteristic 
sediment  sizes 


example  output  of 
composition  of 
surface  layer 


example  header 


6 height  of  left  and  right  banks  of  cross  section  (data  type  4096), 

BANK  HEIGHTS 

LEFT  RIGHT 

(M)  (M) 

2.841  4.258 


7 characteristic  sediment  particle  sizes  of  surface  layer  of  streambed  (data  type 
16384),  and 


BED  MATERIAL 

CHARACTERISTICS 

D16 

D50 

D84 

D90 

DMEAN 

(MM) 

(MM) 

(MM) 

(MM) 

(MM) 

0.236 

4.566 

17.784 

22.164 

3.607 

8 composition  of  surface  layer  of  streambed  (data  type  32768). 


SIZE  FRACTION  DISTRIBUTION 
DIAMETER  FRACTION 


(MM) 

(-) 

0. 

004 

0. 

.000 

0. 

016 

0. 

.000 

0. 

040 

0. 

.000 

0. 

127 

2. 

.126 

0. 

458 

18. 

.194 

1. 

297 

2. 

.553 

2. 

594 

0. 

. 647 

4 . 

362 

9. 

.574 

7. 

336 

14  . 

.498 

12. 

338 

17. 

.373 

20. 

749 

20. 

.463 

32. 

000 

11. 

.673 

43. 

713 

2. 

.898 

Time-Series  Output  at  a Certain  Location 

If  the  user  requests  time-series  output,  CONCEPTS  generates  an  output  file  for 
each  requested  cross  section.  Table  4.8  lists  the  parameters  the  user  can  request 
as  output.  The  output  file  comprises  initial  data  (header)  followed  by  n + 1 
columns  of  data,  where  n is  the  number  of  outputted  variables. 

Initial  Data 

The  initial  data  comprise  the  name  of  the  cross  section,  river  kilometer,  type  of 
output  data,  and  the  start  and  end  dates  of  each  time  series. 

Output  time-series  file  for  cross  section  "First  cross  section"  at 
river  kilometer  0.248 

Output  type  = 595973 

Output  periods  are:  from  10/01/1981  00:00:00  to  12/31/1591  24:00:00 
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CONCEPTS  then  prints  the  header  for  the  multi-column  output,  the  first 
column  of  which  is  time  followed  by  the  name  of  the  outputted  variable  and  its 
unit.  Following  is  a list  of  all  possible  columns  in  the  header: 

• time, 

TIME 


• discharge  (data  type  l), 

DISCHARGE 

(CMS) 


• flow  velocity  (data  type  2), 

VELOCITY 

(M/S) 


• flow  depth  (data  type  4), 

DEPTH 

(M) 


• stage  (data  type  8), 

STAGE 

(M) 


• flow  area  (data  type  16), 

AREA 

(M2) 


• flow  top  width  (data  type  32), 

TOP  WIDTH 
(M) 


• wetted  perimeter  (data  type  64), 

PERIMETER 

(M) 


• hydraulic  radius  (data  type  128), 

RADIUS 

(M) 


• conveyance  (data  type  256), 

CONVEYANCE 

(CMS) 


• friction  slope  (data  type  512), 

F. SLOPE 


• energy  head  (data  type  1024), 

HEAD 

(M) 


• Froude  number  (data  type  2048), 

FROUDE 
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• bed  shear  stress  (data  type  4096), 

BED  SHEAR 
(Pa) 


• sediment  discharge  (silt/sand/gravel/total)  (data  type  8192), 

SILT  DIS  SAND  DIS  GRAVEL  DIS  TOTAL  DIS 
(CMS)  (CMS)  (CMS)  (CMS) 


• cumulative  sediment  yield  (data  type  16384), 

SILT  YLD  SAND  YLD  GRAVEL  YLD  TOTAL  YLD 
(TONS)  (TONS)  (TONS)  (TONS) 


• cumulative  change  in  bed  elevation  (data  type  32768), 

CUM  BED 
(M) 


• thalweg  elevation  (data  type  65536), 

THALWEG 

(M) 


• cumulative  lateral  erosion  (data  type  131072), 

CUM  LAT 
(M) 


• factor  of  safety  of  left  and  right  banks  (data  type  524288), 

SAFETY  L SAFETY  R 


• apparent  cohesion  of  left  and  right  banks  (data  type  1048576), 

AP  COH  L AP  COH  R 
(Pa)  (Pa) 


• pore-water  force  of  left  and  right  banks  per  unit  width  (data  type  2097152), 

POPE  L POPE  P. 

(N/M)  (N/M) 

• matric  suction  force  of  left  and  right  banks  per  unit  width  (data  type 
4194304), 

MATPIO  L M ATP  10  P 
(N/M)  (N/M) 

• weight  of  failure  block  of  left  and  right  banks  per  unit  width  (data  type 
8388608), 

W ELK  L W ELK  R 
(N/M)  (N/M) 

• weight  of  water  on  the  bank  of  left  and  right  banks  per  unit  width  (data  type 
16777216), 

W WATER  L W WATEP  R 
(N/M)  (N/M) 

• horizontal  component  of  the  confining  force  of  left  and  right  banks  per  unit 
width  (data  type  33554432), 

HYD  PR  L HYD  PR  R 
(N/M)  (N/M) 
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• groundwater  elevation  of  left  and  right  banks  (data  type  67108864),  and 

PHPEA  L PHPEA  P 

(M) (M) 

• elevation  of  bank  top  of  left  and  right  banks  (data  type  134217728). 

BANKTOP  L EANKTOP  F 
(M)  (M) 


Continuous  Data 

After  every  completed  time  step  CONCEPTS  checks  if  time  falls  between  the 
start  and  end  time  of  all  requested  time  series.  When  the  model  time  is  within 
the  boundaries  of  a time  series,  CONCEPTS  prints  the  requested  parameters  to 
the  corresponding  output  file. 

Output  for  a Certain  Runoff  Event  along  a Section  of  the 
Modeling  Reach 

If  the  user  requests  output  for  a certain  runoff  event  along  a section  of  the 
modeling  reach,  referred  to  as  a profile , CONCEPTS  generates  an  output  file 
for  each  requested  profile.  Table  4. 1 0 lists  the  parameters  that  the  user  can 
request  as  output.  The  output  file  comprises  initial  data  (header)  and  data  for 
each  runoff  event.  CONCEPTS  writes  the  output  data  following  a multi- 
column  layout.  The  first  three  lines  are  header  information  and  the  first  column 
is  river  kilometer. 


Initial  Data 

The  initial  data  consist  of  a header  and  data  that  depends  on  the  requested 
output  parameters.  The  header  comprises  the  name  of  the  cross  section,  its 
river  kilometer,  type  of  data,  and  dates  of  the  runoff  events  for  which  output  is 
requested.  For  example: 

example  header  Output  profile  file  for  cross  section  "First  cross  section"  at  river 
kilometer  0.676 

to  cross  section  "Last  cross  section"  at  river  kilometer  2.906 


Output  type  = 135 

Output  dates  are:  12/27/1982  16:00:00 
05/18/1983  16:45:00 
09/20/1983  08:30:00 
11/19/1983  19:30:00 


CONCEPTS  also  produces  initial  data  for  the  following  parameters  if 
requested: 
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• thalweg  elevation  (data  type  4), 


example  initial 

thalweg  elevation 

RIVER  KM 

THALWEG 

profile 

(KM) 

(M) 

0.248 

73.270 

0.499 

72.110 

0.819 

71.080 

1.002 

70.620 

1.348 

69.640 

1.593 

69.160 

1.831 

68.420 

2.148 

68.120 

2.421 

67.330 

2.606 

67.100 

2.726 

66.950 

2.906 

66.800 

• in-bank  top  width  of  cross  section  (data  type  16), 

example  initial  in- 

bank  channel  top 

RIVER  KM 

TOP  WIDTH 

width  profile 

(KM) 

(M) 

0.248 

88.660 

0.499 

78.970 

0.819 

27.740 

1.002 

29.080 

1.348 

30.780 

1.593 

76.320 

1.831 

41.790 

2.148 

40.110 

2.421 

27.500 

2.606 

27.700 

2.726 

51.940 

2.906 

34.840 

• height  of  left  and  right  bank  of  cross  section  (data  type  32),  and 

example  initial 

BANK  HEIGHT 

bank-height  profile 

RIVER  KM 

LEFT  RIGHT 

(KM) 

(M) 

(M) 

0.248 

2.450 

3.420 

0.499 

2.350 

2.500 

0.819 

3.840 

4 . 510 

1.002 

1.560 

3.990 

1.348 

4.720 

2.860 

1.593 

4.510 

3.320 

1.831 

5.190 

1.550 

2.148 

3.660 

2.560 

2.421 

4.750 

3.900 

2.606 

5.230 

3.260 

2.726 

3.870 

1.150 

2.906 

3.110 

3.170 
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characteristic  sediment  particle  sizes  of  surface  layer  of  streambed  (data  type 
128). 


example  profile  of 
initial  characteris- 
tic sediment  sizes 


BED 

MATERIAL  < 

CHARACTERISTIC  DIAMETERS 

RIVER 

KM 

D16 

D50 

D84 

D90 

DMEAN 

(KM) 

(MM) 

(MM) 

(MM) 

(MM) 

(MM) 

0. 

248 

0. 

227 

3. 

317 

16. 

740 

20. 

889 

3. 

245 

0. 

499 

0. 

227 

3. 

317 

16. 

740 

20. 

889 

3. 

245 

0. 

819 

0. 

227 

O 

. 

317 

16. 

740 

20. 

889 

3. 

245 

1. 

002 

0. 

194 

0. 

981 

13. 

179 

17  . 

268 

2. 

200 

1. 

348 

0. 

194 

0. 

981 

13. 

179 

17  . 

268 

2. 

200 

1. 

593 

0. 

194 

0. 

981 

13. 

179 

17  . 

2 68. 

2. 

200 

1. 

831 

0. 

194 

0. 

981 

13. 

179 

17. 

268 

2. 

200 

2. 

148 

0. 

194 

0. 

981 

13. 

179 

17  . 

268 

2. 

200 

2. 

421 

0. 

189 

0. 

617 

8. 

795 

12. 

087 

1. 

676 

2. 

606 

0. 

189 

0. 

617 

8. 

795 

12. 

087 

1. 

676 

2. 

726 

0. 

189 

0. 

617 

8. 

795 

12. 

087 

1. 

67  6 

2. 

906 

0. 

189 

0. 

617 

8. 

795 

12. 

087 

1. 

676 

Runoff  Event  Related  Data 

For  each  requested  runoff  event  CONCEPTS  writes  the  storm  date  and 
requested  profiles  to  their  output  files.  The  first  line  in  the  output  file  for  each 
runoff  event  is  the  date  of  the  storm  event: 


example  storm 
date  output  line 


OUTPUT  FOR  THE  STORM  EVENT  STARTED  ON  12/25/1982  11:58:  0 

The  layout  further  consists  of  a header  followed  by  multiple  columns  of  which 
the  first  column  is  river  kilometer.  Following  is  a list  of  possible  columns  in 
the  header: 


• river  kilometer, 


RIVER  KM 
(KM) 

• peak  discharge  (data  type  l), 


DISCHARGE 

(CMS) 

• peak  stage  (data  type  2), 


STAGE 

(M) 


• thalweg  elevation  (data  type  4), 


THALWEG 

(M) 


• change  in  bed  elevation  (data  type  8) 


BED  CHNG 
(M) 
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• in-bank  top  width  of  cross  section  (data  type  16), 


TOP  WIDTH 
(M) 


• height  of  left  and  right  bank  of  cross  section  (data  type  32), 

BANK  HEIGHT 
LEFT  RIGHT 

(M)  (M) 


• storm  event  generated  sediment  yield  (data  type  64),  and 

STORM  EVENT  GENERATED  SEDIMENT  YIELD 
SILT  YLD  SAND  YLD  GRAVEL  YLD  TOTAL  YLD 
(TONS)  (TONS)  (TONS)  (TONS) 

• characteristic  sediment  particle  sizes  of  surface  layer  of  streambed  (data  type 
128). 

BED  MATERIAL  CHARACTERISTIC  DIAMETERS 
D16  D50  D84  D90  DMEAN 

(MM)  (MM)  (MM)  (MM)  (MM) 


Processing  of  Output  Data  Files 

CONCEPTS  does  not  produce  graphical  output,  such  as  plots  of  hydrographs. 
However,  the  data  columns  in  the  output  files  have  a fixed  width  of  1 0 
characters,  hence  it  is  easy  to  import  them  as  text  files  into  ‘spreadsheet’ 
applications  such  as  Microsoft®  Excel  or  Corel®  Quattro®  Pro.  Below  we 
will  show  how  you  can  produce  graphs  from  the  data  in  the  three  types  of 
output  files  using  Microsoft®  Excel  2000. 

Output  at  a Certain  Location  and  for  a Certain  Runoff 
Event 


You  can  produce  charts  of  cross-sectional  geometry,  sediment  particle  size 
distribution,  sediment  yield,  etc.  using  the  outputted  data  at  a certain  location 
and  for  a certain  runoff  event.  This  section  outlines  the  steps  you  need  to  take 
to  import  the  data  into  Microsoft®  Excel.  You  should  refer  to  the  Microsoft® 
Excel  manuals  or  online  help  for  additional  support. 

1 Start  Microsoft®  Excel. 

2 Select  “File>Open”  and  then  “All  Files”  in  the  “Files  of  type:”  select  box. 

3 Navigate  to  the  folder  that  contains  the  output  file.  Select  this  file  and  open 
it.  This  will  bring  up  the  “Text  Import  Wizard  - Step  1 of  3”  window. 
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Figure  5.1  “Text  Import  Wizard  - Step  1 of  3”  window. 


4 Select  “Fixed  width”  as  the  “Original  data  type”  and  click  on  the  “Next  >” 
button.  This  will  bring  up  the  “Text  Import  Wizard  - Step  2 of  3”  window. 

Figure  5.2  “Text  Import  Wizard  - Step  2 of  3”  window. 


5 Scroll  down  to  get  a proper  view  of  the  data  you  will  import.  Create  break 
lines  at  10,  20,  30,  40,  etc.  depending  on  the  output  data  type. 
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Figure  5.3  “Text  Import  Wizard  - Step  2 of  3”  window  with  break  lines  added. 


6 Click  on  the  “Finish”  button.  You  have  imported  the  output  data  in 
Microsoft®  Excel  (see  Figure  5.4). 

Now  you  can  plot  the  data  using  the  graphing  capabilities  of  Microsoft® 
Excel.  Figures  5.6  and  5.7  show  examples  of  temporal  changes  in  cross-section 
geometry  and  bed-material  composition. 

Time-Series  Output  at  a Certain  Location 

You  can  produce  graphs  showing  the  progression  of  discharge,  stage,  thalweg 
elevation,  factor  of  safety,  etc.  (see  Table  4.8)  using  the  outputted  time-series 
data  at  a certain  location.  This  section  outlines  the  steps  you  have  to  take  to 
import  the  data  into  Microsoft®  Excel. 

To  import  the  time-series  output  text  files  you  need  to  take  the  steps  presented 
in  the  previous  section  “Output  at  a Certain  Location  and  for  a Certain  Runoff 
Event”.  However  in  step  5,  you  need  to  place  the  first  break  line  at  position  20 
because  the  time  column  has  a width  of  20  characters.  The  next  break  lines 
need  to  be  placed  at  positions  30,  40,  50,  etc.  (see  Figure  5.5). 

Figure  5.8  shows  the  imported  output  file  in  Microsoft®  Excel.  Figure  5.9 
shows  an  XY  plot  of  factor  of  safety  of  the  right  streambank  against  time. 

Output  for  a Certain  Runoff  Event  along  a Section  of  the 
Modeling  Reach 


You  can  produce  graphs  showing  the  thalweg  profile,  distribution  of  £/50  along 
the  channel,  etc.  (see  Table  4.10)  at  various  points  of  time  using  the  profile 
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Figure  5.4  View  of  Microsoft®  Excel  spreadsheet  after  importing  output  data  file. 
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Figure  5.5  “Text  Import  Wizard  - Step  2 of  3”  window  with  break  lines  added. 
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Figure  5.6  Example  plot  of  temporal  changes  in  cross-section  geometry. 
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Figure  5.7  Example  plot  of  temporal  changes  in  bed-material  composition. 
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Figure  5.8  View  of  Microsoft®  Excel  spreadsheet  after  importing  time-series  output 
data  file. 
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Figure  5.9  Example  plot  of  the  progression  of  factor  of  safety  of  the  right  streambank. 
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Figure  5.10  “Text  Import  Wizard  - Step  2 of  3”  window  with  break  lines  added. 


output.  This  section  outlines  the  steps  you  have  to  take  to  import  the  data  into 
Microsoft®  Excel. 

To  import  the  profile  output  text  files  you  need  to  take  the  steps  presented  in 
the  prior  section  “Output  at  a Certain  Location  and  for  a Certain  Runoff 
Event”.  In  step  5,  you  need  to  place  break  lines  at  every  10th  position,  that  is 
10,  20,  30,  40,  50,  etc.  (see  Figure  5.10). 

Figure  5.1 1 shows  the  imported  output  file  in  Microsoft®  Excel.  Figure  5.12 
shows  an  example  plot  of  temporal  changes  in  thalweg  profile. 
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Figure  5.11  View  of  Microsoft®  Excel  spreadsheet  after  importing  profile  output  data 
file. 


E2  Microsoft  Excel 

- Goodwin  68.txt 

S]  File  Edit  View 

Insert  Format  Tools  Data  Window  Help 

D 13  [r^l 

iay  ife © «-> * 

\la  f"  zt 

U 100% 

- (?)  ” Arial 

- 10  E 

/ u m w m td°  .-to 

■S»  - “ Prompt  ” 

A1 

= Output  pr 

A 

B C 

D 

E 

F 

G 

H 

i 

J K 

L M — 

Output  pr 

ofile  file  from  cros 

s section 

"Goodwin  C 

reek  cross 

section  T 

7A"  at  riv 

er  mile 

248  1 meters 

— 

2 

3 

to  cross 

section  "G  oodwin  Cre 

ek  cross  s 

ection  C4- 

8"  at  rive 

r mile  29 

05.5  meter 

s 

4 

Output  ty 

pe  = 255 

5 

Output  da 

tes  are  1 2/27/82 

16  00  00 

6 

1 2/27/83 

20  00  00 

7 

1 2/5/84 

17  00  00 

8 

1 2/15/85 

15  45:00 

9 

1 2/23/86 

7 00:00 

10 

1 2/31/1987 

16  00  00 

- 

11 

1 2/31/1988 

11  00  00 

12 

1 2/30/1989 

22  30  00 

13 

1 2/30/1990 

14  30.00 

14 

0 9/9/91 

13  15  00 

a 

15 

o a 

16 

1 

17 

1 

18 

Til 

19 

I 

20 

INITIAL  D 

ATAAT 10/  17/198115 

00  00 

i 

21 

■ 

22 

1 

23 

BANK 

HEIGHT 

BED 

MATERIAL 

CHARACTERI 

STIC  DIAME 

TERS 

24 

RIVER  KM 

THALWEG  TOP  WIDTH 

LEFT 

RIGHT 

D16 

D50 

D84 

D90 

DMEAN 

3 

25 

(KM) 

(M)  (M) 

(M) 

(M) 

(MM) 

(MM) 

(MM) 

(MM) 

(MM) 

3 

26 

0 248 

73  27  88  6 

2.45 

3 42 

0 227 

3 317 

1674 

20  889 

3 245 

n 

‘27 

0 306 

72  9 37 

3 95 

4 1 

0 227 

3 317 

16  74 

20  889 

3 245 

3 

28 

0 364 

72  48  59.71 

4 45 

2.1 

0.227 

3317 

16  74 

20  389 

3.245 

3 

*291 

0 432 

72  3 70 

4 45 

4 

0 227 

3 317 

16  74 

20  889 

3 245 

3 

30 

0 499 

72.11  78  97 

2.35 

2.5 

0 227 

3 317 

16  74 

20  889 

3245 

3 

31 

0 676 

71.65  54  89 

1.49 

4 02 

0.227 

3 317 

16  74 

20  889 

3 245 

9 

32 

0 819 

71.08  27,74 

3 84 

4 51 

0 227 

3.317 

16  74 

20  889 

3 245 

if 

Figure  5.12  Example  plot  of  temporal  changes  in  thalweg  profile. 
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Goodwin  Creek,  Mississippi 

This  chapter  presents  an  application  of  CONCEPTS  to  the  lower 
3.5  km  of  Goodwin  Creek,  Mississippi.  You  will  find  the  physical 
characteristics  of  the  channel,  the  input  data  necessary  to  perform 
the  simulation,  and  results  showing  the  effects  of  streambank 
erosion  on  the  channel’s  morphology. 


Study  Area 

Goodwin  Creek  is  located  in  north-central  Mississippi;  its  watershed  is  situated 
on  the  eastern  side  of  the  Bluff  or  Loess  Hills  physiographic  subprovince  (see 
Figure  6.1).  In  the  upper  reaches  of  the  watershed,  streamflow  only  occurs 
immediately  after  precipitation.  In  the  lower  reaches,  a slight  baseflow  usually 
persists  for  most  of  the  year.  Storm  flows  carry  heavy  loads  of  coarse  sediment 
(sand  and  gravel)  derived  from  upland  tributaries  and  gullies  and  from  valley 
alluvium.  The  stream  channels  are  greatly  enlarged  in  places,  are  highly 
variable  in  form  and  are  deeply  incised.  Terrain  elevation  ranges  from  71  to 
128  m above  mean  sea  level,  with  an  average  channel  slope  of  0.004. 

The  Goodwin  Creek  watershed  is  highly  instrumented  to  provide  data  needed 
to  investigate  the  impact  of  landscape  attributes  and  watershed  processes  on 
sediment  yield,  to  test  concepts  developed  in  the  study  of  these  processes,  and 
to  validate  models  developed  in  the  research  (Alonso,  1997).  The  US 
Department  of  Agriculture- Agricultural  Research  Service-National 
Sedimentation  Laboratory  (NSL)  operates  the  watershed.  The  watershed  is 
divided  into  fourteen  nested  subcatchments  with  a flow  measuring  flume 
constructed  at  each  of  the  drainage  outlets.  Twenty  nine  standard  recording 
rain  gages  are  uniformly  located  within  and  just  outside  the  watershed.  The 
average  annual  rainfall  during  1982-1992  from  all  storms  was  1440  mm  and 

the  mean  annual  runoff  measured  at  the  watershed  outlet  was  14xl06  nr\ 
Measuring  flumes  are  designed  to  operate  in  the  supercritical  regime  to  prevent 
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Figure  6.1  Map  of  the  Goodwin  Creek  Watershed. 


MEMPHIS 


deposition  of  sediment  in  the  flumes.  Instrumentation  at  each  gaging  site 
includes  an  electronic  data  acquisition  and  radio  telemetry  system  that  collects, 
stores  and  transmits  the  data  to  a central  computer  at  the  NS L for  processing 
and  archival.  Measurements  collected  at  each  site  and  transmitted  through  the 
telemetry  system  include  water  stage,  accounting  of  automatically  pumped 
sediment  samples,  air  and  water  temperature,  precipitation,  and  climatological 
parameters.  Sampling  of  total  sediment  loads  is  carried  out  during  storm  events 
at  selected  stations  using  Helley-Smith  bedload  samplers  and  DH-48  depth- 
integrating  suspended  sediment  samplers.  NSL  has  made  available  a 
comprehensive  documentation  of  the  Goodwin  Creek  Experimental  Watershed 
and  the  database  compiled  for  the  period  1 982- 1 993  for  downloading  from  the 
NSL  computer  system  at: 

http://www.sedlab.olemiss.edu/cwD  unit/qcwd  ftp.html 

The  US  Army  Corps  of  Engineers  prepared  a detailed  2-ft  contour  interval 
topographic  survey  of  all  the  main  channel  and  primary  tributaries  of  Goodwin 
Creek  in  1977.  In  1982  NSL  resurveyed  30  of  the  1977  cross  sections  in  the 
lower  3.86  km  of  Goodwin  Creek  (see  Figure  6.2;  and  Murphey  and 
Grissinger,  1986).  NSL  has  repeated  these  cross  section  surveys  26  times,  and 
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Figure  6.2  Plan  map  of  the  3.86  km  study  length  of  Goodwin  Creek  showing  cross 
sections  and  flume  locations. 


used  them  to  document  changes  in  channel  geometry  and  to  verify  the 
sediment  delivery  budget  of  the  channel  system. 

The  bed  material  was  sampled  in  1994  (Kuhnle,  1996).  The  main  channels 
were  divided  into  a series  of  reaches  with  a mean  length  of  1 ,000  m.  Figure  6.3 
shows  the  median  size  of  the  bed  material  of  14  sampled  reaches.  Median  sizes 
ranged  from  0.5  to  7.4  mm.  They  are  generally  fine  in  the  upstream  portions  of 
the  watershed,  coarsen  toward  the  central  part  of  the  watershed,  and  become 
finer  again  towards  the  outlet  of  the  watershed. 

In  the  uplands  three  lithologic  units  overly  the  Citronelle  deposits  of  sand  and 
gravel:  (a)  an  early-Holocene  sequence  of  channel  lag,  organic  bog,  and  gray 
silt  overlain  by  a massive  silt  deposit,  (b)  a mid-Holocene  deposit  of  well- 
bedded  relatively  coarse-textured  channel  fill,  and  (c)  a late-Holocene  deposit 
of  fine  to  coarse  meander-belt  alluvium  (Murphey  and  Grissinger,  1986). 
Superimposed  over  these  three  deposits  are  massive  amounts  of  finely-layered 
silt  and  fine  sand  that  were  washed  from  the  hills  following  European 
settlement.  Simon  and  Darby  (1997)  determined  the  following  geotechnical 
characteristics  of  the  bank  material: 
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Figure  6.3  The  14  reaches  along  Goodwin  Creek  for  which  bed-material  samples 

were  collected.  Bed  material  median  sizes  for  the  reaches  are:  1)  1 .0  mm, 
2)  1 .4  mm,  3)  4.4  mm,  4)  6.9  mm,  5)  3.5  mm,  6)  6.7  mm,  7)  3.9  mm,  8)  1 .1 
mm,  9)  1 .4  mm,  10)  7.3  mm,  1 1)  1 .4  mm,  12)  0.5  mm,  13)  2.4  mm,  and  14) 
7.4  mm. 


• Early-Holocene  unit:  ca  « 5 kPa  and  <j)'  « 35° . 

• Late-Holocene  unit:  ca  ~ 35  kPa  and  <j)'  « 28° . 

Application 

We  have  applied  CONCEPTS  to  the  lower  3.5  km  of  Goodwin  Creek  between 
flumes  1 and  2 (Figure  6.2)  for  water  years  1982  through  1991 . We  have  used 
the  data  presented  in  the  previous  section.  The  enclosed  CDROM  contains  the 
input  files  to  perform  the  simulation. 

The  modeling  reach  extends  from  cross  section  T7A-T5  to  cross  section  C4-8 
with  a length  of  2.7  km  (Figure  6.2),  and  contains  22  surveyed  cross  sections. 
However,  in  this  configuration  numerical  problems  occurred  near  the  inlet 
(upstream  boundary,  cross  section  T7A-T5)  when  discharge  at  the  inlet 
increased  rapidly  causing  the  flow  depth  at  the  next  cross  section  to  drop 
significantly.  We  solved  this  problem  by  reducing  the  space  steps  near  the 
inlet;  that  is,  inserting  a cross  section  halfway  between  cross  sections  T7A-T5 
and  T5-T5A  and  a cross  section  halfway  between  cross  sections  T5-T5A  and 
T3-5.  The  geometries  of  the  inserted  cross  sections  assure  a smooth  transition 
from  the  surveyed  cross  section  upstream  to  that  downstream  of  the  inserted 
cross  sections.  The  simulation  does  not  take  into  account  the  effects  of  flume  1 
on  channel  hydraulics  and  morphology.  We  fixed  the  thalweg  elevation  at  the 
downstream  boundary  of  the  modeling  reach.  We  imposed  the  discharges 
measured  at  flume  2 at  the  upstream  boundary  of  the  modeling  reach. 

We  have  simulated  the  response  of  the  modeling  reach  to  the  following  two 
scenarios: 

1 Channel  can  only  adjust  vertically.  We  only  turned  on  the  sediment  transport 
submodel  of  CONCEPTS,  which  includes  variations  in  streambed  elevation. 
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2 Channel  can  adjust  both  vertically  and  horizontally.  CONCEPTS  simulates 
both  streambed  and  streambank  mechanics. 

Figure  6.4  shows  the  Run  Control  Data  File , Figure  6.5  shows  an  example  of  a 
Cross  Section  Data  File , and  Figure  6.6  shows  part  of  the  Discharge  Data  File. 

Figure  6.4  Run  control  data  file  used  in  the  simulation  of  the  lower  3.5  km  of  Goodwin 
Creek,  Mississippi. 

Main  Input  File 


case  name 
Goodwin 

! project  title 

Goodwin  Creek  channel  evolution  simulation  between  1982  and  1991 

! Run  Control  Data 

! upstream  flow  discharge  file 
GChydrography . txt 
! flow-related  data 
0.25  0.00  0 

! Sediment  discharge  at  upstream  end  of  the  channel 
0 0.00  0.00  0.00  0.975  0.975  0.99  0.99  1.00  1.00  1.00  1.00  1.00  1.00 

! Silt  Fraction  and  downstream  bed  control 
0.30  1.0 

! bank  failure  analysis 

0 4 

! type  of  flow  resistance  formulation 
1 

! water  temperature 

10.0 

! sediment  and  streambank  mechanics  options 

1 1 1 

i Simulation  Times  

! start  end  time  step 

10/17/1981  15:00:00  09/31/1991  23:59:00  100 

! Makeup  of  Modeling  Reach  

! number  of  links 
1 

! linktypes  for  the  above  number  of  links 
1 

i Link  x 

! REACH:  number  of  cross  sections  and  their  data  filenames 
24 

GCxsectOl . txt 
GCxsectOla . txt 
GCxsect02 . txt 
GCxsect02a . txt 
GCxsect03 . txt 
GCxsect04 . txt 
GCxsect05 . txt 
GCxsect06 . txt 
GCxsect07 . txt 
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Figure  6.4  Run  control  data  file  used  in  the  simulation  of  the  lower  3.5  km  of  Goodwin 
Creek,  Mississippi,  (continued) 


GCxsect08 

GCxsect09 

GCxsectlO 

GCxsectll 

GCxsectl2 

GCxsectl3 


txt 

txt 

txt 

txt 

txt 

txt 


GCxsectl4 

GCxsectl5 

GCxsectl6 

GCxsectl7 

GCxsectl8 

GCxsectl9 


txt 

txt 

txt 

txt 

txt 

txt 


GCxsect20 

GCxsect21 


txt 

txt 


GCxsect22 


txt 


! single  point  and  time 
1 

17536 
1 12 
5 

12/27/1983  20:00:00 
12/15/1985  15:45:00 
12/31/1987  16:00:00 
12/30/1989  22:30:00 
09/09/1991  13:15:00 


Output 


! single  point,  continuously  in  time 


7 

595973 
1 1 
1 

10/01/1981  00:00:00 
595973 
1 5 

1 

10/01/1981  00:00:00 
1 9 

1 

10/01/1981  00:00:00 
595973 
1 13 

1 

10/01/1981  00:00:00 
1 17 

1 

10/01/1981  00:00:00 
595973 
1 21 
1 

10/01/1981  00:00:00 


12/31/1991  24:00:00 

12/31/1991  24:00:00 
12/31/1991  24:00:00 

12/31/1991  24:00:00 
12/31/1991  24:00:00 

12/31/1991  24:00:00 
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Figure  6.4  Run  control  data  file  used  in  the  simulation  of  the  lower  3.5  km  of  Goodwin 
Creek,  Mississippi,  (continued) 

595973 
1 24 

1 

10/01/1981  00:00:00  12/31/1991  24:00:00 

! profile  at  specific  time  points 
1 

255 

1 1 1 24 

5 

12/27/1983  20:00:00 
12/15/1985  15:45:00 
12/31/1987  16:00:00 
12/30/1989  22:30:00 
09/09/1991  13:15:00 


Figure  6.5  A sample  cross  section  data  file  used  in  the  simulation  of  the  lower  3.5  km 
of  Goodwin  Creek,  Mississippi. 


! Input  file  of  cross  section  14. 

i 

! Name  of  xsection  and  rivermile  in  (km) 

Goodwin  Creek  cross  section  C10-1 
2.1477 

! friction  factor 
0.06 

! discharge  fraction  due  to  tributary  inflow 

0.00 

i Left  FloodPlain 

! number  of  nodes 
6 

! station  and  elevation  for  above  number  of  coordinates  in  (m) 


0. 

.00 

74. 

.21 

1. 

.00 

73. 

.21 

98. 

.63 

73. 

.21 

100. 

.00 

73. 

.21 

100. 

.27 

73. 

.21 

102. 

.07 

73. 

.21 

! Manning's  n of  left  floodplain 

0.10 

! Left  Bank 

! number  of  nodes 
2 

! station  and  elevation  for  above  number  of  coordinates  in  (m) 
102.07  73.21 

106.77  69.55 

! bank  material 

9830.0  17.1  15.0  21000.0 

! bank  sediment  composition 
30.700 
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Figure  6.5  A sample  cross  section  data  file  used  in  the  simulation  of  the  lower  3.5  km 
of  Goodwin  Creek,  Mississippi,  (continued) 

37.300 

19.300 
1.600 

4.400 

1.400 
1.200 
2.500 
1.600 
0.000 
0.000 
0.000 
0.000 

! critical  shear  stress  for  erosion 
4.0 

! groundwatertable,  hydraulic  gradient,  and  angle  of  seepage  force 
70.55  1.00  0.00 

! Manning ' s n 
0.06 

I Channel  Bed 

! number  of  nodes 
7 

! station  and  elevation  for  above  number  of  coordinates  in  (m) 


106. 

.77 

69. 

.55 

109. 

.20 

69. 

.16 

120. 

.76 

69. 

. 15 

124  . 

.26 

69. 

.09 

125. 

.18 

68. 

.73 

130. 

.17 

68. 

.49 

135. 

.11 

68. 

.12 

! Elevation  of  bedrock  (m) 

0.00 

! porosity 
0.4 

! hiding  factors 
0.00  0.95  0.70 

! Surface  layer  and  Substratum  data 

! number  of  sediment  layers  composing  the  bed 
1 

! Layer  1,  layer  depth  below  bed  surface 

0.00 

! bed  composition 
0.000 
0.000 
0.000 
2.684 
40.237 
9.675 
4.464 
5.939 
8.881 


1 16  CONCEPTS  Manual 


Application 


Figure  6.5  A sample  cross  section  data  file  used  in  the  simulation  of  the  lower  3.5  km 
of  Goodwin  Creek,  Mississippi,  (continued) 

10.656 

11.541 

5.043 

0.880 

! critical  shear  stresses  for  erosion  of  and  deposition  on  cohesive 
beds 

! and  erodibility  coefficient 
0.10  4.00  5.00E-06 

! Manning's  n 
0.04 

i Right  Bank 

! number  of  nodes 
6 

! station  and  elevation  for  above  number  of  coordinates  in  (m) 


135. 

.11 

68. 

. 12 

135. 

.42 

68. 

.45 

136. 

.36 

68. 

.85 

138. 

.25 

69. 

.31 

139. 

.10 

70. 

.13 

142. 

.18 

70. 

.77 

! bank  material 

9830.0  17.1  15.0  21000.0 

! bank  sediment  composition 

30.700 

37.300 

19.300 
1.600 

4.400 

1.400 
1.200 
2.500 
1.600 
0.000 
0.000 
0.000 
0.000 

! critical  shear  stress  for  erosion 

4.0 

! groundwatertable,  hydraulic  gradient,  and  angle  of  seepage  force 
69.12  1.00  0.00 

! Manning's  n 
0.06 

! Right  FloodPlain 

! number  of  nodes 
10 

! station  and  elevation  for  above  number  of  coordinates  in  (m) 
142.18  70.77 

146.11  70.77 

147.88  70.77 

152.12  72.08 
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Figure  6.5  A sample  cross  section  data  file  used  in  the  simulation  of  the  lower  3.5  km 
of  Goodwin  Creek,  Mississippi,  (continued) 


153.58 

73.24 

154.37 

73.55 

154.98 

73.64 

156.08 

73.64 

256.08 

73.64 

257.08 

74.64 

! Manning's 
0.10 

n of  right  floodplain 

Figure  6.6  Part  of  the  discharge  data  file  used  in  the  simulation  of  the  lower  3.5  km  of 
Goodwin  Creek,  Mississippi. 


11/30/1981 

17 

40 

00 

0 

0 

11 

11/30/1981 

18 

10 

00 

1 

0 

20 

11/30/1981 

18 

30 

00 

0 

0 

28 

11/30/1981 

18 

40 

00 

0 

0 

34 

11/30/1981 

19 

00 

00 

0 

0 

48 

11/30/1981 

19 

08 

00 

0 

0 

63 

11/30/1981 

19 

16 

00 

0 

0 

95 

11/30/1981 

19 

20 

00 

0 

1 

17 

11/30/1981 

19 

24 

00 

0 

1 

54 

11/30/1981 

19 

26 

00 

0 

1 

78 

11/30/1981 

19 

28 

00 

0 

2 

03 

11/30/1981 

19 

32 

00 

0 

2 

20 

11/30/1981 

19 

40 

00 

0 

2 

63 

11/30/1981 

19 

52 

00 

0 

3 

18 

11/30/1981 

20 

00 

00 

0 

3 

59 

11/30/1981 

20 

10 

00 

0 

3 

62 

11/30/1981 

20 

20 

00 

0 

3 

67 

11/30/1981 

20 

30 

00 

0 

3 

67 

11/30/1981 

20 

36 

00 

0 

3 

62 

11/30/1981 

20 

48 

00 

0 

3 

54 

11/30/1981 

21 

00 

00 

0 

3 

50 

11/30/1981 

21 

11 

00 

0 

3 

30 

11/30/1981 

21 

20 

00 

0 

3 

14 

11/30/1981 

21 

22 

00 

0 

3 

11 

11/30/1981 

21 

34 

00 

0 

2 

85 

11/30/1981 

21 

44 

00 

0 

2 

64 

11/30/1981 

21 

56 

00 

0 

2 

41 

11/30/1981 

22 

00 

00 

0 

2 

34 

11/30/1981 

22 

07 

00 

0 

2 

20 

11/30/1981 

22 

20 

00 

0 

1 

96 

11/30/1981 

22 

30 

00 

0 

1 

79 

11/30/1981 

22 

42 

00 

0 

1 

59 

11/30/1981 

22 

54 

00 

0 

1 

41 

11/30/1981 

23 

00 

00 

0 

1 

33 

11/30/1981 

23 

04 

00 

0 

1 

29 

11/30/1981 

23 

16 

00 

0 

1 

16 

11/30/1981 

23 

26 

00 

0 

1 

07 

11/30/1981 

23 

30 

00 

0 

1 

03 
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Figure  6.6  Part  of  the  discharge  data  file  used  in  the  simulation  of  the  lower  3.5  km  of 
Goodwin  Creek,  Mississippi,  (continued) 


11/30/1981 

23 

38 

00 

0 

0 

96 

11/30/1981 

23 

48 

00 

0 

0 

86 

12/01/1981 

00 

00 

00 

0 

0 

75 

12/01/1981 

00 

01 

00 

0 

0 

73 

12/01/1981 

00 

12 

00 

0 

0 

69 

12/01/1981 

00 

20 

00 

0 

0 

66 

12/01/1981 

00 

22 

00 

0 

0 

66 

12/01/1981 

00 

29 

00 

0 

0 

62 

12/01/1981 

00 

36 

00 

0 

0 

59 

12/01/1981 

00 

52 

00 

0 

0 

53 

12/01/1981 

01 

02 

00 

0 

0 

49 

12/01/1981 

01 

10 

00 

0 

0 

47 

12/01/1981 

01 

21 

00 

0 

0 

44 

12/01/1981 

01 

34 

00 

0 

0 

42 

12/01/1981 

01 

49 

00 

0 

0 

38 

12/01/1981 

02 

16 

00 

2 

0 

33 

12/01/1981 

03 

10 

00 

0 

0 

26 

12/01/1981 

03 

30 

00 

0 

0 

24 

12/01/1981 

03 

47 

00 

0 

0 

22 

12/01/1981 

05 

08 

00 

0 

0 

16 

12/01/1981 

05 

32 

00 

0 

0 

14 

12/01/1981 

06 

58 

00 

0 

0 

11 

01/02/1982 

20 

44 

00 

0 

0 

11 

01/02/1982 

20 

45 

00 

0 

0 

13 

01/02/1982 

20 

47 

00 

0 

0 

15 

01/02/1982 

20 

49 

00 

1 

0 

18 

01/02/1982 

20 

52 

00 

0 

0 

24 

01/02/1982 

20 

54 

00 

0 

0 

34 

01/02/1982 

20 

55 

00 

0 

0 

43 

01/02/1982 

20 

56 

00 

0 

0 

56 

01/02/1982 

20 

58 

00 

0 

0 

76 

01/02/1982 

21 

00 

00 

0 

1 

05 

01/02/1982 

21 

02 

00 

0 

1 

29 

01/02/1982 

21 

03 

00 

0 

1 

39 

01/02/1982 

21 

04 

00 

0 

1 

53 

Hydraulics 


Figure  6.7  shows  changes  in  the  shape  of  the  rating  curve  at  the  outlet  of  the 
modeling  reach  during  1982  and  at  the  end  of  the  simulation  period.  We 
obtained  these  results  for  the  second  scenario;  however,  the  predicted  changes 
in  rating-curve  shape  are  similar  for  scenario  one.  Figure  6.7  shows  that: 

• flow  depth  at  base  flow  increases; 

• average  slope  of  the  rating  curve,  that  is  dQ/dh , decreases;  and 

• difference  in  flow  depths  for  the  same  discharge  on  the  rising  and  falling 
limb  of  the  hydrograph  increases. 
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Figure  6.7  Simulated  storm  rating-curves  at  the  outlet  of  the  modeling  reach. 
35 
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We  can  relate  the  above  points  to  adjustments  of  channel  geometry,  which  is 
discussed  in  the  next  section  on  “Morphology”.  The  changes  in  shape  of  the 
rating  curve  and  the  larger  flow  depth  at  base  flow  infer  that  channel  slope  has 
reduced  over  the  simulation  period.  The  types  of  flood  waves  propagating 
through  the  channel  have  changed  from  diffusion  to  dynamic  waves. 

The  mean  slope  of  the  rating  curve,  dQ/dh  , of  the  stonn  event  occurring  on 

October  18,  1981,  is  16.0  m2/s,  whereas  it  is  5.5  m2/s  for  the  storm  event 
occurring  on  April  27,  1991.  Assuming  that  changes  in  the  derivative  of 
conveyance  with  respect  to  flow  depth,  dK/ dh , are  negligible,  the  channel 
slope  has  reduced  by  a factor  of  8.5.  Figure  6.8  shows  that  channel  slope  has 
approximately  reduced  from  0.0012  to  0.0002  along  the  downstream  most  0.5 
km.  This  agrees  well  with  the  change  in  mean  slope  of  the  rating  curve. 
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Figure  6.8  Simulated  changes  in  profiles  of:  thalweg  elevation  ((a)  scenario  1 and  (b) 
scenario  2)  and  d50  ((c)  scenario  1 and  (d)  scenario  2). 
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Morphology 

Figure  6.8  shows  simulated  profiles  of  thalweg  elevation  and  d50  for  the  two 
scenarios.  Changes  in  the  thalweg  profile  occur  mainly  in  the  first  two  years  of 
the  simulation  period.  This  agrees  well  with  the  changes  in  the  shape  of  the 
rating  curve,  which  mainly  occur  in  the  first  year  of  the  simulation  period  (see 
Figure  6.7).  The  modeling  reach  incises  more  in  the  case  where  it  can  adjust 
only  vertically  (scenario  1,  Figure  6.8a)  than  in  the  case  where  it  can  adjust 
both  vertically  and  horizontally  (scenario  2,  Figure  6.8b).  At  the  end  of  the 
simulation  period  the  modeling  reach  continues  to  adjust  for  scenario  1 (Figure 
6.8a),  whereas  it  has  reach  a quasi-dynamic  equilibrium  for  scenario  2 (Figure 
6.8b). 

The  composition  of  the  bed  surface  is  coarser  for  scenario  1 (Figure  6.8c)  than 
for  scenario  2 (Figure  6.8d),  because: 

1 Incision  is  larger  for  scenario  1 than  scenario  2.  Hence,  more  ‘finer’ 
particles  are  removed  from  the  bed  in  scenario  1 . 

2 Bank  failures  occurring  in  scenario  2 add  fine  material  to  the  sediment 
transport. 

Figure  6.9  shows  simulated  changes  in  geometry  of  cross  section  C45-1 . The 
cross  section  is  incising.  There  is  some  minor  deposition  on  the  bar  near  the 


121 


Goodwin  Creek,  Mississippi 


Figure  6.9  Simulated  changes  in  geometry  of  cross  section  C45-1 : (a)  scenario  1 and 
(b)  scenario  2. 
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right  bank.  In  scenario  2,  three  left-bank  failures  occur,  mainly  because  the 
flow  removes  bank  material  from  the  bank  toe. 
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Goodwin  Creek  Bendway,  Mississippi 


This  chapter  presents  the  application  of  CONCEPTS  to  a bendway 
in  Goodwin  Creek,  Mississippi.  You  will  find  the  physical 
characteristics  of  the  bendway,  the  input  data  necessary  to  perform 
the  simulation,  and  results  showing  the  capability  of  CONCEPTS 
to  accurately  predict  streambank  erosion. 


Study  Area 

Since  February  1996  the  US  Department  of  Agriculture- Agricultural  Research 
Service-National  Sedimentation  Laboratory  (NSL)  has  conducted  extensive 
research  on  streambank  failure  mechanics  along  a bendway  of  Goodwin  Creek, 
northern  Mississippi  (Simon  et  al.,  1999).  The  following  data  are  being 
collected: 

• cross  section  geometry, 

• water  surface  elevations, 

• bank-material  properties, 

• bank-material  shear- strength  parameters, 

• pore- water  pressures  in  the  bank, 

• root  mapping  and  tensile  strength,  and 

• plant  stem  flow. 

Two  flow  measuring  flumes  (flumes  #3  and  #4,  see  Figure  6.1)  in  upstream 
tributaries  provide  continuous  discharge  and  fine  sediment  data.  Figures  7.1 
and  7.2  show  a photo  and  a plan  view  of  the  bendway  with  locations  of  ten 
surveyed  cross  sections,  respectively.  The  flow  is  from  top  to  bottom. 

Bank  material  consists  of  about  2 m of  brown,  clayey-silt  of  late  Holocene  age 
(LH  unit)  overlying  1 .5  m of  early  Holocene  gray,  blocky  silt  of  lower 
permeability  (EH  unit).  These  units  are  separated  by  a thin  (0. 1 to  0.2  m)  layer 
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Figure  7.1  Photo  of  the  Goodwin  Creek  Bendway  study  site.  Flow  is  from  right  to  left. 


containing  manganese  nodules  and  characterized  by  very  low  permeability, 
which  perches  water.  These  materials  overlie  1 m of  sand  and  1 .5  m of  packed 
sandy  gravel. 

Apparent  cohesion  and  effective  angle  of  internal  friction  were  measured  in 
situ  using  an  Iowa  Borehole  Shear  Tester  (Luttenegger  and  Hallberg,  1981). 
For  the  LH  unit,  results  revealed  an  apparent  cohesion  varying  between  0.0  and 
8.4  kPa  with  an  average  friction  angle  of  28.1°.  The  underlying  EH  unit  has 
apparent  cohesion  values  varying  between  0.0  and  37.9  kPa  with  an  average 
friction  angle  of  28.5°  for  the  period  between  June  1996  and  July  1998  (Simon 
et  al .,  1999).  Simon  et  al.  (1999)  provide  the  following  values  for  effective 

cohesion  and  angle  §b  : 

• LH  unit,  c'  = 2.67  kPa  and  (j/  - 10.4°;  and 

• EH  unit,  c'  = 13.0  kPa  and  <\>b  = 17.5°. 

Four  major  failure  episodes  occurred  at  the  research  site  between  February  and 
December  1 996,  resulting  in  up  to  2 m of  top-bank  retreat.  This  rate  is  greater 
than  the  30-year  average  of  about  0.5  m/yr,  and  is  attributed  to  persistent 
precipitation,  manifest  in  1 0 discharge  peaks  with  a 1 -year  recurrence  interval 
or  greater  (Simon  and  Darby,  1997).  Planar  and  cantilever  failures  were 
relatively  common  along  the  steepest  section  of  the  4.7  m high  banks. 
Cantilevers  were  formed  by: 

1 preferential  erosion  of  sands  and  silts  by  fluvial  undercutting  about  3.0  to 
3.5  m below  the  top  bank,  and 
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Figure  7.2  Plan  view  of  the  Goodwin  Creek  Bendway  field  study  site. The  surveyed 
cross  sections  are  shaded  red. 


EASTING,  in  meter 

2 by  sapping  and  small  pop-out  failures  in  the  region  of  contrasting 

permeabilities  in  the  Holocene  units  about  1.6  to  2 m below  the  top  bank. 

Both  processes  resulted  in  oversteepening  at  the  base  of  the  EH  unit  and 
subsequent  collapse  during  wet  periods. 


Application 

We  employed  CONCEPTS  to  simulate  the  above  failure  processes  between 
March  1,  1996  and  March  31,  1997.  The  modeling  reach  contains  10  surveyed 
cross  sections  (see  Figure  7.2)  and  one  inserted  cross  section  20  m downstream 
of  cross  section  10.  This  cross  section  is  a copy  of  cross  section  10,  and  was 
inserted  to  prevent  the  downstream  boundary  condition  from  affecting  the 
simulated  results.  Observed  bed  elevations  did  not  vary  much  during  this 
period,  so  we  turned  off  the  sediment  transport  and  bed  adjustment  module  to 
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fix  the  bed.  We  imposed  the  combined  discharges  measured  at  flumes  3 and  4 
at  the  upstream  boundary  of  the  modeling  reach. 

Because  CONCEPTS  assumes  that  bank-material  properties  are  homogeneous, 
we  used  the  following  average  shear-strength  parameters:  c'  = 4.5  kPa, 

(f>'  = 28.3° , and  = 13.5° . The  phreatic  surface  was  located  3.5  m below 
the  top  of  the  bank.  CONCEPTS  follows  the  excess  shear-stress  method  of 
Osman  and  Thorne  (1988),  see  the  section  on  “Fluvial  Erosion  Process”  in 
Chapter  3,  to  calculate  the  lateral  erosion  at  the  bank  toe.  Critical  shear  stresses 
varied  along  the  right  bank  from  8 Pa  at  cross  section  1 to  1 .5  Pa  at  the  apex  of 
the  bend  (cross  section  6)  to  4 Pa  at  cross  section  10.  Because  CONCEPTS  is 
only  applicable  to  straight  channels,  we  have  to  vary  the  critical  shear  stresses 
to  simulate  channel  widening  due  to  bend  migration. 

Figure  7.3  shows  the  Run  Control  Data  File,  Figure  7.4  shows  an  example  of  a 
Cross  Section  Data  File,  and  Figure  7.5  shows  part  of  the  Discharge  Data  File. 

Figure  7.3  Run  control  data  file  used  in  the  simulation  of  the  Goodwin  Creek 
Bendway  field  study  site. 


! Main  Input  File 

i 

! case  name 
GCB 

! project  title 

Goodwin  Creek  Bendway  1 bank  stability  analysis 

! Run  Control  Data 

! upstream  flow  discharge  file 
GCBHydrography . txt 
! flow-related  data 
0.10  0.0  0 

! sediment  discharge  at  upstream  end  of  the  channel 
0 0.10  0.10  0.10  0.05  0.00  0.00  0.00  0.00  0.00  0.00  0.00  0.00  0.00 

! silt  fraction  and  downstream  bed  control 
0.30  1.0 

! bank  failure  analysis 
7 10 

! type  of  flow  resistance  formulation 
1 

! water  temperature 
10.0 

! sediment  and  streambank  mechanics  options 
0 1 1 

! Simulation  Times 

! start  end  time  step 

03/03/1996  00:00:00  03/31/1997  24:00:00  100 

i Makeup  of  Modeling  Reach  

! number  of  links 
1 

! linktypes  for  the  above  number  of  links 
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Figure  7.3  Run  control  data  file  used  in  the  simulation  of  the  Goodwin  Creek 
Bendway  field  study  site,  (continued) 


1 

! Link  1 

! REACH:  number  of  cross  sections  and  their  data  filenames 
11 

CrossSectionOl . txt 
CrossSection02 . txt 
CrossSection03 . txt 
CrossSection04 . txt 
CrossSection05 . txt 
Cross Sect ionO 6 . txt 
CrossSection07 . txt 
CrossSection08 . txt 
CrossSection09 . txt 
Cross Sect ionlO . txt 
CrossSectionll . txt 

I Output 

! single  point  and  time 
6 

17536 
1 1 
14 

03/06/1996  00:00:00 
03/25/1996  11:45:00 
04/22/1996  22:30:00 
06/12/1996  08:00:00 
07/24/1996  23:15:00 
09/21/1996  05:00:00 
11/07/1996  12:45:00 
11/30/1996  10:15:00 
12/16/1996  18:00:00 
12/26/1996  18:00:00 
01/24/1997  00:30:00 
02/04/1997  03:15:00 
03/02/1997  04:00:00 
03/18/1997  22:45:00 
17536 
1 3 

14 

03/06/1996  00:00:00 
03/25/1996  11:45:00 
04/22/1996  22:30:00 
06/12/1996  08:00:00 
07/24/1996  23:15:00 
09/21/1996  05:00:00 
11/07/1996  12:45:00 
11/30/1996  10:15:00 
12/16/1996  18:00:00 
12/26/1996  18:00:00 
01/24/1997  00:30:00 
02/04/1997  03:15:00 
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Figure  7.3  Run  control  data  file  used  in  the  simulation  of  the  Goodwin  Creek 
Bendway  field  study  site,  (continued) 


03/02/1997 

04:00:00 

03/18/1997 

22:45:00 

17536 

1 5 

14 

03/06/1996 

00:00:00 

03/25/1996 

11:45:00 

04/22/1996 

22:30:00 

06/12/1996 

08:00:00 

07/24/1996 

23:15:00 

09/21/1996 

05:00:00 

11/07/1996 

12:45:00 

11/30/1996 

10:15:00 

12/16/1996 

18:00:00 

12/26/1996 

18:00:00 

01/24/1997 

00:30:00 

02/04/1997 

03:15:00 

03/02/1997 

04:00:00 

03/18/1997 

22:45:00 

17536 

1 6 

14 

03/06/1996 

00:00:00 

03/25/1996 

11:45:00 

04/22/1996 

22:30:00 

06/12/1996 

08:00:00 

07/24/1996 

23:15:00 

09/21/1996 

05:00:00 

11/07/1996 

12:45:00 

11/30/1996 

10:15:00 

12/16/1996 

18:00:00 

12/26/1996 

18:00:00 

01/24/1997 

00:30:00 

02/04/1997 

03:15:00 

03/02/1997 

04:00:00 

03/18/1997 

22:45:00 

17536 

1 8 

14 

03/06/1996 

00:00:00 

03/25/1996 

11:45:00 

04/22/1996 

22:30:00 

06/12/1996 

08:00:00 

07/24/1996 

23:15:00 

09/21/1996 

05:00:00 

11/07/1996 

12:45:00 

11/30/1996 

10:15:00 

12/16/1996 

18:00:00 

12/26/1996 

18:00:00 

01/24/1997 

00:30:00 
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Figure  7.3  Run  control  data  file  used  in  the  simulation  of  the  Goodwin  Creek 
Bendway  field  study  site,  (continued) 


02/04/1997  03:15:00 
03/02/1997  04:00:00 
03/18/1997  22:45:00 
17536 
1 10 
14 

03/06/1996  00:00:00 
03/25/1996  11:45:00 
04/22/1996  22:30:00 
06/12/1996  08:00:00 
07/24/1996  23:15:00 
09/21/1996  05:00:00 
11/07/1996  12:45:00 
11/30/1996  10:15:00 
12/16/1996  18:00:00 
12/26/1996  18:00:00 
01/24/1997  00:30:00 
02/04/1997  03:15:00 
03/02/1997  04:00:00 
03/18/1997  22:45:00 


! single  point,  continuously  in  time 


6 

602125 
1 1 
1 

03/01/1996  00:00:00 
602125 
1 3 

1 

03/01/1996  00:00:00 
602125 
1 5 

1 

03/01/1996  00:00:00 
602125 
1 6 
1 

03/01/1996  00:00:00 
602125 
1 8 
1 

03/01/1996  00:00:00 
602125 
1 10 
1 

03/01/1996  00:00:00 
! profile  at  specific 
1 
7 

1 1 1 10 


03/31/1997  23:59:59 


03/31/1997  23:59:59 


03/31/1997  23:59:59 


03/31/1997  23:59:59 


03/31/1997  23:59:59 


03/31/1997  23:59:59 
time  points 
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Figure  7.3  Run  control  data  file  used  in  the  simulation  of  the  Goodwin  Creek 
Bendway  field  study  site,  (continued) 


14 

03/06/1996 

00 

00 

00 

03/25/1996 

11 

45 

00 

04/22/1996 

22 

30 

00 

06/12/1996 

08 

00 

00 

07/24/1996 

23 

15 

00 

09/21/1996 

05 

00 

00 

11/07/1996 

12 

45 

00 

11/30/1996 

10 

15 

00 

12/16/1996 

18 

00 

00 

12/26/1996 

18 

00 

00 

01/24/1997 

00 

30 

00 

02/04/1997 

03 

15 

00 

03/02/1997 

04 

00 

00 

03/18/1997 

22 

45 

00 

Figure  7.4  A sample  cross  section  data  file  used  in  the  simulation  of  the  Goodwin 
Creek  Bendway  field  study  site. 


1 

! Input  file  of  cross  section  5. 

i 

! Name  of  xsection  and  rivermile  in  (km) 

Cross  Section  5 Goodwin  Creek  Bendway  1 

0.0337 

! friction 

:actor 

0.06 

! discharge 

fraction  due  to  tributary  inflow 

0.00 

1 

Left  FloodPlain 

! number  of 

nodes 

5 

! station  and  elevation  for  above  number 

of  coordinates  in 

(m) 

-100.00 

85.37 

-100.00 

84.37 

0.00 

84.37 

0.54 

84.36 

1.97 

84.11 

! Manning's 

n of  left  Floodplain 

0.10 

I 

Left  Bank  

! number  of 

nodes 

6 

! station  and  elevation  for  above  number 

of  coordinates  in 

(m) 

1.97 

84.11 

3.39 

83.17 

7.58 

82.65 

11.10 

81.83 

13.24 

81.41 

13.77 

81.17 
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Figure  7.4  A sample  cross  section  data  file  used  in  the  simulation  of  the  Goodwin 
Creek  Bendway  field  study  site,  (continued) 

! bank  material 

4500.0  28.5  10.4  19400.0 

! bank  sediment  composition 
28.540 
18.240 
31.420 
21.370 
0.330 
0.040 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 

! critical  shear  stress  for  erosion 
7.7 

! groundwatertable,  hydraulic  gradient,  and  angle  of  seepage  force 
(degrees) 

81.00  1.00  0.00 
! Manning's  n 
0.04 

! Channel  Bed 

! number  of  nodes 
13 

! station  and  elevation  for  above  number  of  coordinates  in  (m) 


13. 

.77 

81. 

. 17 

15. 

.32 

81. 

.05 

18. 

.19 

80. 

.93 

18. 

.70 

80. 

.76 

20. 

.79 

80. 

.54 

21. 

.85 

80. 

.11 

22. 

. 14 

80. 

.04 

23. 

.39 

79. 

.93 

24, 

. 61 

79. 

.71 

25. 

.99 

79. 

.60 

27. 

.03 

79. 

.69 

27. 

.66 

79. 

.94 

27. 

.94 

80. 

.10 

! Elevation  of  bedrock  (m) 

0.00 

! porosity 
0.4 

! hiding  factor 
0.25  0.95  0.70 

! Surface  layer  and  Substratum  data 
! number  of  sediment  layers  composing  the  bed 
1 

! Layer  1,  layer  depth  below  bed  surface 
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Figure  7.4  A sample  cross  section  data  file  used  in  the  simulation  of  the  Goodwin 
Creek  Bendway  field  study  site,  (continued) 

0.00 

! bed  composition 
21.860 
11.390 
28.360 
37.690 
0.520 
0.060 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 

! critical  shear  stresses  for  erosion  of  and  deposition  on  cohesive 
beds, 

! and  erodibility  coefficient 
0.10  7.70  2.18E-07 

! Manning's  n 
0.03 

I Right  Bank 

! number  of  nodes 
3 

! station  and  elevation  for  above  number  of  coordinates  in  (m) 

27.94  80.10 

28.17  81.16 

28.90  84.39 

! bank  material 

8500.0  28.5  13.5  19400.0 

! bank  sediment  composition 
28.540 
18.240 
31.420 
21.370 
0.330 
0.040 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 

! critical  shear  stress  for  erosion 
2.3 

! groundwatertable,  hydraulic  gradient,  and  angle  of  seepage  force 
(degrees ) 

81.00  1.00  0.00 
! Manning's  n 
0.03 
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Figure  7.4  A sample  cross  section  data  file  used  in  the  simulation  of  the  Goodwin 
Creek  Bendway  field  study  site,  (continued) 

; Right  FloodPlain 

! number  of  nodes 
5 

! station  and  elevation  for  above  number  of  coordinates  in  (m) 

28.90  84.39 

31.11  84.52 

31.20  84.67 

131.20  84.67 

131.20  85.67 

! Manning's  n of  right  Floodplain 
0.05 


Figure  7.5  Part  of  the  discharge  data  file  used  in  the  simulation  of  the  Goodwin  Creek 
Bendway  field  study  site. 


03/20/1996 

11 

59 

00 

0 

0 

13 

03/20/1996 

21 

20 

00 

2 

0 

10 

03/24/1996 

23 

39 

00 

1 

0 

10 

03/24/1996 

23 

41 

00 

0 

0 

11 

03/24/1996 

23 

44 

00 

0 

0 

13 

03/24/1996 

23 

45 

00 

0 

0 

14 

03/24/1996 

23 

46 

00 

0 

0 

19 

03/24/1996 

23 

48 

00 

0 

0 

28 

03/24/1996 

23 

51 

00 

0 

0 

56 

03/24/1996 

23 

52 

00 

0 

0 

74 

03/25/1996 

00 

00 

00 

0 

2 

15 

03/25/1996 

00 

01 

00 

0 

2 

42 

03/25/1996 

00 

05 

00 

0 

3 

33 

03/25/1996 

00 

08 

00 

0 

3 

93 

03/25/1996 

00 

10 

00 

0 

4 

52 

03/25/1996 

00 

14 

00 

0 

5 

52 

03/25/1996 

00 

16 

00 

0 

6 

09 

03/25/1996 

00 

19 

00 

0 

6 

88 

03/25/1996 

00 

20 

00 

0 

7 

39 

03/25/1996 

00 

25 

00 

0 

9 

33 

03/25/1996 

00 

27 

00 

0 

10 

27 

03/25/1996 

00 

28 

00 

0 

10 

80 

03/25/1996 

00 

30 

00 

0 

11 

83 

03/25/1996 

00 

37 

00 

0 

14 

25 

03/25/1996 

00 

38 

00 

0 

14 

60 

03/25/1996 

00 

47 

00 

0 

17 

33 

03/25/1996 

00 

49 

00 

0 

17 

93 

03/25/1996 

00 

50 

00 

0 

18 

11 

03/25/1996 

00 

59 

00 

0 

20 

11 

03/25/1996 

01 

01 

00 

0 

20 

60 

03/25/1996 

01 

02 

00 

0 

20 

80 

03/25/1996 

01 

07 

00 

0 

21 

89 

03/25/1996 

01 

11 

00 

0 

22 

77 

03/25/1996 

01 

12 

00 

0 

22 

94 

03/25/1996 

01 

14 

00 

0 

23 

26 
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Figure  7.5  Part  of  the  discharge  data  file  used  in  the  simulation  of  the  Goodwin  Creek 
Bendway  field  study  site,  (continued) 


03/25/1996 

01 

15 

00 

0 

23 

41 

03/25/1996 

01 

19 

00 

0 

23 

92 

03/25/1996 

01 

24 

00 

0 

24 

35 

03/25/1996 

01 

25 

00 

0 

24 

46 

03/25/1996 

01 

33 

00 

0 

24 

98 

03/25/1996 

01 

35 

00 

0 

25 

02 

03/25/1996 

01 

36 

00 

0 

24 

98 

03/25/1996 

01 

42 

00 

0 

24 

99 

03/25/1996 

01 

43 

00 

0 

24 

99 
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01 

47 
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0 

24 

90 

03/25/1996 

01 

51 

00 

0 

24 

80 

03/25/1996 

01 

55 

00 

0 

24 

53 

03/25/1996 

01 

58 

00 

0 

24 

35 

03/25/1996 

02 

01 

00 

0 

24 

21 

03/25/1996 

02 

04 

00 

0 

24 

03 

03/25/1996 

02 

09 

00 

0 

23 

75 

03/25/1996 

02 

10 

00 

0 

23 

71 

03/25/1996 

02 

13 

00 

0 

23 

42 

03/25/1996 

02 

19 

00 

0 

22 

93 

03/25/1996 

02 

21 

00 

0 

22 

76 

03/25/1996 

02 

31 

00 

0 

21 

96 

03/25/1996 

02 

43 

00 

0 

20 

94 

03/25/1996 

02 

52 

00 

0 

20 

31 

03/25/1996 

02 

57 

00 

0 

20 

04 

03/25/1996 

03 

11 

00 

0 

19 

25 

03/25/1996 

03 

21 

00 

0 

18 

64 

03/25/1996 

03 

31 

00 

0 

17 

87 

03/25/1996 

03 

33 

00 

0 

17 

70 

03/25/1996 

03 

48 

00 

0 

16 

48 

03/25/1996 

04 

18 

00 

0 

13 

75 

03/25/1996 

04 

39 

00 

0 

11 

42 

03/25/1996 

04 

41 

00 

0 

11 

20 

03/25/1996 

04 

47 

00 

0 

10 

50 

03/25/1996 

04 

51 

00 

0 

10 

06 

03/25/1996 

04 

57 

00 

0 

9 

53 

03/25/1996 

05 

02 

00 

0 

8 

97 

03/25/1996 

05 

10 

00 

0 

8 

34 

03/25/1996 

05 

15 

00 

0 

7 

92 

03/25/1996 

05 

24 

00 

0 

7 

24 

03/25/1996 

05 

28 

00 

0 

6 

94 

03/25/1996 

05 

36 

00 

0 

6 

46 

03/25/1996 

05 

46 

00 

0 

5 

87 

03/25/1996 

05 

49 

00 

0 

5 

68 

03/25/1996 

05 

56 

00 

0 

5 

32 

03/25/1996 

06 

04 

00 

0 

5 

04 

03/25/1996 

06 

11 

00 

0 

4 

81 

03/25/1996 

06 

18 

00 

0 

4 

56 
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06 

30 

00 

0 

4 

19 

03/25/1996 

06 

40 

00 

0 

3 

87 

03/25/1996 

06 

45 

00 

0 

3 

73 
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Figure  7.5  Part  of  the  discharge  data  file  used  in  the  simulation  of  the  Goodwin  Creek 
Bendway  field  study  site,  (continued) 


03/25/1996 

06 

56: 

00 

0 

3. 

49 

03/25/1996 

07 

11: 

00 

0 

3. 

16 

03/25/1996 

07 

13. 

00 

0 

3 

12 

03/25/1996 

07 

27- 

00 

0 

2 

86 
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07 

42 
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2 

58 
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07 
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2 

49 
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07 

57 
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0 

2 

47 

03/25/1996 

08 

28 

00 

0 

2 

32 

03/25/1996 

08 

50 

00 

0 

2 

26 

03/25/1996 

09 

04 

00 

0 

2 

20 

03/25/1996 

09 

34 

00 

0 

2 

08 

03/25/1996 

09 

59 

00 

0 

2 

01 

03/25/1996 

10 

48 

00 

0 

1 

87 
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11 

41 
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0 

1 

76 
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12 

17 

00 

0 

1 

68 
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49 
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34 
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29 
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09 
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56 
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87 
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62 
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56 
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0 

0 

53 
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01 
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0 

0 

53 
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00 
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48 
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0 

0 

48 
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09 

13 
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0 

0 

39 
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09 

18 

00 

0 

0 

39 
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00 

00 

00 

0 

0 

25 
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00 

01 

00 

0 

0 

23 
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13 

05 
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0 

0 

17 
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00 

00 

0 

0 

12 
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12 

20 

00 

2 

0 

10 
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09 
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0 

10 
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09 

21 
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Figure  7.6  compares  simulated  and  surveyed  cross  sections  6 and  8.  The  right 
bank  of  cross  section  6 underwent  two  successive  failures  on  December  1, 
1996  and  midwinter  1997  (Figure  7.6a).  The  retreat  of  the  top  bank  was  0.6  m 
and  1 m,  respectively.  CONCEPTS  simulated  failure  of  the  bank  on  February 
4,  1997,  with  the  top  bank  retreating  2.7  m (Figure  7.6a).  CONCEPTS 
underpredicted  the  rate  of  basal  scour.  Also,  the  observed  failure  plane  angle  is 
greater  than  that  computed.  Cantilever  failures  will  generally  occur  along 
steeper  slip  surfaces  than  planar  failures.  The  cantilever  failure  is  not  yet 
implemented  in  CONCEPTS. 

Figure  7.6b  shows  that  the  right  bank  of  cross  section  8 fails  in  late  April,  1996 
(retreat  of  the  top  bank  is  1 .0  m)  and  again  on  February  3,  1997  (retreat  of  top 
bank  is  1.6  m).  CONCEPTS  simulated  failure  of  the  bank  on  March  6,  1997. 
The  slip  surface  intersects  the  bank  3.1m  below  the  top  bank.  CONCEPTS 
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Figure  7.6  Comparison  of  observed  and  simulated  bank  failures  at  the  Goodwin 

Creek  Bendway  field  study  site:  (a)  cross  section  6 and  (b)  cross  section  8. 
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was  unable  to  predict  the  ‘February  1997’  failure  because  the  flow  did  not 
remove  all  of  the  slumped  material  from  the  ‘March  1996’  failure.  On  March 

31,  1997,  0.05  m^/m  of  the  modeled,  slumped  material  remained.  In 
CONCEPTS  the  next  simulated  flow  event  after  March  31,  1997  should  be 
able  to  erode  the  remaining  slumped  material  and  continue  to  steepen  the  bank, 
which  may  then  fail  in  the  spring  of  1997. 
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Figure  7.7  Predicted  factor  of  safety  for  the  right  banks  of  cross  sections  6 and  8. 
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Figure  7.7  shows  how  the  factor  of  safety  of  the  right  bank  of  cross  section  6 
reduces  after  each  runoff  event.  Scour  of  bank  material  near  the  bank  toe 
steepens  the  bank  and  reduces  the  factor  of  safety.  The  spikes  in  the  factor  of 
safety  represent  the  increase  in  factor  of  safety  due  to  rising  flow  stage  and 
thus  an  increase  in  confining  pressure. 
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Little  Salt  Creek,  Nebraska 

This  chapter  presents  the  application  of  CONCEPTS  to  the  middle 
reach  of  Little  Salt  Creek  in  eastern  Nebraska.  You  will  find  (1)  the 
physical  characteristics  of  the  channel,  (2)  the  input  data  necessary 
to  perform  the  simulation,  and  (3)  results  showing  the  capability  of 
CONCEPTS  to  simulate  the  morphology  of  a channel  made  up  of 
cohesive  bed  and  bank  material  and  the  effects  of  grade  control 
structures  on  channel  stability. 


Study  Area 

The  Little  Salt  Creek  is  a tributary  of  Salt  Creek  and  drains  approximately  1 19 

km  (see  Figure  8.1).  Bed  and  bank  materials  are  cohesive.  The  average 
channel  slope  is  about  0.0015.  Saline  wetlands  in  the  Salt  Creek  basin  have 
diminished  from  1 6,000  acres  to  approximately  1 ,200  acres.  They  hold  a 
diversity  of  wildlife  and  contain  rare  plant  communities.  A saline  wetland 
north  of  raymond  Road  (Figure  8.2)  is  endangered  by  channel  incision,  which 
may  lead  to  draining  of  the  wetland.  Channel  incision  may  also  lead  to 
additional  streambank  failures  and  loss  of  wetland  area  (Figure  8.3). 

Scientists  of  the  US  Department  of  Agriculture- Agricultural  Research  Service- 
National  Sedimentation  Laboratory  (NSL)  have  collected  data  in  cooperation 
with  the  US  Geological  Survey  (USGS)  on  streambank  shear- strength 
parameters  and  streambed  erodibility.  In  addition,  the  USGS  surveyed  the 
study  reach  (Raymond  Road  to  Bluff  Road,  see  Figure  8.1),  took  bed  and  bank- 
material  samples,  and  collected  discharge  data. 
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Figure  8.1  Map  of  Little  Salt  Creek,  Lancaster  County,  Nebraska. 
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Study  Area 


Figure  8.2  View  of  Little  Salt  Creek  and  wetlands  upstream  of  Raymond  Road. 


Channel  Geometry 


The  USGS  surveyed  seven  cross  sections  upstream  of  Raymond  Road,  six 
cross  sections  upstream  of  Mill  Road,  and  seven  cross  sections  near  Bluff 
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Road.  There  are  two  reaches  that  could  not  be  surveyed  because  they  were 

inaccessible: 

1 A reach  extending  from  raymond  Road  to  the  first  surveyed  cross  section 
upstream  of  Mill  Road.  The  length  of  this  reach  is  1 .7  km.  Seven  simulated 
cross  sections  were  inserted  along  this  reach  with  equidistant  spacing.  Their 
geometries  gradually  varied  from  that  of  the  most  downstream  surveyed 
cross  section  at  Raymond  Road  to  that  of  the  most  upstream  surveyed  cross 
section  at  Mill  Road. 

2 A reach  extending  from  Mill  Road  to  the  most  upstream  surveyed  cross 
section  near  Bluff  Road.  The  length  of  this  reach  is  3.6  km.  Seventeen 
simulated  cross  sections  were  similarly  inserted  along  this  reach  between 
Mill  and  Bluff  Roads. 

Hydrographs 


A USGS  gaging  station  is  located  at  Arbor  Road  on  Little  Salt  Creek.  The 
USGS  provided: 

• daily  discharges  for  water  years  1969  through  1998, 

■j 

• peak-flow  data  above  a base  discharge  of  1 5.6  m /s  for  the  same  time 
period,  and 

• half-hourly  instantaneous  discharges  for  water  years  1991-1 994  and  1 996- 
1998. 

Mean-daily  discharges  may  not,  however,  be  a good  representative  of  actual 
flow  because  they  cannot  account  for  the  peak  discharge  and  the  commonly 
rapid  rise  from  base  flow  to  peak  flow.  The  largest  shear  stresses  exerted  by  the 
flow  on  the  streambed  often  occur  near  peak  flow.  It  is  important  to  properly 
calculate  the  bed  shear  stress  to  accurately  predict  the  evolution  of  the 
streambed.  The  following  procedure  was  used  to  convert  daily  discharges  to 
unsteady  flow  hydrographs: 

1 We  selected  a base-flow  discharge  of  0. 1 2 m /s  after  analysis  of  the  daily 
discharge  values. 

2 Using  the  available  peak-flow  data,  we  determined  a ratio  between  peak 
discharge  and  its  corresponding  daily  discharge  for  each  of  the  90  peak 
discharges.  The  average  of  all  ratios  was  3.9. 

3 We  traced  runoff  events  from  the  daily  discharges.  We  set  discharges 
smaller  than  the  base-flow  discharge  equal  to  the  base-flow  discharge.  We 
then  multiplied  the  maximum  daily  discharge  in  a runoff  event  by  3.9  to 
obtain  peak  discharge.  We  replaced  calculated  peak  discharge  by  that 
observed  for  known  peak  discharges.  This  resulted  in  breakpoint  data  for 
each  runoff  event.  We  assumed  that  each  discharge  value  occurred  at  noon 
of  that  day.  We  added  two  breakpoints  to  the  day  of  peak  flow  to  preserve 
runoff  volume. 
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CONCEPTS  requires  discharge  hydrographs  at  the  upstream  boundary  of  the 
modeling  reach.  Using  drainage-area  analysis,  we  converted  discharges  at  the 
gaging  station  to  those  at  the  upstream  boundary  near  Raymond  Road. 
Similarly,  we  can  determine  discharges  at  the  downstream  boundary  near  Bluff 
Road.  Hence,  we  can  calculate  the  coefficient  Cq  in  (4. 1 ) as 


u = z 


DA 


Bluff 


DA 


Raymon 


- 1 


(8.1) 


where  L is  channel  length  between  Raymond  and  Bluff  Roads,  DA  is  drainage 
area,  and  the  exponent  \ varies  between  0.7  and  1.0  (Leopold,  1994).  Here, 

CQ  = 9.07  x 10“5  using  £U*Raymond  = 43.8  km2,  £MBluff  = 85.8  km2, 

= 0.8,  and  L = 7.86  km. 


Bed-Material  Properties 

We  performed  in  situ  field  experiments  using  a submersible  jet  device 
(Hanson,  1990)  to  measure  the  erodibility  coefficient  and  critical  shear  stress 
of  the  cohesive  streambed  material  (cf.  (2.7)).  The  apparatus  consists  of  a 
submerged  jet  with  a nozzle  diameter  of  13  mm,  set  at  a height  of  0.22  m above 
the  initial  soil  surface.  Monitoring  of  the  depth  of  scour  during  a test  yields  the 
erodibility  coefficient  and  critical  shear  stress.  The  USGS  collected  bed- 
material  samples  to  determine  particle-size  distributions. 

Critical  shear  stress  is  fairly  constant  along  the  modeling  reach  (Langendoen 
and  Simon,  2000).  The  average  critical  shear  stress  is  7.7  Pa.  The  average 

erodibility  (defined  here  as  e/xc ) is  0.28x  1 O'6  m/Pa-s.  The  streambed  material 
is  a silt  loam  with  a clay  content  of  approximately  14%. 

Bank-Material  Properties 

We  performed  a series  of  in  situ  field  experiments  using  a borehole  shear  test 
(BST)  device  (Luttenegger  and  Hallberg,  1 98 1 ) to  determine  the  shear  strength 
of  the  cohesive  streambank  material.  The  BST  determines  consolidated, 
drained  apparent  cohesion  and  effective  angle  of  internal  friction.  Using  the 
known  elevation  of  the  phreatic  surface,  we  calculated  average  effective 
cohesion  as  4.3  kPa  and  average  angle  of  internal  friction  as  29°. 

The  USGS  collected  bank-material  samples  at  the  BST  locations  to  determine 
bank-material  composition  and  unit  weight.  The  average  bulk  density  is 

1.58  g/cm3. 
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Application 

We  employed  CONCEPTS  to  study  the  long-term  stability  of  Little  Salt  Creek 
between  Raymond  and  Bluff  Roads  using  the  30-year  discharge  breakpoint 
data  between  1 969  and  1 998.  The  modeling  reach  contains  20  surveyed  cross 
sections  and  26  simulated  cross  sections.  We  inserted  two  of  the  simulated 
cross  sections  between  the  third  and  fourth  surveyed  cross  section  near  Bluff 
Road  to  represent  a drop  in  the  streambed. 

Data  on  the  critical  shear  stress  to  entrain  bank-material  particles  was 
unavailable.  We  assumed  it  to  equal  that  of  bed-material  particles,  that  is 

7.7  Pa.  We  assumed  the  bulk  unit  weight  of  the  bank  material  to  be  17  kN/m3, 
which  is  a 1 0%  increase  over  the  average  dry  bulk  unit  weight  to  include  water 

content  of  the  bank  material.  The  angle  (f)^  indicating  the  increase  in  shear 
strength  for  an  increase  in  matric  suction  was  assumed  to  be  17°.  The 
groundwater  table  at  each  streambank  was  1 m above  the  elevation  of  the  bank 
toe. 

Hereafter,  we  refer  to  the  distance  downstream  of  the  upstream  boundary  of  the 
modeling  reach  as  model  kilometer. 

Figure  8.4  shows  the  Run  Control  Data  File , Figure  8.5  shows  an  example  of  a 
Cross  Section  Data  File , and  Figure  8.6  shows  part  of  the  Discharge  Data  File. 


Figure  8.4  Run  control  data  file  used  in  the  simulation  of  the  long  term  stability  of 
Little  Salt  Creek,  Nebraska. 


Main  Input  File 


! case  name 
Salt 

! project  title 

Little  Salt  Creek  channel  evolution  simulation 

i Run  control  data 

! upstream  flow  discharge  file 
LSChydrography . txt 
! flow-related  data 
0.12  9.07E-5  0 

! sediment  discharge  at  upstream  end  of  the  channel 
0 0.10  0.10  0.10  0.05  0.00  0.00  0.00  0.00  0.00  0.00  0.00  0.00  0.00 

! silt  fraction 
0.30 

! bank  failure  analysis 
0 3 

! type  of  flow  resistance  formulation 
1 

! water  temperature 

10.0 


start 


Simulation  Times  

end  time  step 
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Figure  8.4  Run  control  data  file  used  in  the  simulation  of  the  long  term  stability  of 
Little  Salt  Creek,  Nebraska,  (continued) 

02/25/1969  00:00:00  09/30/1998  23:59:00  100 

! Makeup  of  Modeling  Reach 

! number  of  links 
1 

! linktypes  for  the  above  number  of  links 
1 

i Link  1 

! REACH:  number  of  cross  sections  and  their  data  filenames 
46 

Ray-xs-1 . txt 
Ray-xs-2 . txt 
Ray-xs-3 . txt 
Ray-xs-4 . txt 
Ray-xs-5 . txt 
Ray-xs-6 . txt 
Ray-xs-7 . txt 
RM-xs-1 . txt 
RM-xs-2 . txt 
RM-xs-3 . txt 
RM-xs-4 . txt 
RM-xs-5 . txt 
RM-xs-6 . txt 
RM-xs-7 . txt 
Mil-xs-1 . txt 
Mil-xs-2 . txt 
Mil-xs-3 . txt 
Mil-xs-4 . txt 
Mil-xs-5 . txt 
Mil-xs-6 . txt 
MW-xs-1 . txt 
MW-xs-2 . txt 
MW-xs-3 . txt 
MW-xs-4 . txt 
MW-xs-5 . txt 
MW-xs-6 . txt 
MW-xs-7 . txt 
MW-xs-8 . txt 
MW-xs-9 . txt 
MW-xs-10 . txt 
MW-xs-11 . txt 
WB-xs-1 . txt 
Wav-xs-1 . txt 
WB-xs-2 . txt 
WB-xs-3 . txt 
WB-xs-4 . txt 
WB-xs-5 . txt 
Blu-xs-1 . txt 
Blu-xs-2 . txt 
Blu-xs-3 . txt 
Blu-xs-3B. txt 
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Figure  8.4  Run  control  data  file  used  in  the  simulation  of  the  long  term  stability  of 
Little  Salt  Creek,  Nebraska,  (continued) 


Blu-xs-3C.txt 
Blu-xs-4 . txt 
Blu-xs-5 . txt 
Blu-xs-6 . txt 
Blu-xs-7 . txt 

i 

Output  

single  point  and  time 
3 

17536 
1 7 

5 

03/30/1976  06:00:00 
11/10/1983  06:00:00 
01/22/1989  06:00:00 
12/22/1993  22:00:00 


06/20/1998  17 
17536 
1 19 

5 

O 

o 

o 

o 

03/30/1976  06:00:00 
11/10/1983  06:00:00 
01/22/1989  06:00:00 


12/22/1993  22 
06/20/1998  17 
17536 
1 39 

5 

o o 
o o 

o o 
o o 

03/30/1976  06:00:00 
11/10/1983  06:00:00 
01/22/1989  06:00:00 


12/22/1993  22 
06/20/1998  17 
! single  point, 
19 

602125 
1 1 
1 

00:00 

00:00 

continuously  in  time 

01/01/1969  00 
602125 
1 7 

1 

00:00  09/30/1998  23:59:59 

01/01/1969  00 
602125 
1 12 
1 

00:00  09/30/1998  23:59:59 

01/01/1969  00 
602125 
1 17 

1 

00:00  09/30/1998  23:59:59 

01/01/1969  00 

00:00  09/30/1998  23:59:59 
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Figure  8.4  Run  control  data  file  used  in  the  simulation  of  the  long  term  stability  of 
Little  Salt  Creek,  Nebraska,  (continued) 


602125 
1 22 
1 

01/01/1969 
602125 
1 28 
1 

01/01/1969 
602125 
1 33 

1 

01/01/1969 
602125 
1 39 

1 

01/01/1969 
602125 
1 46 

1 

01/01/1969 


00:00:00 


00:00:00 


00:00:00 


00:00:00 


00:00:00 


09/30/1998  23:59:59 


09/30/1998  23:59:59 


09/30/1998  23:59:59 


09/30/1998  23:59:59 


09/30/1998  23:59:59 


! profile  at  specific  time  points 
1 
7 

1 1 1 46 

5 

03/30/1976  06:00:00 
11/10/1983  06:00:00 
01/22/1989  06:00:00 
12/22/1993  22:00:00 
06/20/1998  17:00:00 


Figure  8.5  A sample  cross  section  data  file  used  in  the  simulation  of  the  long-term 
stability  of  Little  Salt  Creek,  Nebraska. 


! Input  file  of  cross  section  19. 

i 

! Name  of  xsection  and  rivermile  in  (km) 

Cross  Section  5 U/S  Mill  Road 
2.9263 

! friction  factor 
0.06 

! discharge  fraction  due  to  tributary  inflow 
0.00 

i Left  FloodPlain 

! number  of  nodes 
4 

! station  and  elevation  for  above  number  of  coordinates  in  (m) 
-69.33  355.85 

-69.33  354.85 
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Figure  8.5  A sample  cross  section  data  file  used  in  the  simulation  of  the  long-term 
stability  of  Little  Salt  Creek,  Nebraska,  (continued) 

-19.33  354.85 

-12.28  354.85 

! Manning's  n of  left  Floodplain 
0.15 

! Left  Bank 

! number  of  nodes 
6 

! station  and  elevation  for  above  number  of  coordinates  in  (m) 
-12.28  354.85 

-10.71  353.92 

-4.73  353.39 

-3.56  352.44 

-2.80  352.40 

-2.09  351.51 

! bank  material 

4500.0  30.0  17.0  17000.0 

! bank  sediment  composition 
28.540 
18.240 
31.420 
21.370 
0.330 
0.040 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 

! critical  shear  stress  for  erosion 
7.7 

! groundwatertable,  hydraulic  gradient,  and  angle  of  seepage  force 
(degrees) 

352.51  1.00  0.00 

! Manning's  n 
0.06 

i Channel  Bed 

! number  of  nodes 
3 

! station  and  elevation  for  above  number  of  coordinates  in  (m) 

-2.09  351.51 

0.00  351.33 

0.69  351.38 

! Elevation  of  bedrock  (m) 

0.00 

! porosity 
0.4 

! hiding  factor 
0.25  0.95  0.70 
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Figure  8.5  A sample  cross  section  data  file  used  in  the  simulation  of  the  long-term 
stability  of  Little  Salt  Creek,  Nebraska,  (continued) 

! Surface  layer  and  Substratum  data 
! number  of  sediment  layers  composing  the  bed 
1 

! Layer  1,  layer  depth  below  bed  surface 

0.00 

! bed  composition 
21.860 
11.390 
28.360 
37.690 
0.520 
0.060 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 

! critical  shear  stresses  for  erosion  of  and  deposition  on  cohesive 
beds, 

! and  erodibility  coefficient 
0.10  7.70  2.18E-07 

! Manning ' s n 
0.025 

I Right  Bank 

! number  of  nodes 
8 

! station  and  elevation  for  above  number  of  coordinates  in  (m) 


0. 

.69 

351. 

.38 

0. 

.88 

352. 

.04 

1. 

.47 

352. 

.18 

2. 

.59 

353. 

.00 

3. 

.39 

353. 

.08 

5. 

.19 

353. 

.55 

6. 

.47 

354. 

.37 

6. 

.58 

354. 

.86 

! bank  material 

4500.0  30.0  17.0  17000.0 

! bank  sediment  composition 
28.540 
18.240 
31.420 
21.370 
0.330 
0.040 
0.000 
0.000 
0.000 
0.000 
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Figure  8.5  A sample  cross  section  data  file  used  in  the  simulation  of  the  long-term 
stability  of  Little  Salt  Creek,  Nebraska,  (continued) 

0.000 

0.000 

0.000 

! critical  shear  stress  for  erosion 
7.7 

! groundwatertable,  hydraulic  gradient,  and  angle  of  seepage  force 
(degrees) 

352.38  1.00  0.00 

! Manning's  n 
0.06 

i Right  FloodPlain 

! number  of  nodes 
4 

! station  and  elevation  for  above  number  of  coordinates  in  (m) 

6.58  354.86 

15.08  354.94 

65.08  354.94 

65.08  355.94 

! Manning's  n of  right  Floodplain 
0.15 


Figure  8.6  Part  of  the  discharge  data  file  used  in  the  simulation  of  the  long-term 
stability  of  Little  Salt  Creek,  Nebraska. 
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Figure  8.7  Simulated  evolution  of  the  thalweg  profile,  Little  Salt  Creek,  for  a 30-year 
period. 


MODEL  KILOMETER,  in  kilometer 

Figure  8.7  shows  the  predicted  thalweg  profiles  of  the  modeling  reach  at 
various  points  of  time.  Figure  8.8  shows  the  simulated  changes  in  cross- 
sectional  geometry  at  cross  sections  upstream  of  Raymond  Road,  upstream  of 
Mill  Road,  and  near  Bluff  Road. 

The  erodibility  of  the  streambed  material  governs  the  morphology  of  the 
modeling  reach.  In  situ  field  experiments  suggest  that  erodibility  of  the 
streambed  material  is  similar  along  the  modeling  reach.  As  a result,  channel 
incision  increases  along  the  modeling  reach  because  of  increasing  discharge 
and,  therefore,  shear  stress.  Channel  incision  varies  from  about  4.3  m near  the 
downstream  end  of  the  channel  (model  kilometer  7. 1 ) to  1 .5  m at  the  upstream 
end  of  the  channel.  The  increase  in  bank  height  associated  with  incision  leads 
to  bank  failures  along  the  entire  modeling  reach.  The  average  widening  is  0.6 
m upstream  of  Raymond  Road,  1.5  m upstream  of  Mill  road,  0.4  m upstream  of 
North  14th  Street  and  Waverly  Road,  and  1.2  m near  Bluff  Road. 

Given  the  observed  erodibility  of  the  surface  of  the  streambed  and  types  of 
flows  occurring  between  1969  and  1998.  Little  Salt  Creek  will  incise  further. 
Figure  8.7  shows  that  incision  progresses  from  the  middle  and  downstream  end 
of  the  modeling  reach  to  its  upstream  end.  Hence,  controlling  channel  grade  at 
selected  locations  may  deter  incision. 

We  performed  a simulation  in  which  we  controlled  the  grade  of  all  bridge 
crossings  along  the  modeling  reach.  Figure  8.9  shows  the  resulting  simulated 
evolution  of  the  thalweg  profile.  With  the  addition  of  the  simulated  grade- 
control  structures  at  the  bridge  crossings  there  is  no  incision  upstream  of 
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Little  Salt  Creek,  Nebraska 


Figure  8.8  Computed  changes  in  cross  section  geometry  at  selected  cross  sections 
along  Little  Salt  Creek  for  a 30-year  period:  (a)  cross  section  at  model 
kilometer  0.29  upstream  of  Raymond,  (b)  cross  section  at  model  kilometer 
2.93  upstream  of  Mill  Road,  and  (c)  cross  section  at  model  kilometer  6.98 
near  Bluff  Road. 
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Figure  8.9  Simulated  evolution  of  the  thalweg  profile,  Little  Salt  Creek,  for  a 30-year 
period  and  the  streambed  controlled  at  each  bridge  crossing. 


MODEL  KILOMETER,  in  kilometer 

Raymond  Road  at  the  saline  wetlands.  The  simulated  structures  reduce  incision 
between  North  1st  Street  and  Bluff  Road  by  approximately  60  percent.  Figure 
8.9  also  shows  that  the  second  grade-control  structure,  bridge  crossing  at  North 
1 st  Street,  may  be  omitted. 
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